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Chapter 1 



Introduction 



Active galactic nuclei, or AGN, have received considerable attention during 
the last 40 years. When Maarten Schmidt recognized in 1963 that the 
quasi-stellar object 3C273 has a redshift of almost 0.16 (Schmidt 1963), it 
became immediately clear that there were objects outside our galaxy with 
tremendous luminosities. Zeldovich (1964), Salpeter (1964) and Lynden-Bell 
(1969) suggested that extragalactic radio sources were mainly driven by 
gas accreted into a disc around super-massive black holes, and Blandford 
& Konigl (1979) suggested that the radio emission from AGN is produced 
by a relativistic outflow of plasma along magnetic field lines. Their model, 
with numerous modifications, is still thought to be valid, and hence the basic 
foundations of what we know about AGN today are almost 25 years old. 

Like in all branches of astronomy, the progress in AGN science was 
tightly correlated with technical improvements. Radio astronomy started 
out with very low angular resolution due to relatively small dishes and long 
wavelengths. This made XjD, governing the resolution of any observing 
instrument, very poor compared to optical instruments. Almost simultane- 
ously, however, astronomers started to experiment with radio interferometers 
with progressively longer baselines and finally, to combine interferometric 
measurements to simulate a larger dish (Ryle & Hewish 1960). This technique 
is called "aperture synthesis" and was transferred from directly linked, locally 
distributed antennas to spatially widely separated radio telescopes in the 
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1.1. The AGN Zoo 



1970s. What is known as VLBI, or Very Long Baseline Interferometry, 
today, is a combination of single-dish radio astronomy, interferometry and 
aperture synthesis, and therefore certainly one of the most advanced technical 
acliievements of the 20th century. Consequently, in 1974, Sir Martin Ryle was 
awarded the Nobel prize in physics for his contribution. 

AGN science and VLBI are tightly interrelated: no other instrument yields 
the angular resolution necessary to spatially resolve the innermost regions of 
AGN and details in the jets, and AGN are almost the only objects that are 
observable with VLBI, with a few exceptions like masers or extremely hot 
gases. Thus, one cannot live without the other, but the connection has turned 
out to be very fruitful. 



1.1 The AGN Zoo 

The first AGN-related phenomena, although not recognized as such, were re- 
ported as early as 1908 (Fath 1908) and in the following decades (Shpher 1917, 
Humason 1932, Mayall 1934), when observers noticed bright emission lines 
from nuclei of several galaxies, and the first extragalactic jet in M87 was ob- 
served (Curtis 1918). The first systematic study of these objects was carried 
out by Carl Seyfert (Seyfert 1943). He noticed that emission lines from the nu- 
clei of six galaxies were unusually broadened (several lOOOkms"^), but sim- 
ply stated that this phenomenon was "probably correlated with the physical 
properties of the nucleus", without any further interpretation. In the 1950s and 
1960s, the first radio sky surveys were completed and catalogues, like the Third 
Cambridge Catalogue, were published. When interferometric techniques were 
further developed, all kinds of differences between the AGN turned up, and 
following ancient habits, astronomers started to classify what they observed. 
Meanwhile, they have created a colourful collection of mostly phenotypical 
classes and an equally colourful bunch of acronyms to describe these classes. 

1 .1 .1 The General Picture 

The current idea of what AGN are is shown in Fig. 1.1. A supermassive black 
hole (M, between 10^ and 10^ M0) accretes surrounding material that 
settles in a circumnuclear accretion disc. Some of the material's potential 
energy in the gravitational field of the AGN is turned into radiation by 
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Figure 1 .1 : Schematic of AGN (not to scale). Different AGN types are believed to be due to different 
viewing angles. 
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viscous friction in the accretion disc, but most of both matter and energy 
ends up in the black hole. Some of the material is expelled into two jets in 
opposite directions. Whether the jets are made of an electron-proton or an 
electron-positron plasma is still controversial, but as the observed radiation is 
certainly synchrotron radiation, given the high brightness temperatures and 
degrees of polarization, it must be ionized material circulating in magnetic 
fields. The jets are not smooth, and in many cases new components are 
observed as they are ejected from the AGN and travel outwards into the jet 
direction, and jet bends of any angle are observed. Though not very well 
determined, around 10 % of the material's rest mass is turned into energy, the 
process thus being incredibly efficient compared to hydrogen burning in stars, 
where only 0.7 % of the hydrogen mass is turned into energy in the production 
of helium. Thermal gas from the jet surroundings may also be entrained into 
the jet. The jets can propagate very large distances, e.g. more than 100 kpc in 
the case of Cygnus A, thus forming the largest physically connected structures 
in the universe (after cluster galaxies). At distances of less than a parsec, 
high-density (10^ cm^) gas clouds orbit the AGN and form the broad line 
region, or BLR. These clouds have speeds of several thousand kilometres per 
second and hence cause the linewidths observed by Seyfert in 1943. Further 
out, slower gas clouds constitute the narrow hne region, or NLR, with lower 
densities (10^ cm~^) and speeds of only a few hundred kilometres per second. 
Surrounding the AGN and its constituents in the polar plane is a toroidal 
agglomeration of material that can deeply hide the AGN and its activity. 

After two decades of hassle with AGN phenotypes it became clear in the 
1990s that there are basically two separate kinds of AGN classification: each 
object either is radio loud or radio quiet, and either belongs to the type 1 or 
2 AGN. The former classification is based on the ratio of radio luminosity to 
optical luminosity, R, with the dividing line being at around R = l (Kellermann 
et al. 1989). The latter classification is based on whether the optical emission 
fines are broad (type 1) or narrow(type 2). 

The cause of the bimodality in radio luminosity is still unknown, but sev- 
eral suggestions have been made, involving either "intrinsic" differences in the 
central engine or "extrinsic" differences in the surrounding medium. Intrinsic 
differences that have been suggested include (1) systematically lower black 
hole masses (Laor 2000), (2) lower black hole spins (Wilson & Colbert 1995), 
(3) a "magnetic switch" that was identified by Meier et al. (1997) during nu- 
merical modelling of jets, (4) the production of buoyant plasmons that bubble 
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up through the density gradient of the NLR instead of a coUimated relativis- 
tic jet (Pedlar et al. 1985, Whitde et al. 1986, Taylor et al. 1989), (5) a large 
thermal plasma fraction in the jet (Bicknell et al. 1998), or (6) radiative ineffi- 
ciency (Falcke & Biermann 1995). Extrinsic differences generally invoke the 
rapid deceleration of initially relativistic jets by collisions in a dense surround- 
ing BLR or interstellar medium (e.g. Norman & Miley 1984). Unfortunately, 
there are still only very few observational constraints, especially because the 
radio-weak objects are difficult to observe, and with so many possible causes, 
the question remains open why they have such low absolute luminosities. 

Unlike the difference between radio-loud and radio-quiet objects, the 
separation into type 1 and 2 AGN is understood as being due to an orientation 
effect. It depends on whether one can look into the central regions and see the 
innermost few tenths of a parsec, where the BLR clouds are, or whether the 
circumnuclear material shadows the BLR, in which case only narrow emission 
lines from the NLR are observed. Although each individual object has its 
peculiarities, virtually all AGN belong to one of the radio loud/quiet and type 
1/2 classes. 



Before going into the details of classification, a brief description is needed 
of three relativistic effects which are important in the understanding of AGN: 
relativistic beaming, Doppler boosting and apparent superluminal motion. The 
former two terms describe a directional anisotropy of synchrotron emission 
arising from charges moving at large fractions, j3, of the speed of light, c, 
towards the observer. If an object is moving at an angle 6 towards the observer 
with speed jSc so that the Lorentz factor, y, with 

(...) 



is >> 1, and if the object emits an isotropic flux density 5o(v), then the ra- 
diation is confined to a cone with half opening angle 1 /y (Fig. 1 .2) and the 
observer measures 

S{v)=S{v/d)d'=So{v)d'-<', (1.2) 
where a is the source spectral index defined here as 

a = i?#4^ (1.3) 

log(Vi/V2) 
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and 5 is the Doppler factor 

5 = y-\l-l5cosd)-\ (1.4) 

Thus, the radiation is not only confined to a smaller cone (beaming) but the 
observer also measures an increased flux density (Doppler boosting) because 
the source is moving nearly as fast as its own radiation. Both effects increase 
the flux density radiated in the forward direction and decrease it in the opposite 
direction by the same amount, so that even in relatively modest relativistic 
sources (yss 4), the jet to counter-jet ratio of flux densities can reach 10^. 

In the same geometric source configuration (small angle between the line 
of sight and the jet, relativistic speed of jet material), components can be seen 
moving at speeds exceeding the speed of light. This is a purely geometric 
effect: because the component is travelling almost as fast as the radiation it 
emits, the time interval between the emission of two photons appears shortened 
to us observers, and the apparent transverse speed can therefore exceed c. 




Figure 1 .2: Sketch of an electron moving in a magnetic field (reproduced from Peacoclc 1999). 

1.1.2 Seyfert Galaxies 

Seyfert galaxies belong to the class of radio-quiet AGN. Unlike their power- 
ful siblings, they rarely show clearly defined, linear radio jet structures on pc 
scales. They are mostly spirals and, although supermassive black holes have 
been found in some of them (e.g., 3.6 x 10^ M0 in NGC4258, Miyoshi et al. 
1995), their luminosities are only tiny fractions of that observed in quasars and 
radio galaxies. The sub-division into type 1 and type 2 Seyferts indicates the 
viewing angle: when the angle between the accretion disc normal and the line 
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of sight is intermediate or small (type 1), one has a direct view onto the BLR, 
and broad emission lines are visible. If the angle is large, circumnuclear mate- 
rial blocks the direct view onto the central engine and the BLR, and only NLR 
clouds are seen in the optical. Compelling evidence for this picture comes 
from galaxies like NGC 1068 which are classified as type 2 Seyferts, but show 
typical type 1 spectra in polarized light (Antonucci et al. 1994a; Antonucci & 
Miller 1985). In these cases, the type 1 spectrum is hidden by a foreground 
absorber, but reflections outside the absorber direct the light into our sight. 

1 .1 .3 Radio Galaxies 

Radio galaxies mostly exhibit single-sided, and only rarely two-sided jets on pc 
scales, but they almost always show two-sided radio emission on kpc scales, 
mostly in the form of bubble-like, irregularly shaped "radio lobes". Fifteen 
percent to 20 % of radio galaxies are radio-loud AGN (e.g., Kellermann et al. 
1989). On sub-pc to tens of pc scales, material is transported outwards in 
collimated jets which exhibit bright knots. On scales of kpc to hundreds of 
kpc, the jet flow becomes unstable, and either gently fades or abruptly stops 
in hot spots. Here, huge radio lobes evolve, and the jet material is slowly 
flowing back into the host galaxy. Because the material in the lobes is no 
longer relativistic, its emission is isotropic and the visibility of radio lobes is 
(to zeroth order) independent of the inclination angle. The single-sidedness 
on pc scales is due to relativistic effects, and only in a few sources where the 
jets have angles close to 90° with the line of sight, a double-sided structure on 
pc scales is seen. Another interpretation of the double-sidedness is that jets 
are always intrinsically single-sided, and that double-sidedness occurs when 
the direction of the jet rapidly flips from one side of the accretion disc to the 
other (e.g., Rudnick & Edgar 1984, Feretti et al. 1993). Similar to Seyfert 
galaxies, radio galaxies are divided into type 1 and type 2 based on optical 
appearance and as a result of the same geometric configuration. However, the 
terms Narrow Line Radio Galaxies (NLRG) and Broad Line Radio Galaxies 
(BLRG) are also established. A further sub-classification was established for 
radio galaxies based on their kpc-scale appearance by Fanaroff & Riley (1974). 
Analysing a sample of 57 radio galaxies and quasars from the 3CR catalogue, 
they discovered that the relative positions of regions of high and low surface 
brightness in the lobes of extragalactic radio sources are correlated with their 
radio luminosity. Fanaroff and Riley divided the sample into two classes using 
the ratio /?fr of the distance between the regions of highest surface brightness 
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on opposite sides of the central galaxy or quasar, to the total extent of the 
source up to the lowest brightness contour in the map. Sources with /?fr < 0.5 
were placed in class 1 and sources with 7?fr > 0.5 in class 11. It was found that 
nearly all sources with luminosity 

L(178MHz) < 2 X lO^^/z^o^jWRz"^ sr'^ (1.5) 

were of class 1 while the brighter sources were nearly all of class 11. The bound- 
ary between them is not very sharp, and there is some overlap in the luminosi- 
ties of sources classified as FR-1 or FR-II on the basis of their structures. The 
physical cause of the FR-1/11 dichotomy probably lies in the type of flow in the 
jets. FR-1 jets are thought to be subsonic, possibly due to mass entrainment, 
which makes them amenable to distortions in the interaction with the ambient 
medium, while the jets in FR-II sources are expected to be highly supersonic, 
allowing them to travel large distances. 



1 .1 .4 Quasars, BL Lacs, OVVs 



The objects in this section all belong to the radio-loud class. If the angle be- 
tween the jet axis and the Une of sight is small, one can see the BLR and the 
accretion disc directly. Relativistic effects are now dominating the phenotyp- 
ical properties of the AGN. The jet emission is focused into a narrow cone, 
and rapid (intra-day) variability might occur. Quasars reveal single-sided pc- 
scale jets and apparent superluminal motion of knots that travel down the jets. 
BL Lac objects are highly beamed, they show strong variability from radio to 
optical wavelengths and they have almost no optical emission lines, neither 
broad nor narrow. This happens when, at very small inclination angles, the 
thermal emission from the AGN is superimposed with the optical synchrotron 
emission from the jet base. At larger inchnation angles, the optical synchrotron 
emission is mostly beamed away from the observer, revealing the underlying 
thermal emission, and emission lines start to show up. This effect does not 
occur in the radio regime due to the absence of radio emission lines. The lack 
of optical emission lines makes a distance determination of BL Lacs difficult. 
The Optically Violently Variables, or OVVs, are a subclass of the BL Lacs, 
showing broad emission lines and rapid optical variability. 
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1 .1 .5 Compact Symmetric Objects (CSO) 

CSOs are double-sided, sub-kpc scale, symmetric radio sources whose spectra 
frequently have a peak in the GHz regime. Based on size and proper motion 
measurements, they are commonly regarded as being young (< 10"^ yr) radio 
galaxies whose jets have not yet drilled their way through the host galaxy's 
interstellar medium. 

1 .2 The Gas Around AGN 

The circumnuclear material in AGN is thought to settle in a probably rotation- 
ally symmetric body around the black hole. This gas is commonly referred to 
as the circumnuclear torus, although the toroidal shape has been established in 
only a few cases. Width, height and radius of this "torus" are not necessarily 
well constrained and are probably very different from object to object. In any 
case, however, the torus is expected to shield the central engine from the ob- 
server's view if the inclination angle is right. Observational evidence for tori 
comes from radio observations of total intensity and H I and molecular absorp- 
tion and molecular emission, from H2O maser observations, from the shape of 
ionized [OIII] and Ha-i-[Nn] regions, from X-ray observations and from the 
detection of thermal emission. A description of the physical properties of the 
circumnuclear gas can be found in, e.g., Krolik & Lepp (1989). 

1.2.1 Radio Absorption l\1easurements 

In those few cases where double-sided pc-scale radio jets are observed, a 
rather narrow "gap" in the emission across the source is frequently detected 
towards lower frequencies (e.g., in NGC 1052, Vermeulen et al. 2003, and in 
NGC4261, Jones et al. 2001). These gaps are mostly due to free-free ab- 
sorption, as identified by its characteristic frequency-dependence of the optical 
depth: 

Tv°=V-^-^ (1.6) 

In these cases, the UV radiation from the AGN ionizes the inner parts of 
a circumnuclear absorber, and the ionized medium then gives rise to free-free 
absorption. This situation is seen in, e.g., NGC 1052 (Kameno et al. 2001; 
Vermeulen et al. 2003), NGC 4261 (Jones et al. 2000, 2001) and Centaurus A 
(Tingay et al. 2001). A reliable detection of free-free absorption requires 
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flux density measurements at three frequencies at least. If, however, such 
measurements exist at only two frequencies, the spectral index a can be 
used to exclude synchrotron self-absorption (SSA). In SSA, the synchrotron 
radiation produced by the relativistic electrons is absorbed in the source, 
and if the distribution of the electron energies follows a power law, the 
spectral index cannot exceed +2.5. Thus, whenever a > 2.5 is observed at cm 
wavelengths, the absorption process is most likely free-free absorption (in the 
cases presented here, the Razin-Tsytovich effect has no effect, see Chapter 3). 

H I absorption has been found on sub-pc scales in, e.g., NGC 1052 (Ver- 
meulen et al. (2003)), Cygnus A (Conway & Blanco 1995) and in NGC 4261 
(van Langevelde et al. 2000). In the latter, the H I absorption was found no 
closer than 2.5 pc away from the core, supporting the idea that the material 
is ionized at smaller distances to the AGN. This was proposed by GalUmore 
et al. (1999), who detected H i absorption in a number of Seyfert galaxies 
almost exclusively towards off-core radio components. 

A variety of molecular lines is also detected in AGN: NGC 1052 shows 
OH in absorption and emission (Vermeulen et al. (2003)), in NGC 4261, Jaffe 

& McNamara (1994) detected CO in absorption and Fuente et al. (2000) found 
CO+ absorption in the core of Cygnus A. This list is by no means complete 
and many other Une measurements exist. 



1 .2.2 Masers 

The most compelling evidence for the partly molecular nature of at least parts 
of the circumnuclear material comes from the detection of H2O "megamasers" 
in the vicinity of AGN. A few prominent examples are NGC 3079 (Henkel 
et al. 1984, and, more recently, Kondratko 2003), NGC 2639 (Wilson et al. 
1995) and NGC 4258 (Miyoshi et al. 1995). In the last case, the velocity dis- 
tribution of the masers yielded a model-independent measure for the enclosed 
mass with high accuracy. In NGC 1068, both the alignment and the velocity 
gradient of the masers found by Gallimore et al. (1996) are oriented perpen- 
dicular to the radio jet axis, hence suggesting the presence of a circumnuclear, 
relatively dense region of material. The physical properties of H2O in AGN 
have been described by Neufeld et al. (1994). 



1.3. Feeding Gas Into the AGN 



21 



1 .2.3 Ionization Cones 

In the last 15 years, optical observations of AGN showed cone-shaped regions 
of Une emission (Pogge 1988), frequently in two opposite directions (Storchi- 
Bergmann et al. 1992; Wilson et al. 1993) and sometimes aligned with linear 
radio structures (Falcke et al. 1998). The shape of the line emitting regions was 
immediately interpreted as being due to shadowing by circumnuclear material, 
i.e., the nuclear UV emission can only escape towards the poles of the torus, 
where it ionizes the gas. 

1 .2.4 X-ray Observations 

Evidence for circumnuclear absorbers also comes from X-ray observations. 
Marshall et al. (1993) showed that the X-ray spectrum of NGC 1068 is 
best modelled as a continuum source seen through Compton scattering. In 
Cygnus A, Ueno et al. (1994) found evidence for an absorbed power-law spec- 
trum. Seyfert 2 galaxies are mostly heavily absorbed in the X-ray regime, with 
column densities of the order of > lO^'^cm"^ (Krohk & Begehnan 1988). 

1.2.5 Thermal Emission 

In NGC 1068, Gallimore et al. (1997) have discovered a region of flat-spectrum 
radio continuum emission with brightness temperatures too low to be due to 
self-absorbed synchrotron emission. The shape and orientation of the region is 
suggestive of a circumnuclear disc or torus, and they conclude that the emis- 
sion is either due to unseen self-absorbed synchrotron emission that is reflected 
by a torus into the Une of sight or free-free emission from the torus itself. 

1 .3 Feeding Gas Into the AGN 

1.3.1 Magnetic Fields 

A lot of material is seen in the vicinity of AGN, but how is the material moved 
closer to and into the black hole and into the jets? Shlosman et al. (1989) 
suggested that a stellar bar in a galaxy sweeps gas into the central few hundred 
pc. The gas forms a disc which also develops a bar potential and funnels the 
gas further in to scales of a hundred parsecs. Processes that transport the gas 
further in to parsec scales are largely unknown, and no theory exists that could 
be tested by observations. To move gas into smaller radius orbits requires a 



22 



1.3. Feeding Gas Into the AGN 



mechanism to shed angular momentum, and the candidates are viscosity and 
magnetic fields. Viscosity, responsible for angular momentum transport in 
the innermost regions of AGN, would form the gas into a disc, which is not 
generally observed. Only in few objects have such discs with diameters of 
~ lOOpc been found (e.g., in NGC4261, Ferrarese et al. 1996; Jaffe et al. 
1996 and in Mrk231, Klockner et al. 2003), but the role that the discs in these 
objects play in gas transport is unclear. 

Magnetic fields, on the other hand, play an essential role in models of tori, 
accretion discs and jet formation in the central parsec and sub-parsec scale 
regions (e.g., Koide et al. 2000; Krolik & Begelman 1988; Meier et al. 2001). 
The basic idea is that the accretion disc is threaded with magnetic fields, and 
that differential rotation of the disc twists the fields to a spiral structure, which 
applies a breaking torque on the inner disc and accelerates material in the outer 
disc, thus transporting angular momentum outwards. 

The origin of the required magnetic fields is not clear. They are proba- 
bly frozen into the accreted material on scales of > lOpc, although parts of 
the field may be lost in reconnection. Dynamo processes in the disc can also 
produce magnetic fields (e.g., von Rekowski et al. 2003) which participate in 
the generation of an outflow, but the details are not well known. Estimates 
of magnetic field strengths in AGN and their surroundings mainly come from 
equipartition arguments, assuming that the energy in particles equals that in 
magnetic fields. Whether this assumption is justified or not is not known, and 
so equipartition arguments are unsatisfactory. 

Magnetic fields are routinely observed on the largest scales in galaxies and 
are modelled on the smallest, but are largely unknown on intermediate scales 
of a few parsecs to several tens of parsecs. As they are probably involved in the 
transport of gas into the AGN, their strength and orientation is of considerable 
interest. 

1 .3.2 Faraday Rotation and Free-Free Absorption 

One way to measure magnetic field strengths is by means of Faraday rotation. 
When an electromagnetic wave travels through an ionized medium that is 
interspersed with a magnetic field with a component parallel to the direction 
in which the wave is travelling, then the plane of polarization of the wave is 
rotated. The amount of Faraday rotation cannot be measured directly because 
the intrinsic position angle of the polarization is not known and the rotation 
has ambiguities of n. Instead, one exploits the frequency dependence of 
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the Faraday rotation, and the change of the effect with frequency yields the 
constant of proportionality, the rotation measure RM. 

The physics behind this effect is the birefringence of the magnetized 
plasma. A linearly polarized beam of radiation with electric vector position 
angle 6 can be considered as the superposition of two circularly polarized 
waves with equal amplitudes but opposite senses of rotation. The circular po- 
larizations have different indices of refraction in the plasma which causes one 
polarization to be retarded with respect to the other, and the plane of linear 
polarization, composed of the two circular polarizations, rotates. 

Faraday rotation is of particular importance in this thesis so we give a brief 
derivation of the formula here, following Kraus & Carver 1973, p. 737f. 

In a magnetized plasma, where the magnetic field is parallel to the direction 
of propagation of the waves, the phase constants jS"*" and j3~ of the two waves 
are given by 

I5^ = (j0y/^{en±£n), (1-7) 

where co is the angular frequency of the wave, /io is the vacuum permeability, 
and en and en are elements of the permittivity tensor e, with 



£11 



£l2^ 



1 + - 



(On 



(0^ 



£o 



(1.8) 



Here, coo = ey^n^Jsom denotes the angular plasma frequency (e is the par- 
ticle charge in C, «e the particle density in m^^ and m the particle mass in 
kg), cOg = {e/m)B the angular gyrofrequency and Eq the vacuum permittivity. 
Travelling a distance dl through the plasma changes Q by 

dd = f—^dl. (1.9) 

Inserting the expressions for (Oq and cOg into Eq. 1.7 and integrating over 
the line of sight then yields 
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Sn^c^eom^ Jo " 
= X^x 2.63 X 10-^^ J^n^B\idl ^^'^^^ 
= X RM, 

where A is the wavelength in m, e is the elementary charge in C, c is the speed 
of light in ms~^ £0 is the vacuum permittivity in Fm^^ m is the electron mass 
in kg, rie is the electron density in m~^, S|| is the line-of-sight component of 
the magnetic flux density in T, and / is the path length in m. 

Observations of Faraday rotation yield a measure of the integral of the 
product ?iefi||- To separate the magnetic field from the electron density 
requires an independent measurement of the electron density and the path 
length through the ionized gas. 



Ionized gases also produce free-free absorption which causes an exponen- 
tial decrease of intensity with path length through the absorber. The spectral 
energy distribution of synchrotron radiation which is absorbed by free-free ab- 
sorption, is given by 



5v=5oV«»xexp(-<), (1.11) 

where 



= 8.24 X IQ-^r-i-^^ v"^-^ J N+N-ds (1.12) 

(e.g., Osterbrock 1989, eq. 4.32). Here, 5v is the observed flux density in mJy, 
^0 is the intrinsic flux density (before the radiation passes the absorber) in mJy, 
V is the observing frequency in GHz, Oq is the dimensionless intrinsic spectral 
index, T is the gas temperature in K, and are the number densities of 
positive and negative charges, respectively, in cm~^, and s is the path length in 
pc. 

Faraday rotation depends linearly on the electron density, n^, whereas the 
optical depth in free-free absorption goes as nl (assuming = N-= n^). An 
analysis of the frequency-dependence of Ty , together with a diameter measure- 
ment of the absorber from VLBl images, yields n^. Putting and the diameter 
into Eq. 1.10 then allows one to solve for the magnetic field strength, Sy. 
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1.3.3 Jet Collimation 

As was mentioned in the section §1.3.1, models exist for the innermost re- 
gions of AGN where the jets are launched and coUimated. Some possible 
mechanisms for the acceleration of jets are the "magnetic slingshot" model 
(Blandford & Payne 1982), through the extraction of Poynting flux from the 
black hole, the so-called Blandford-Znajek mechanism (Blandford & Znajek 
1977), or by radiation or thermal pressure (Livio 1999; Odell 1981). 

Unfortunately, although today's radio interferometers routinely yield an- 
gular resolutions of 0. 1 mas, these regions are still not resolved, except for a 
few nearby objects with high black hole masses, like M87 (Junor et al. 1999). 
Jet collimation is expected to happen on scales of 10 to 1000 Schwarzschild 
radii (/?s), and in a moderately distant radio galaxy with redshift z = 0. 1 and 
a typical black hole with 10^ Mq, these scales are still factors of 20 to 2000 
smaller than the synthesized beam. Hence, only little observational evidence 
exists so far to constrain jet formation models. In a few nearby AGN, however, 
the highest resolution VLBI observations at 86 GHz can in principle resolve 
scales of < 1000 /?s- These observations are challenging because most nearby 
sources are weak at high frequencies and the antenna sensitivities decrease, 
and most attempts to detect the targets have failed. 

1 .4 The Aim of This Thesis 

Nearby AGN provide unique opportunities to study the circumnuclear envi- 
ronments of supermassive black holes. VLBI observations have long been 
concentrated on the brightest, but unfortunately more distant objects, yielding 
rather low linear resolutions of more than one parsec. In this thesis, I present 
VLBI observations primarily of nearby objects, yielding some of the highest 
linear resolution images ever made. The goal was to investigate the distribu- 
tion of gas and magnetic fields in those objects, and to probe the jet collimation 
region of a nearby radio galaxy with highest Unear resolution. This required 
the development of fast frequency switching for phase calibration of 86 GHz 
VLBI observations of weak sources. 

1 .4.1 Polarimetric Observations of Six Nearby AGN 

Six nearby AGN were selected because they show good evidence for circum- 
nuclear free-free absorbers that shadow the radiation of the AGN. If such ab- 
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sorbers are interspersed with magnetic fields and the radio emission from the 
AGN is polarized, Faraday rotation is expected to occur. A joint analysis of 
the free-free absorption and Faraday rotation then allows one to determine the 
magnetic field strength in the absorber and so yields a measurement of a quan- 
tity which is difficult to determine by other means. I present the results of 
pilot VLBI observations that were made to look for polarized emission before 
making time-consuming Faraday rotation measurements. 

1 .4.2 Case Study of NGC 3079 

In addition to the polarimetric observations of NGC 3079, 1 present an analysis 
of multi-epoch, multi-frequency observations of this nearby Seyfert 2 galaxy to 
investigate the nature and origin of the radio emission. Furthermore, I present 
a statistical analysis of VLBI observations of Seyfert galaxies reported in the 
hterature to compare NGC 3079 to other Seyferts and to compare Seyferts to 
radio-loud objects to investigate the difference between powerful and weak 
AGN. 

1 .4.3 Phase Calibration Strategies at 86 GHz 

We have explored the feasibiUty of a new phase calibration strategy for VLBI 
observations, in which one cycles between a lower reference frequency, at 
which the source is strong enough for self calibration, and the target frequency. 
The phase solutions from the reference frequency are scaled by the frequency 
ratio and interpolated onto the target frequency scans to remove the atmo- 
spheric phase fluctuations. The result is a phase-referenced image at the target 
frequency, and indefinitely long coherent integrations can be made on sources 
that are too weak for self-calibration. The primary use of the technique is to 
image nearby, weak AGN at 86 GHz to obtain highest linear resolutions. Using 
data from a pilot project, we have improved the observing and data calibration 
strategy to a ready-for-use level, and we have obtained the first detection of 
NGC 4261 which could not be imaged previously at 86 GHz because it is too 
weak. 



Chapter 2 



Special VLBI Techniques 



Very Long Baseline Interferometry, or VLBI, has become a standard observ- 
ing technique in the radio astronomy community during the last 20 years. Al- 
though frequently referred to as "experiments", VLBI observations have left 
the experimental stage, and especially observations with the U.S. 10-element 
Very Long Baseline Array, or VLB A, are relatively easy to prepare and anal- 
yse, and now the first automated data calibration methods exist (Sjouwerman 
et al. 2003). A good description of VLBI principles and the data calibration 
steps required has recently been given by Klare (2003), and details can be 
found in, e.g., Thompson et al. (1986). I therefore will not go into the details. 

There are, however, non-standard VLBI observations that require special 
observing procedures and/or data calibration steps. Two of these techniques, 
phase-referencing and polarimetry, have been used to gather data for this thesis 
and therefore are described in detail. Furthermore, a short section is dedicated 
to a description of ionospheric phase noise in phase referencing and how it can 
be caUbrated. 

2.1 Phase Referencing 

In VLBI, the true visibility phase is altered by errors from numerous sources. 
They can be divided into geometric, instrumental and atmospheric errors. Geo- 
metric errors arise from inaccurate antenna positions, motion of antennas with 
tectonic plate motion, tidal effects, ocean loading and space curvature due to 
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the mass of the sun and the planets close to the targeted position. Instrumental 
errors include clock drifts, changes in antenna geometry due to gravity forces, 
cable length changes, and electronic phase errors due to temperature variations. 
Most of these errors are sufficiently well known or are measured continuously 
and are accounted for in the correlator model of the array. Atmospheric and 
ionospheric phase noise, however, is difficult to predict and is therefore the 
largest source of error in VLBl observations. The following description of 
the effect of atmospheric and ionospheric phase noise on VLBI observations 
closely follows Beasley & Conway (1995). 



2.1 .1 Tropospheric Phase Noise 

At cm wavelengths, the largest source of error in VLBI observations comes 
from fluctuations of the tropospheric water vapour content along the line of 
sight of the telescopes. Changes in the water vapour content cause phase 
changes of the observed visibilities and thus limit the atmospheric coherence 
time. This is the time over which data can be coherently averaged, and is taken 
to be the average time it takes for the phase to undergo a change by one radian 
57°). 

The tropospheric excess delay, i.e., the additional time it takes for the 
waves to travel through the atmosphere compared to vacuum, can be divided 
into two components, the dry troposphere and the wet troposphere. This means 
that the mixture of air and water vapour would have the same effect on the 
visibility phases as a layer of dry air and a layer of water vapour, both of thick- 
nesses equivalent to the tropospheric content of air and water vapour, respec- 
tively. The zenith excess delay due to the dry troposphere, Iz, can be modelled 
quite accurately using measurements of pressure and antenna latitude and alti- 
tude (Davis et al. 1985): 

^ 0.228 m/mbarPp 

^ 1 -0.00266 cos (2A)- 0.00028 km- i/i' ^ ' ^ 

Here, is the zenith excess path length in cm, is the the total pres- 
sure at the surface in mb and h is the altitude of the antenna above the geoid 
in km. With this equation, the dry troposphere, contributing around 2.3 m, 
can be modelled to an accuracy of 0.5 mm. The equation can be extended 
to include the tropospheric water vapour, yielding the Saastamoinen model 
(Ronnang 1989) 
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_ 0.228m/mbar [Po + Pw(^ + 0.05)] 

~ 1 - 0.00266 cos (2A ) - 0.00028 km- 1 h ' ^^'^^ 

Here, P„ is partial pressure of water vapour and T is the temperature in 
K. However, the effects of the wet troposphere are much more difficult to de- 
scribe because the water vapour is not well mixed with the dry air, and turbu- 
lence makes the delay due to the wet troposphere highly variable. In general, 
tropospheric water vapour contributes up to 0.3 m of excess path delay, but to 
determine the exact value requires precise measurements of the water vapour 
in front of each telescope. Water vapour radiometers therefore have become in- 
creasingly popular, especially at radio telescopes operating at mm wavelength. 
The new Effelsberg Water Vapour Radiometer (Roy et al. 2003) is now able to 
measure the delay to an accuracy of 0.12 mm (1 a), corresponding to, e.g., 2.2° 
at 15.4 GHz. Further improvements will shortly increase the accuracy to about 
0.04 mm and hence allow one to use the water vapour radiometer to caUbrate 
VLBl data at frequencies of up to 86 GHz. 

The path length through the atmosphere scales as sin^'(j8), where j3 is 
the antenna elevation. Thus, the amount of troposphere along the Une of 
sight has doubled at an elevation of 30°, and has increased to five times its 
zenith value at 12° elevation. As a consequence, the phase noise dramatically 
increases towards low antenna elevations. At frequencies above ~ 5 GHz, 
tropospheric phase noise is the dominant source of error in VLBI observations. 



2.1.2 Ionospheric Phase Noise 

The ionosphere is a region of free electrons and protons at altitudes of 60 km 
to 10000 km above the earth's surface. It adds phase and group delays to the 
waves and causes Faraday rotation of linearly polarized waves. The excess 
zenith path in m is (Thompson et al. 1986, eq. 13.128) 

40.3 Hz^m^^ 
k ^ -2 TEC, (2.3) 

where v is the frequency in Hz and TEC is the vertical total electron con- 
tent in m~^. The TEC changes on various timescales. Long-term variations 
are caused predominantly by changing solar radiation (solar cycle, seasonal 
and diurnal variations), but short-term variations are mostly caused by atmo- 
spheric gravity waves through the upper atmosphere, oscillations of air caused 
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by buoyancy and gravity. The biggest effect on radio observations are medium- 
scale travelling ionospheric disturbances (MSTIDs) with horizontal speeds of 
lOOms^' to 300 ms^\ periods of lOmin to 60min and wavelengths of several 
hundred km. The ionosphere introduces a delay that scales as v^^, whilst the 
phase scales as v~\ and this contribution is the dominant source of phase er- 
ror in VLBl observations at frequencies below ~ 5 GHz. Ionospheric Faraday 
rotation is negligible at cm wavelengths, with a maximum of 15 turns of the 
electric vector position angle at 100 MHz (Evans & Hagfors 1968) and hence 
at most 9 turns, 3 turns and 1 turn of phase at 1 .7 GHz, 5.0 GHz and 15.4 GHz, 
respectively. 

2.1 .3 Phase Referencing 

In standard VLBI observations at cm wavelengths, the so-called phase 
self-calibration is used to solve for phase errors that cannot be accounted for 
by the correlator model and mostly are tropospheric. Starting with a point 
source model in the field centre, and refining the model iteratively, correction 
phases are derived that make the visibility phases compliant with the model. 
This procedure is known as self-calibration, or hybrid mapping (Comwell 
& Wilkinson 1981). But in two cases, this procedure does either not work 
or is not desirable. In the first case, the target source is weak and cannot be 
detected reliably within the atmospheric coherence time. As an illustration, 
in a typical 15 GHz VLB A experiment with ten stations in moderate weather, 
the coherence time is 60s and the 7a detection limit follows to 21mJy, 
and weaker sources cannot be detected (a detailed calculation of sensitivity 
limits is given in Chapter 6). In the second case, one is interested in absolute 
astrometry of the target source. A priori phase self-calibration destroys that 
information because the visibility phases are initially adjusted to fit a point 
source at a position provided by the observer. Literally, one can "move 
around" the source in the field of view and one is always able to properly 
adjust the visibility phases. Phase-referencing solves both of these problems. 

In brief, phase referencing uses interleaved, short observations ("scans") 
of a nearby calibrator to measure the tropospheric phase noise which is sub- 
tracted from the target source visibilities. It works only if the interval between 
two calibrator scans is shorter than the atmospheric coherence time and the 
calibrator structure is well known. The first demonstration of this technique 
was published by Alef (1988) (a similar approach had already been proven to 
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work by Marcaide & Shapiro 1984, but in their case, the calibrator and tar- 
get were in the telescope primary beams, and no source switching was needed). 

Consider an observing run in which each target source scan of several min- 
utes is sandwiched between short scans on a nearby caUbrator. Let us further 
assume that all a priori calibration information has been applied, i.e., the ampU- 
tudes are calibrated and the bandpass shape and instrumental delay and phase 
offsets are corrected for. The measured visibility phases can then be described 
with 

</>cal (fl ) = <^c (f 1 ) + Cs + <^>pos 1 ) + 0ant 1 ) + 0atm 1 ) + Cn (h ) 

(f>tM = 0t(?2) + 0ins(^2) + 0pos(?2) + 0L(?2) + KM + Konitl) (2.4) 

$cdt3) = Uh) + 0ins(?3) + 0pos(^3) + Ct(?3) + ^.tjt^) + Cnl^s)- 

(^c and are the true visibility phases on the calibrator and the target 
source, (^ins is the residual instrumental phase error due to clock drifts and other 
electronics, (^os and are geometric errors arising from source and antenna 
position errors, and 0atm and ^lon are tropospheric and ionospheric phase noise 
contributions. Using self-calibration, the measured visibility phase is decom- 
posed into the source structure phase and the difference in antenna-based phase 
errors at ti and tj, and interpolation yields the calibrator visibihties at t2: 

M2) = ^c{h) + &(^2) + ^pos(^2) + &t(^2) + ^atm(^2) + #on(^2) (2.5) 

(a tilde denotes interpolated values). Subtracting the interpolated calibrator 
visibility phases from the observed target visibility phases at time t2 gives 

- 0cal = (^t - ^c) + (Cs - Cs) + (0pos - 0pos) 

+ (•/'Int - Ct) + (tim - 0atm) + (Cn " Cn)- (2-6) 

In this equation, most terms are zero. Residual instrumental errors vary 
slower than the lag between the calibrator and target source scan, and therefore 

'^'ins ~ ^kis ~ ^- Atmospheric and ionospheric noise will also be the same in 
adjacent scans if the separation between calibrator and target is less than a 
few degrees and the lines of sight pass through the same isoplanatic patch, the 
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region across which tropospheric and ionospheric contributions are constant. 
This yields 0^^^ — ^^f^ = and (pl^^ — 0-^^^ = 0. Antenna position errors have 
the same effect on both the calibrator and the target if their separation on the 
sky is small, and hence (^^, — = 0. The calibrator visibility phases change 
slowly, and hence 0cai = 0c- Also, when using either a compact calibrator or a 
good calibrator model in phase self-calibration, the phase errors derived from 
phase self-caUbration are those of a point source, and hence = 0. This yields 

<?>tar - ^cal = + (^pos " 0pos) + 0int, (2.7) 

in which (pjm denotes interpolation errors. This equation expresses that after 
interpolation, the difference between the calibrator and the target source 
phases is the target source structural phases plus the position error Hence, 
phase-referencing not only allows one to calibrate the visibility phases, 
but also to precisely measure the target source position with respect to the 
calibrator position. 

In general, the ionospheric phase component can be calibrated in self- 
calibration, together with the tropospheric component. But especially at low 
elevations and frequencies below 5 GHz, the assumption that ^j'^^^^ — ^^^^ = is 
no longer valid, and a correction has to be applied. One approach to correct 
for the ionospheric delays is to use TEC models derived from GPS data. They 
are provided by several working groups, yielding global TEC maps every 
two hours, and giving the TEC on a grid with 2.5° spacings in latitude and 
5° spacings in longitude. From this grid, the TEC at each antenna can be 
interpolated. The error in these maps, however, can be quite high, up to 20 % 
when the TEC is as high as a few tens of TEC units (1 TECU=10^^e^m~^), 
and up to 50 % or higher when the TEC is of the order of a few TEC units. 
Details on how the TEC is derived from GPS data are discussed by Ros et al. 
(2000), and a set of tests is described in Walker & Chatterjee (2000). 



2.2 Polarimetry 

Measuring the linear or circular polarization of radio waves is a relatively 
young technique in VLBI. The first polarization-sensitive VLBI observation 
of an AGN jet was published by Cotton et al. (1984) (unfortunately, their 
3C 454.3 image was rotated by 180° because they had mistaken the phase 



2.2. Polarimetry 



33 



signs). But VLBI polarimetry has become increasingly popular, leading to ex- 
tensive surveys in the last few years (e.g., Pollack et al. 2003; Zavala & Taylor 
2003). This section describes technical details and calibration of polarization- 
sensitive VLBI observations. 



2.2.1 Stokes Parameters 




Figure 2.1 : The Poincar^ sphere, reproduced from Scott et al. (2001). 



The polarization of any source can be described by means of the Stokes 
parameters /, Q, U and V. I denotes the total intensity of the source, Q and 
U the fractions of linear polarization perpendicular to the propagation of the 
wave and in directions that enclose angles of 45°, and V denotes the fraction of 
circular polarization. A convenient way to look at this is the Poincare sphere 
(Fig. 2. 1). Consider a Cartesian coordinate system in which Q is on the x axis, 
U is on the y axis and V is on the z axis. The radius of the sphere is the degree 
of polarization, 

/p = VQ^ + U^ + V\ (2.8) 

and the points on the sphere represent the different states of polarization. 
The intersection points of the Q, U and V axes with the sphere represent the 
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following states of polarization, respectively: linear, horizontal (+Q), linear, 
vertical (-0, linear, +45° {+U), linear -45° (-U), circular, left-handed 
i+V) and circular, right-handed {—V). In general, the polarization is a 
combination of all three parameters, and the polarization is elliptical. 



2.2.2 Interferometer Response to a Polarized Signal 

Determination of an extended source's polarization characteristics requires 
measuring the Stokes parameters over the source region. An interferometer 
with coordinates (m,v) in the plane perpendicular to the line of sight measures 
the Fourier transform of the sky brightness distribution, given in angular coor- 
dinates ( C 5 ^ ) • Similar to the Fourier transform of the total intensity brightness 
distribution, /, the interferometer response to the sky distribution of the other 
Stokes parameters, Q, U and V can be defined: 



Starting with the interferometer response to a point source, 1 now describe 
how the quantities /, Q, U and V are restored from the visibility measure- 
ments. I then describe the calibration of polarization-sensitive interferometer 
data. The polarization calibration procedure used in this thesis was developed 
by Leppanen et al. (1995). I follow their notation to explain the interferometer 
response and the correlation to the Stokes parameters; the full derivation can 
be found in Leppanen (1995). 




J^(m,v) 



y{u,v) 




(2.9) 



Response to a Point Source 

The voltages A present at an antenna output can be described by 
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(2.10) 



Here, the g's are complex quantities proportional to the antenna gains, £"s 
are the RCP and LCP parts of the electric field, respectively, D's are the com- 
plex leakage of power from one polarization into the other, and e^^" describe 
the rotation of the wave with respect to the feed horns due to the parallactic 
angle of the source, CC. These terms are commonly referred to as D-terms. Be- 
cause virtually all sources exhibit fractional polarizations that are of the same 
order as the D-terms on the scales probed with VLBl, the significance of the 
detection of polarized emission depends critically on the D-terms being prop- 
erly caUbrated. All quantities, except for D, are functions of time. 

The cross-correlation function, p, between two antennas is given by 



This equation expresses an average of the multiplication of the antenna 
voltages. A, from two antennas, m and n, with a certain delay, T. The indices p 
and q represent either LCP or RCP, yielding four possible combinations for 
each baseline. The functions A are real, and hence their cross-correlation 
is equal to their convolution. One gets the cross-power spectrum, r(co), by 
Fourier-transforming Eq. 2.11. Applying the convolution theorem then yields 



where the ^ (ffl) are the Fourier transforms of the terms in Eq. 2.10. Eq. 2.12 

is valid for any antenna pair because we still assume a point source in the field 
centre, and the interferometer amplitude therefore is independent of basehne 
length and orientation. 

When Eq. 2.10 is inserted into Eq. 2.11 before Fourier transformation, 
Eq. 2.12 contains the terms E^E^* E^E^* E^E^*, which we denote 

by RR, LL, RL and LR, and their Fourier transforms by and 
££M. The four cross-power spectra then take on the form 



(2.11) 



(2.12) 
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+ D^* X J^^ + D^D^*e-^'(""'-"") X ^M] 

+£)«*g;(«»-«„) X + X ^^]. 

This particular representation of cross-power spectra, electric vectors and 
leakage factors is called the "leakage-term model" (Cotton 1993). The equa- 
tions become a lot simpler considering that sources observed with VLBI are 
only weakly polarized (thus and terms are small) and that the D- 
terms are small (thus DD terms are even smaller). Then, in Eq. 2.13, all terms 
involving products of or MJ/f and D-terms, as well as those involving 
DD terms, can be neglected, yielding 



flL — „L „L* j(a,„-a„) (/J (/> 
' mn Smon ^ ^ =^ ^ 

rin = sign* X I^^+D^ej^"--""^ x 

+ D^*e-j^""^-""^ x^£g] 

+ £)fV(«--'^)x^Jf]. 



Response to an Extended Source 

For an extended source, the cross-power spectrum is the superposition of in- 
coherent contributions from all parts of the source, and integration over the 
source region is required: 



rP%u,v,(o)=J J F'^ii;,'n,(0)e--'^''^^"+'^''^di;dri. (2.15) 
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This equation shows that the cross-power spectra of extended sources ob- 
served by an interferometer are the Fourier transforms of the brightness distri- 
bution in sky coordinates: 

rP''{u,v,co) = ^{I'"'{^,ri,(o)}. (2.16) 
It is straightforward to enhance Eq. 2. 14 for extended sources, yielding 

+Dfe^'(«'"-«")x^^]. 

We now need to establish a relation between the cross-power spectra in 
Eq. 2.17 and the Stokes parameters. The Stokes parameters and the cross- 
correlation of the electric field components are connected through 





= RR 


= /+y 




= LL 


= 1-V 


E^E^* 


= RL 


= Q + jU 


E^E^* 


= LR 


= Q-jU 



(2.18) 



(Thompson et al. 1986), from which the Stokes parameters can be separated 
into 



I = -{RR+LL) 
Q = ^{RL + LR) 
U = j^{LR-RL) 
V = ^{RR-LL). 



(2.19) 
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Applying the inverse Fourier transform to Eq. 2.17 then yields the desired 
sky distribution of the Stokes parameters. 

2.2.3 Calibration of Instrumental Effects 

Given the usually low degrees of polarization in extragalactic sources, mea- 
suring the four Stokes parameters requires an accurate calibration of the array. 
There are six sources of error in VLBl polarization observations: 1) D-term 
calibration errors; 2) errors in the relative phase between RCP and LCP (R-L 
phase offset); 3) thermal noise; 4) gain calibration errors; 5) deconvolution 
errors and 6) closure errors. 4), 5) and 6) are usually small, because the dy- 
namic ranges in polarization images are low, and therefore are neglected, and 
3) cannot be calibrated, but integration time is planned to make this small. 
This leaves 1) and 2) as the dominant sources of error. 2) is easily caUbrated 
(in the case of the VLB A) using the station monitoring data, and 1) requires 
significant effort. 

The R-L Phase Offset 

The relative phase between the two circular polarization receiving channels 
affects the apparent position angle of the linear polarization on the sky, or 
electric vector position angle (EVPA). Also, the relative phases of adjacent 
baseband (or IF) channels need to be Uned up to allow averaging over the ob- 
serving bandwidth. At the VLBA, pulses are injected at the start of the signal 
path with a period of 10^^ s. In the frequency domain, the pulses produce a 
comb of lines spaced by 1 MHz, bearing a fixed, known phase relationship to 
each other. They pass through the receiver and downconversion chain along 
with the radio astronomical signal and the phases of the pulses are measured 
in the back-end before the signal is digitized and written on tape. This allows 
one to derive the time-dependent, instrumental phase changes over the band 
(the instrumental delay) and the offsets between the IFs. Basically, the pulse 
caUbration system moves the delay reference point from the samplers to the 
pulse calibration injection point at the feeds, reducing instrumental phase 
changes almost to zero. This also holds for the phase offsets between the two 
parallel hands of circular polarization and between the cross hands. Leppanen 
(1995) has shown that, after application of the pulse caUbration, the residual 
R-L phase errors are of the order of 1°. 
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2.2.1 Parallel-hand spectrum of 3C345 at 5.0GHz without pulse calibration applied, showing a 
1 min-average of the visibiHty phase (deg, upper panels) and ampHtude (Jy, lower panels) in 14 
out of 16 observed frequency channels. The two highest frequency channels in each IF have been 
flagged due to the bandpass limitations. Instrumental delays dominate the phase slopes across the 
pass band and the phase offsets between the IFs. 
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2.2.2 The same spectrum as above with phase corrections derived from the pulse calibration sys- 
tem. The phase slopes have been removed and the phases across the band have been aligned, 
allowing coherent integration over the observing bandwidth. 



Figure 2.2: Illustration of the effect of the pulse calibration system 
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In practice, applying the pulse calibration is straightforward. A pulse cal- 
ibration, or PC table is generated at the correlator and is attached to the (m,v) 
data. The AIPS task PCCOR is then used to generate a calibration table with 
phase corrections. PCCOR takes two phase measurements per IF, at the upper 
and lower edge of the band, and computes the delay. To resolve phase ambi- 
guities, one needs to specify a short scan on a strong source from which the 
delay is measured using a Fourier transform with a finer channel spacing than 
that used by the pulse calibration system. Once the ambiguity is resolved and 
assuming that the instrumental phase changes are small throughout the obser- 
vation (which is virtually always true), PCCOR computes phase corrections 
from the PC table. The data are then prepared for averaging in frequency to in- 
crease the signal-to-noise ratio, and if polarization measurements are desired, 
the EVPA is calibrated. 



Calibration of D-terms 

In brief, the calibration of D-terms uses the effect that the linearly polarized 
emission from the source rotates in a different way with respect to the feed 
horns than does the leakage polarization as the earth rotates and the antennas 
track the source. Let us multiply each line in Eq. 2.17 by a power of e such 
that those terms unaffected by the D-terms are unrotated and the power of the 
e-factor is zero (e.g., fQj. i^j^g third Une). This yields 



r^nM =gi^*[^^ + Die-^i<^ X ^M + D^^e-^i"^ x ^^], 

in which only the leakage terms rotate twice as fast with the parallactic angle 
of the source, and the visibility is independent of parallactic angle. This 
situation can be understood as a vector of a certain length and position angle, 
representing the linearly polarized emission from the source, to which the 
vector of the polarization leakage, rotating with the parallactic angle of the 
source, is added. Separation of these vectors yields the D-terms. 



Two direct consequences can be drawn from this picture: the smaller 
the D-terms, the more difficult they are to determine because the relative 
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contribution of thermal noise to the polarization leakage vector increases, 
and the larger the range of parallactic angles over which a source has been 
observed is, the more accurate the Z)-term calibration is possible. 

Separation of the true linear polarization and the leakage polarization is 
implemented in the AIPS task LPCAL, which works as follows. One can think 

of the leakage polarization as the convolution of a "leakage beam" with the 
total intensity structure. In such an image, unphysical features would appear 
as linear polarization. Unfortunately, unlike in the total intensity images, one 
cannot use a positivity constraint for Stokes Q and U because they can have 
either sign. However, polarized emission should usually appear only at loca- 
tions in the image where total intensity emission is seen, and LPCAL uses this 
constraint (the "support") to derive the D-terms. The support is provided to 
the algorithm in terms of a total intensity model derived from self-calibration. 
In general, the polarized intensity structure differs considerably from the total 
intensity structure, and so the model must be divided into sub-models, small 
enough that the polarization of each can be described by a single complex 
number pc = {Qc + jUc)/Ic. 

The visibilities of each of the sub-models, /, can now be written as pi^i. To 
derive the model visibilities from Eq. 2.20 requires two more steps: the gains 
g can be set to unity (because model component visibilities do not require cal- 
ibration) and the obsen'ed parallel-hand visibilities r^{u,v) and r^„{u,v) are 
used instead of the true visibilities 1%.^ and respectively. The last step 

requires that the visibihties are sufficiently well calibrated in self-calibration. 
The predicted cross-hand visibilities f'^{u,v) and f^^(M, v) can now be written 
as the sum of the sub-model visibilities piJ^{ and the leakage of the observed 
parallel-hand visibihties, rotated by the parallactic angle: 



A least squares fit can now be used to solve for the D and p terms. 

The results from least squares fits are not necessarily unique, and parts of 
the true polarization of the source may be assigned to the polarization leakage 




(2.21) 
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of the feeds. Leppanen (1995) shows that the stability of LPCAL with respect 
to source structure changes is high, the rms of the Z)-term solutions being less 
than 0.002. 



Chapter 3 



The Sample 



The aim of the observations presented here is to study the sub-pc structure of 
jets in AGN, with particular focus on measuring the magnetic field strength in 
the AGN's close vicinity. The tool we chose is Faraday rotation, which occurs 
in magnetized, thermal plasmas. Evidence for Faraday rotation, however, 
requires the observation of the electric vector position angle, or EVPA, of 
linearly polarized emission at at least three, but preferably more, frequencies. 
Hence, the observations require a lot of telescope time, and we decided to look 
for polarized emission in our candidate sources at a single frequency before 
making Faraday rotation measurements. 

The following few considerations led to the selection criteria for the ob- 
serving instrument and the sources. 

• To observe scales of 1 pc or smaller in an AGN at a distance of 100 Mpc 
requires a resolution of 2 mas or better, or baseUnes of 100 MA or longer. 
The VLB A provides an angular resolution of less than 2 mas at 8.4 GHz 
and less than 1 mas at 15.4 GHz. The nominal sensitivities after a 2- 
hour integration are O.lOmJybeam^ at 8.4GHz and O.lSmJybeam^ 
at 15.4 GHz. Any birefringent effects causing depolarization due to 
wave propagation in ionized media, either in the source, in our galaxy 
or in the earth's ionosphere, decrease as v~^, i.e., are reduced by a fac- 
tor of 3.4 at 15.4 GHz compared to 8.4 GHz. From these arguments, we 
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decided to use the VLB A at 15.4 GHz as it provides an extra margin of 
resolution, sufficient sensitivity within a reasonable integration time and 
Uttle depolarization. 

Faraday rotation occurs in magnetized, thermal plasmas, but as our 
project aimed at measuring the magnetic fields and only little is known 
about AGN magnetic field strengths and structures, we confined our 
sample compilation to those sources that have clear evidence for a ther- 
mal plasma in front of the AGN core or jet. An unambiguous signpost 
of thermal plasma is free-free absorption which can be traced by ra- 
dio spectra. It causes an exponential cutoff towards low frequencies, 
and the observed spectral index, a, can be arbitrarily high. In contrast, 
synchrotron self-absorption can produce a maximum spectral index of 
a = +2.5 (e.g., Rybicki & Lightman 1979, chap. 6.8). Another ef- 
fect that can cause exponential cutoffs at low frequencies is the Razin- 
Tsytovich effect. It becomes important at frequencies below 

v« = 20^ (3.1) 

(Pacholczyk 1970, eq. 4.10), where is the particle density in cm~^ 
and B is the magnetic field strength in gauss. Typical values for 
circumnuclear absorbers are = lO^cm"^ and S = O.IG, so that 
Vr = 2 MHz. As the particle density is usually reasonably well defined 
by observations, the magnetic fields would need to be three order of 
magnitude less than currently expected to cause measurable effects 
in the GHz regime. Thus, whenever a spectral index of a > -1-2.5 is 
observed in the central regions of AGN, it is almost inevitably caused 
by free-free absorption. 

This allows us to re-state the selection criteria in the following way. 

The sources must show clear evidence for a foreground free-free ab- 
sorber. 

The sources need to be bright enough to be observed with the VLBA 
at 15.4 GHz within a few hours. Especially the absorbed parts, where 
Faraday rotation is expected to occur, need to be strong enough so that 
polarized emission of the order of 1 % can be observed. 
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Source 


Dist. / Mpc 


Scale / pemas ' 


.S(l5GHz)/mJ> 


1 cm ' 


NGC 3079 


15.0" 


0.07 


50« 




NGC 1052 


19.4* 


0.09 


489* 


3 X 10* to 10= 


NGC 4261 


35.8"^ 


0.17 


156' 


3x 10* 


Hydra A 


216"* 


1.05 


127"* 


8x 10^ 


Centaurus A 


4.2'' 


0.02 


1190^ 


> 3 X 10* 


Cygnus A 


224^ 


1.09 


600* 


10* 



Table 3.1 : Summary of source parameters. References: a - Irwin & Seaquist (1991), b - Tonry et al. (2001), 
c - Nolthenius (1993), d - Taylor (1996), e - Tonry et al. (2001), / - Owen et al. (1997), g - Krichbaum et 
al., priv. comm., h - Kameno et al. (2001), - extrapolated from Jones et al. (2001), j - interpolated from 
Tingay et al. (2001), k - Kellermann et al. (1998). 



• The sources need to be closer than 200 Mpc to achieve a linear resolution 
of less than 1 pc. 

The sample listed in Table 3.1 is, we believe, a complete list of objects 
selected from the literature that have good evidence for a parsec-scale free- 
free-absorber seen against the core or jet. We excluded 3C 84 and OQ208 
because they are unpolarized (Homan & Wardle 1999; Stanghellini et al. 1998) 
and the Compact Symmetric Object (CSO) 1946-1-708 because this class of 
objects is known to exhibit very low polarization (Pearson & Readhead 1988), 
NGC 4258 because the jet is too weak (Cecil et al. 2000), and the Seyferts 
Mrk231, Mrk348 and NGC 5506 because their parsec-scale structures do 
not show a clear continuous jet (Middelberg et al. 2004; Ulvestad et al. 
1999b). For the remaining candidates there have been no VLBI polarization 
observations yet made. 



3.1 NGC 3079 

NGC 3079 is a nearby (15.0 Mpc, de Vaucouleurs et al. 1991, 
Ho = 75kms-iMpc-i) LINER (Heckman 1980) or Seyfert 2 (Sosa-Brito 
et al. 2001) galaxy in an edge-on, dusty spiral. It does not formally satisfy 
our selection criteria as it is relatively weak (Sy < 50mJy at 15 GHz) and 
the location of the core is not known. We have observed NGC 3079 as part 
of another project and will analyse it considering the same aspects as in the 
analysis of the radio galaxies, but a separate section is dedicated to a detailed 
interpretation of the observational results of NGC 3079. 
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3.1. NGC 3079 
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Figure 3.1 : VLBI image of NGC 3079 by Trotter et al. (1998). 



Radio: VLA B-array observations at 5 GHz show a bright core and two 
remarkable radio lobes 2 kpc in length ejected from the nucleus perpendicular 
to the disc of the galaxy (Duric et al. 1983). Pilot VLBI observations by Ir- 
win & Seaquist (1988) resolved the core into two strong components, A and B, 
separated by 20.2 mas. They also found emission along a line between these 
components (component C). These three components do not align with ei- 
ther the kpc-scale lobes or the galaxy disc. The spectra of both A and B were 
found by Trotter et al. (1998) to peak between 5 GHz and 22 GHz, with spec- 
tral indices of a = 1.8 and a = 1.6 between 5 GHz and 8.4 GHz, respectively. 
Their 5 GHz image is reproduced in Fig. 3.1 for convenience. Spectral-line 
and continuum observations by Trotter et al. (1998) and Sawada-Satoh et al. 
(2000) have resolved the H2O maser emission in NGC 3079 into numerous 
aligned spots at an angle of 30° to 35° with respect to the jet axis. Trotter et al. 
(1998) suggest that the location of the AGN is on the A-B axis, where it in- 
tersects with the line of masers, although there is no radio emission from that 
location. They propose that this point is the centre of a circumnuclear disc of 
approximately 2pc diameter. Sawada-Satoh et al. (2000) identify component 
B as the core and propose that the masers are located in a circumnuclear torus. 
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Component D was found by Trotter et al. (1998) further down the A-B axis 
but was not confirmed by Sawada-Satoh et al. (2000). A and B were found by 
Sawada-Satoh et al. (2000) to separate at 0.16c, while B appeared to be sta- 
tionary with respect to the brightest maser clumps. VLA H I line observations 
of NGC 3079 and its neighbouring galaxy NGC 3073 by Irwin et al. (1987) 
revealed a tail of emission in NGC 3073 pointing away from NGC 3079. This 
tail is remarkably aligned with the axis of the radio lobes in NGC 3079 seen by 
Duric et al. (1983), and was explained to be gas stripped away by the outflow 
of NGC 3079. Irwin & Seaquist (1991) detected an H I absorbing column den- 
sity, A^m, of 10^1 cm-2 to lO^^cm'^ (assuming Jspin = 100 K) using the VLA. 
Sawada-Satoh et al. (2000) resolved the HI absorption spatially into two com- 
ponents coincident with what they identified as continuum emission from A 
and B at 1.4 GHz. From observations of the CO (J=l-0) transition, Koda et al. 
(2002) identify, among other features, a rotating disc of gas 600 pc in diameter. 
Recent VLBI observations of the red- and blueshifted H2O maser emission by 
Kondratko (2003) indicate a black hole mass of 2 x lO^M©. 

IR / Optical / UV: NGC 3079 is host to one of the clearest examples of 
a galactic superwind. Using long-split spectra, FiUppenko & Sargent (1992) 
found evidence for gas moving at speeds of up to 2000km s~^ with respect 
to the systemic velocity at distances of a few hundred pc from the galactic 
nucleus. This is considerably higher than in any other galaxies of a sample 
by Heckman et al. (1990), and therefore might not exclusively be starburst, 
but also AGN driven. This is supported by Cecil et al. (2001), who find a gas 
filament at the base of the wind that ahgns with the axis of the VLBI scale 
radio jet. 

X-ray: Cecil et al. (2002) report on X-ray observations with the Chandra 
satellite and on / band and [N II] XX 6548, 6583-i-H a line emission. They find 
striking similarities in the patterns of the X-ray and H a emission which they 
ascribe either to standoff bow shocks in the wind or to the conducting surface 
between a hot, shocked wind and cooler ISM. 

3.2 NGC1052 

NGC 1052 is a nearby (Z) = 19.4 Mpc, Tonry et al. 2001, independent of red- 
shift) elliptical LINER galaxy (Heckman 1980). 

Radio: VLBI observations at 2.3 GHz to 43.2 GHz (Kameno et al. 2001, 
Vermeulen et al. 2003) show a two-sided, continuous jet, inclined at « 90° to 
the line of sight (Vermeulen et al. 2003). VLBI observations by Kellermann 
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et al. (1998) at 15.4 GHz show a well-defined double-sided jet structure, the 
western jet of which was later found to be receding (Vermeulen et al. 2003). 
In a trichromatic VLBl observation, Kameno et al. (2001) found the spec- 
trum of the brightest component at 15.4 GHz to be a^^^'^ > +3 (Sy °< v"). 
This exceeds the theoretical maximum of -1-2.5 for optically-thin synchrotron 
self-absorption in uniform magnetic fields, and Kameno et al. (2001) con- 
clude that the absorption is due to free-free absorption. They estimate the 
path length through the absorber, L, to be « 0.7 pc, the electron density, ng, 
to be > 3.1 X lO^cm"^ and the magnetic field strength, B, to be larger than 
7 X 10~^G. These values yield a maximum RM of the order of lO^radm"-^. 
Observing NGC 1052 at seven frequencies, Vermeulen et al. (2003) find nu- 
merous components with a low-frequency cutoff which they ascribe to free- 
free absorption. They derive = lO^cm"^ if the absorber is uniform and 
has a thickness of 0.5 pc. H I absorption was first reported by Shostak et al. 
(1983) and van Gorkom et al. (1986), and more recent VLBA observations of 
the H I line by Vermeulen et al. (2003) reveal numerous velocity components 
and a highly clumpy structure of H I absorption in front of the pc and sub- 
pc scale radio continuum of both the jet and the counter-jet. They derive HI 
column densities of the order of 10^*^cm~^ to 10^^ cm^^ (T^spin = 100 K) for 
various parts of the source, and hence HI densities of lOOcm^-' to lOOOcm"^. 
A deficit in HI absorption was found to be coincident with the location of 
the strongest free-free absorption, indicating that the gas in the central regions 
is ionized by the AGN. Vermeulen et al. (2003) also report on OH absorp- 
tion and emission. The H2O maser emission found by Braatz et al. (1994) 
using the Effelsberg 100-m telescope was resolved into multiple spots along 
the counter-jet by Claussen et al. (1998). The velocity dispersion of the maser 
spots perfectly agrees with the single-dish data, which were showing an un- 
usual broad line width of 85kms~^ compared to typical hne widths of much 
less than lOkms"^ OH absorption measurements using the VLA (Omar et al. 
2002) yielded velocities equal to the centroid of the 22 GHz maser emission, 
but the connection between OH, H2O and H I is still not understood. 

IR / Optical / UV: NGC 1052 was the first LINER in which broad Ha 
emission lines were detected in polarized light (Barth et al. 1999), indicating 
that LlNERs, similar to Seyfert 2 galaxies (e.g., NGC 1068, Antonucci et al. 
1994b), harbour an obscured type 1 AGN. This is supported by Barth et al. 
(1998), who discovered that LINERs are most likely to be detected in low- 
inclination, high-extinction hosts, suggestive of considerable amounts of gas. 
There is a disagreement whether the gas in NGC 1052 is predominantly shock- 
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ionized (Sugai & Malkan 2000) or photoionized (Gabel et al. 2000). 

X-ray: Using the Chandra X-ray observatory, Kadler et al. (2004) find 
an X-ray jet and an HI column density towards the core of 10^2 cm-2 to 
10^^ cm~^, in agreement with the radio observations. Guainazzi et al. (2000) 
use the bulge B magnitude versus black hole mass relation of Magorrian et al. 
(1998) to derive a black hole mass of 10^^^ Mq. 



3.3 NGC 4261 (3C270) 

NGC 4261 is an elUptical LINER (Goudfrooij et al. 1994) galaxy located at a 
distance of 35.8Mpc (Nolthenius 1993, Hq = 75kms"i Mpc"^). 

Radio: VLBl radio continuum images at frequencies from 1.6 GHz to 
43 GHz by Jones & Wehrle (1997), Jones et al. (2000) and Jones et al. (2001) 
show a double-sided jet structure (the counter-jet pointing towards the east) 
with a clear gap 0.2 pc wide across the counter-jet at low frequencies only. The 
spectral index in the gap exceeds +3 between 4.8 GHz and 8.4 GHz, and so it is 
clearly due to a free-free absorber. Jones et al. (2001) infer «e = 3 x 10"^ cm^^ 
for L = 0.3 pc and = 10^ K. VLBI observations by van Langevelde et al. 
(2000) revealed HI absorption towards the counter-jet, 18 mas (3.1 pc) east 
of the core, with a column density of A'^hi = 2.5 x 10^' cm^^ (T^spin = 100 K), 
yielding a density of nui = 600 cm^^. No H I absorption was found towards the 
core, yielding an upper limit on the H I column density of < 2.2 x 10^*^ cm~^. 
This is much lower than that towards the counter-jet as expected, because the 
gas close to the core is ionized by UV radiation from the AGN, which ionizes 
a Stromgren sphere out to a radius of rs = 31.8pc (assuming an HI density of 
600cm~^ and a temperature of 10^ K). van Langevelde et al. (2000) ascribe 
the inconsistency between the Stromgren sphere radius and the H I absorption 
close to the core to a presumably higher Tspin. causing larger densities and 
hence smaller ionized spheres. No water maser emission was found in obser- 
vations by Braatz et al. (1996) and Henkel et al. (1998). 

IR / Optical / UV: NGC 4261 has gained some fame through the beautiful 
HST images made by Jaffe et al. (1996) showing a well-defined disc around a 
central object with a diameter of 125 pc and a thickness < 40pc. They also find 
symmetrically broadened forbidden Unes which are likely to be caused by ro- 
tation in the vicinity of a 4 x 10^ compact object in the central 0.2 pc. Fer- 
rarese et al. (1996) report on continuum and spectral HST observations from 
which they conclude that this disc has a pseudospiral structure and provides 
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the means by which angular momentum is transported outwards and material 
inwards. They also report that the disc is offset by (16ib 8)pc with respect to 
the isophotal centre and by (4 it 2)pc with respect to the nucleus, which they 
interpret as being remnant of a past merging event after which an equilibrium 
position has not yet been reached. Ferrarese et al. (1996) derive a black hole 
mass of (4.9 ± 1.0) x 10^ M©, 12 times more than Jaffe et al. (1996). 

X-ray: Chiaberge et al. (2003) observed NGC4261 with the Chandra X- 
ray observatory and detected both the jet and the counter-jet, but could not 
establish whether the emission was synchrotron emission from the jet or scat- 
tered radiation from a "misaligned" BL Lac-type AGN. They derive an HI 
column density of (6.33 ±0.78) x 10^^ cm^^. X-ray observations by Sam- 
bruna et al. (2003) using the XMM-Newton satelUte showed significant flux 
variations on timescales of 1 h. This result challenges ADAF accretion models 
because they predict the X-ray emission to come from larger volumes. 

3.4 Hydra A (3C 21 8) 

Hydra A is the second luminous radio source in the local {z < 0.1) Uni- 
verse, surpassed only by Cygnus A. It is hosted by an optically inconspicu- 
ous cD2 galaxy (Matthews et al. 1964) at a distance of 216 Mpc (Taylor 1996, 
Ho = 75kms-iMpc-i). 

Radio: Taylor (1996) observed this source with the VLBA between 
1.3 GHz and 15 GHz and found symmetric parsec-scale jets and a spectral in- 
dex of -1-0.8 between 1.3 GHz and 5 GHz towards the core and inner jets. The 
core was resolved and, for the measured brightness temperature of 10^ K and 
equipartition conditions, synchrotron self-absorption would produce a turnover 
at around 100 MHz. Free-free absorption fits the spectrum more naturally, and 
assuming that L = 15 pc, equal to the width of the absorbed region, yields 
?ie = 800cm~^. Taylor (1996) also found neutral hydrogen absorption with 
A'^Hi = 1.4 X 10^^ cm^^ towards the core and along the jets out to 30 pc, which 
he interpreted as a circumnuclear disc with thickness ~ 30pc. This result 
agrees well with X-ray observations by Sambruna et al. (2000), who discov- 
ered nuclear emission that is best fitted with a heavily absorbed power law 
with an intrinsic H I column density of Nui = 3 x lO^^cm"^ and photon index 
r 1 .7. Hydra A did not show any H2O maser emission in a survey by Henkel 
et al. (1998). 

IR / Optical / UV: Hydra A appears to have a double optical nucleus (De- 
whirst 1959). Observations by Ekers & Simkin (1983) revealed two indepen- 
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dently rotating regimes in Hydra A: a fast rotating, inner disc and an extended, 
infalling component. 

X-ray: Hydra A has an associated type II cooling-flow nebula (Heckman 
et al. 1989), characterized by high Ha and X-ray luminosities, but relatively 
weak Nil and S n and strong OI A 6300 emission lines, usually found in LIN- 
ERs. 

3.5 Centaurus A (NGC 51 28) 

The giant elliptical galaxy NGC 5128, at a distance of 4.2 Mpc (Tonry et al. 
2001), hosts Centaurus A, by far the nearest AGN and a strong radio source. 
Its proximity allowed detailed observations at almost every observable wave- 
length, showing a wealth of detail (see Israel 1998 for a comprehensive re- 
view). 

Radio: In the radio, Centaurus A covers an area of 10° x 5° in the sky 
(Cooper et al. 1965). VLA HI observations by van Gorkom et al. (1990) 
showed the neutral hydrogen to be well associated with the dust lane seen 
in optical wavelengths, and filling the gap between the kpc-scale X-ray emis- 
sion and the inner parts of the galaxy (Karovska et al. 2002). Centimetre radio 
observations show a complex jet-lobe structure from the largest (500 kpc) to 
the smallest (0.1 pc) scales. Tingay et al. (2001) found the compact VLBI 
core to be strongly absorbed below 8.4 GHz with a spectral index of +3.8 be- 
tween 2.2 GHz and 5 GHz, clearly ruling out synchrotron self-absorption, even 
though Tb > 9 X 10^ K at 8.4 GHz at the core. Absorption is seen against the 
core and not against the jet 0.1 pc away, so they let L < 0.016pc (i.e. less than 
their limit on the core diameter) and Te = lO^K and derive > 3 x 10"* cm^^ 
in the free-free absorber. Observations in the millimetre regime by Hawarden 
et al. (1993) revealed a circumnuclear disc perpendicular to the centimetre jet, 
but at an angle to the optical dust lane. They find flat spectrum at wavelengths 
as short as 800 |im and conclude that Centaurus A is a "misaligned" blazar, in 
agreement with Bailey et al. (1986), who made this assessment based on IR 
polarimetry. 

IR / Optical / UV: The appearance of Centaurus A is remarkably different 
from other elliptical galaxies due to its prominent dust lane. Moreover, the dust 
lane is oriented along the minor axis of the optical galaxy and is well inside it, 
as can be seen in the deep images of Haynes et al. (1983). 

X-ray: Observations with the Chandra X-ray satellite by Karovska et al. 
(2002) show a kpc-scale structure, suggestive of a ring seen in projection. It 
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possibly arises from infalling material where it hits cooler dust in the galaxy 
centre. They speculate that arc-shaped, axisymmetric emission on scales 
>> kpc is due to transient nuclear activity 10^ yr ago. They also find a jet- 
like feature which is coincident with the arcsecond-scale radio jet. Assuming 
optically thin thermal emission, they derive Nm = 3.9 x 10^^ cm~-^, in contrast 
to Turner et al. (1997), who found A^hi = 10^^ cm'^ using the ROSAT satellite. 

3.6 Cygnus A (3C 405) 

Cygnus A at a distance of 224 Mpc (Owen et al. 1997, Hq = 75kms"i Mpc~^) 
is the prototypical FR n radio galaxy, hosted by an elliptical cD galaxy with 
complex optical structure, located in a poor cluster. 

Radio: Observations in the radio bands display finest details from kpc to 
sub-pc scales (Krichbaum et al. 1998; Perley et al. 1984). Proper motions of 
jet components seen with VLBI were first reported by Carilli et al. (1994), who 
did not significantly detect a counter-jet. Detailed multi-frequency and multi- 
epoch VLBI imaging by Krichbaum et al. (1998) show a pronounced double- 
sided jet structure with rich internal structure and indicate that the (eastern) 
counter-jet is obscured by a circumnuclear absorber. The spectral index of 
the pc-scale counter-jet is systematically higher than the spectral index of the 
jet, and the frequency dependence of the jet-to-counter-jet ratio, R, shows the 
characteristics expected to arise from a free-free absorber rather than from 
synchrotron self-absorption. Supported by HI column densities inferred in- 
dependently from X-ray observations by Ueno et al. (1994) and from VLA 
HI absorption measurements by Conway & Blanco (1995), Krichbaum et al. 
(1998) derive «e = 10"* cm"^ (using Tg = 10"* K and L < 15pc) for the absorber. 

IR / Optical / UV: The optical emission from Cygnus A exhibits narrow 
Unes in total intensity and very broad (26000 kms"\ Ogle et al. 1997) lines 
in polarized light, supporting the presence of a hidden broad-line quasar core. 
HST near-infrared observations by Tadhunter et al. (1999) revealed an edge- 
brightened biconcical structure around a bright point source. They could not 
establish whether this structure is due to an outflow or a radiation-driven wind. 
HST near-infrared polarization measurements by Tadhunter et al. (2000) sug- 
gest an intrinsic anisotropy in the radiation fields in these cones, and a high 
degree of nuclear polarization (25 %) indicates that most of the radiation is 
re-processed in scattering and that previous works have substantially underes- 
timated the nuclear extinction. 

X-ray: From X-ray spectroscopy with the Ginga satellite, Ueno et al. 
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(1994) conclude that the nucleus in Cygnus A is absorbed by atomic hydrogen 
with a column density of A'^hi = (3.75 ±0.73) x 10^^ cm^^. This is two orders 
of magnitude more than the column density of A'^hi = (2.54 it 0.44) x 10^' cm^^ 
(Tspin = 100 K) inferred from VLA HI line observations by Conway & Blanco 

(1995) . Although Conway & Blanco (1995) cite the Ueno et al. (1994) mea- 
surements, they do not comment on the obvious discrepancy. 
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Chapter 4 



Observations 



The observations of the sample have been carried out with the Very Long 
Baseline Array, or VLB A, partly supplemented by the Effelsberg 100 m 
telescope and a single VLA antenna. The VLB A comprises 10 identical 
antennas of 25 m diameter that have been constructed exclusively for VLBI, 
in contrast to arrays composed of already existing antennas like the European 
VLBI Network. Due to the homogeneous antenna characteristics, simplified 
scheduling and user support, the VLBA has made interferometric observations 
much easier than with any other VLBI array. Some of the key advantages of 
the VLBA antennas are their frequency agility, good polarization properties 
and high slewing rates. The frequencies are changed within 10 s by moving 
the subreflector, allowing almost uninterrupted observing at many frequencies. 
The polarization leakage of the feed horns, i.e., the fraction of LCP power 
received by the RCP receiver channel and vice versa, are frequently less than 
two or three percent. Fast slewing (90° min^^ in azimuth and 30° min^^ in ele- 
vation) enables rapid switching between sources, improving phase-referenced 
observations by decreasing coherence losses due to tropospheric phase noise. 
These three advantages have contributed much to the quality of the data 
presented here. The use of the Effelsberg 100 m antenna during some of the 
observations provided the highest sensitivity on intercontinental baselines. 
This is especially valuable at low frequencies because one can match the 
synthesized beam to that at higher frequencies by using approximately scaled 
arrays, without losing too much data to tapering. A single VLA antenna 
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was also used during parts of the observations for better (m,v) coverage, 
particularly for short spacings. 



4.1 Observations 

The observations were made in five observing campaigns, the details of which 
are summarized in Table 4.1. Polarization calibrators were observed four or 
more times at each frequency, covering a range of parallactic angles larger 
than 90°. 

• NGC3079 was observed with the VLBA, a single VLA antenna (Yl) 
and Effelsberg (5.0 GHz and 15.4 GHz only) at 1.7 GHz, 5.0 GHz and 
15.4 GHz during three epochs: on November 20, 1999, March 6, 2000 
and November 30, 2000. A total bandwidth of 32 MHz was recorded 
using 2-bit samphng. This campaign was intended as a pilot monitoring 
program over a wide range of frequencies. The results - new compo- 
nents at 5.0 GHz and only few detections at 1.7 GHz - led to a different 
experiment design for the fourth epoch. 

• NGC 3079 was re-observed with the VLBA and Effelsberg on Septem- 
ber 22, 2002 at 1 .7 GHz, 2.3 GHz and 5.0 GHz. Effelsberg observed only 
at 1.7 GHz because the receiver is in the prime focus whereas all other 
receivers are in the secondary focus, and to match the 1.7 GHz synthe- 
sized beam to the beams of the higher frequencies. A total bandwidth 
of 64 MHz was recorded using 2-bit sampling. Denser frequency spac- 
ing at lower frequencies allowed us to measure the components' spectra. 
Phase-referencing was used because the source was previously found to 
be weak below 5.0 GHz and provided absolute positions for image reg- 
istration and later proper motion measurements. Polarization calibrator 
scans were inserted to caUbrate the D-terms. The angular separation to 
the phase-referencing source J0957-I-5522 is 41'. 

• NGC 1052 was observed with the VLBA on December 28, 2001 at 
13.4 GHz and 15.4 GHz. A total bandwidth of 32 MHz was recorded 
using 2-bit samphng. The experiment aimed at detecting polarized 
emission, but we were assigned enough observing time to observe at 
two frequencies with sufficient sensitivity. We decided to distribute the 
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time equally among 13.4 GHz and 15.4 GHz so that we could look for a 
change in EVPA with frequency had polarized emission been detected. 

• NGC4261, Hydra A and Centaurus A were observed with the VLB A 
on July 5, 2002, at 13.4 GHz and 15.4 GHz using several combinations 
of recording bandwidths and frequencies, but all used 2-bit sampling. 
NGC4261 was observed with 64 MHz bandwidth due to the weakness 
of the absorbed gap, and Centaurus A was strong enough to be observed 
at 13.4 GHz and 15.4 GHz to look for a possible EVPA change with fre- 
quency. Hydra A was observed at 15.4 GHz only using 32 MHz band- 
width. 

• Cygnus A was observed with the VLBA, Effelsberg and a single VLA 
antenna on July 7, 2002, at 15.4 GHz, and a total bandwidth of 64 MHz 
was recorded using 2-bit sampling. 



4.2 Calibration 

The initial data calibration was carried out in the Astronomical Image Pro- 
cessing System, or AlPS, a software package developed and maintained by the 
U.S. National Radio Astronomy Observatory. 

The visibility amplitudes were calibrated using T^ys measurements and the 
known elevation-dependent antenna gains. T^yg was measured at each antenna 
every two minutes and typically varied on timescales of hours. 

The absolute value of the complex gain of a station (i.e., how much it con- 
tributes to the visibility amplitude) is affected by fluctuations in the thresholds 
of the digitizer samplers, especially when the signal is digitized with two bits 
(four levels) per sample instead of one bit (two levels). These errors can be 
corrected using the autocorrelation data of an antenna, and the amphtudes of 
these corrections were typically less than 5 %. 

At low source elevations, the opacity of the atmosphere becomes important 
and the amplitudes are decreased. This effect needed to be considered in the 
case of Centaurus A, because its declination of —43° caused it to be below 17° 
elevation at the central five stations of the VLBA. The opacity was estimated 
and corrected as follows. System noise at a station can be expressed as a tem- 
perature equivalent (the system temperature), T^ys, which is the sum of many 
contributions from the telescope, e.g., the receiver temperature. Tree, the an- 
tenna temperature, T^i, the atmospheric temperature, T^oa, and the microwave 
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Table 4.1: Sources and observational parameters. (1) - source names, (2) - phase centre (J2000), (3) - position reference, (4) - purpose of source, (5)-(13) - 
frequencies, bandwidths and integration times used in the observations. References: 1 - Taken from Sawada-Satoh et al.'s key file for their 1996 observations, 2 - 
VLBA cahbrator survey, 3 - Gambis (1999), 4 - 43 GHz map centre in Jones et al. (2000), 5 - G. B. Taylor (priv comm.), 6 - Ma et al. (1998), 7 - T. P. Krichbaum, 
priv. comm. Superscript "EE" and "Yl" indicate the participation of Effelsberg and a single VLA antenna, respectively. 
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background temperature, Tcmb- In this simplified treatment, we refer to all of 
these contributions as Tq, except for the atmospheric noise contribution, Tatm- 
The amount of noise from the atmopshere is the product of the physical tem- 
perature of the air in front of the telescope, T^, its opacity, T, and the air mass 
in front of the telescope, which is proportional to 1/ sin(0) and is unity for an 
antenna elevation of 9 = 90° : 

Tsys = To + T^T^rA = To + (4.1) 
^ sin(0) 

Measuring T^y^ at various antenna elevations and plotting it against 
l/sin(0) yields a linear function with intercept Tsyg and slope TairTzen- To 
solve for Tzen hence requires an estimate of T^, but the atmospheric tempera- 
ture profile is a very complicated matter and would require a long treatment. 
Assuming that Tatm is 273.15 K and constant with altitude and time, we have 
derived and applied a first-order approximation of T. 

In the next step, the visibility phases were corrected for parallactic angle 
changes of the sources with respect to the antenna feed horns. This step is of 
particular importance for polarization observations as it affects the caUbration 
of the antenna feed horn polarization leakage, and also for phase referencing. 

To account for ionospheric delays, we derived corrections using the AIPS 
task TECOR from the TEC maps provided by the Center for Orbit Determi- 
nation in Europe, CODE^ We appUed corrections only to the September 2002 
observations of NGC 3079 at all frequencies. We did not apply corrections 
to any other observations because the effects are neghgible at 15.4 GHz and 
the 1999/2000 1.7 GHz observations were not phase-referenced and the source 
was too weak for self-caUbration at 1 .7 GHz. 

To allow coherent averaging over the recorded bandwidth requires one to 
remove the phase and delay offsets between the different baseband (or IF) 
channels. At the VLBA, this alignment of IPs can be done using the pulse 
caUbration system (see §2.2.3). The system provided excellent results for all 
data sets, except for Effelsberg and the VLA antenna, which do not provide 
such pulses, and self-calibration on a short scan of a strong source was used to 
measure the phase offsets. 

The final step in source-independent data calibration is the so-called fringe- 
fitting to solve for and remove residual delays and fringe-rates. The residuals 
remain due to errors in the correlator model and changing weather conditions, 
especially changes in the wet troposphere. Fringe-fitting is the first step in 

'http : / /www . aiub . unibe . ch/ ionosphere . html 
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self-calibration. The visibility data are divided into short segments. Their 

length may not exceed the time scale on which the phase rate does not change, 
typically 30s to 60s at 15.4 GHz, but this needs to be estimated from plots 
of visibility phase with time for each data set separately. These segments are 
Fourier-transformed from lag-time into delay-rate space, where a maximum is 
searched. The derived residual phase, rate and delay errors are subtracted from 
the data and can then be used to form an image by a Fourier transform from 
the (m,v) plane to the image plane. 

Fringe-fitting requires the assumption of a point source in the field centre, 
a condition which is satisfied by virtually all AGN to zeroth order. The sum of 
the visibility phases in each antenna triangle, known as closure phase (Jennison 
1958), is calculated. In the case of a point source, this sum should be zero. 
Under this assumption, deviations from the closure phase must be due to errors, 
and antenna-based errors can be derived. In a closure triangle, two baselines 
determine the phase of the third, and so there are two independent measurables 
for three antenna-based phase errors. To solve this underdetermined set of 
equations requires setting one of the antenna phases to zero and solving for 
the phase at the other two stations. This assumption destroys the absolute 
position information of the visibility phases, and phase-referencing is required 
if astrometric information is to be preserved. 

The September 2002 observations of NGC3079 were cahbrated using 
phase-referencing. Each scan of 4min (5.0 GHz) or 5.5 min (2.3 GHz and 
1.7 GHz) length on NGC3079 was sandwiched between short (Imin at all 
frequencies) scans on a strong, nearby calibrator, the fringe-fit delays, rates 
and phases of which were then interpolated onto the NGC 3079 scans. We had 
to use this observing technique because the 1999 and 2000 observations had 
shown that NGC 3079 is too weak for fringe-fitting below 5.0 GHz. In addi- 
tion, phase-referencing preserves absolute coordinates, allowing one to look 
for absolute proper motion when multiple epochs are observed. 

The data were then exported to Difmap (Shepherd 1997), a software 
package for imaging, calibration and modelling of radio interferometer data. 
Difmap implements hybrid mapping (Cornwell & Wilkinson 1981, equivalent 
to self-calibration) and easy and fast data editing. After initial editing of obvi- 
ous outliers, the data were averaged into 10 s bins to speed up computing and 
to calculate errors based on the scatter in each bin. The data were then Fourier- 
transformed into a first image, the so-called "dirty map". The dirty map is the 
sky brightness distribution convolved with the point spread function of the ar- 
ray, which is the Fourier transform of the sampling function (i.e., the tracks in 
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the (m,v) plane). 

The Difmap implementation of the "Clean" algorithm by Hogbom (1974) 
was used to obtain a first, coarse model, in which the brightness distribution 
is represented with delta components. This model is Fourier transformed back 
to the (m, v) plane to calculate the model visibilities, the phases of which are 
subtracted from the measured visibility phases. The residual antenna-based 
phase errors can then be solved for by using the closure phase condition (self- 
calibration) and a new image is formed from the measured visibilities with the 
refined calibration. Finally, the set of Clean model components is convolved 
with the "clean beam", a two-dimensional Gaussian which has been fitted to 
the Fourier transform of the sampling function. Adding the residuals from the 
Clean algorithm then yields the "clean image". 

This loop - applying Clean and subsequent refinement of visibiUty phases 
to fit the model - is known as phase self-caUbration or hybrid mapping. Be- 
cause the process is underdetermined, many different solutions exist that are 
compatible with the data and so care must be taken not to build non-existing 
components into the model. This is especially difficult with data that have 
low SNR or sparse (m,v) coverage. Difmap also implements amplitude self- 
calibration (Readhead et al. 1980), but this was only used with long (> 30min) 
solution intervals to correct slowly changing station-based gain errors. 

The calibration was considered to be complete when either the hybrid map- 
ping had converged and no further improvements were possible or when the 
rms noise in the image was within a few tens of percent of the theoretically 
calculated image noise. The former limit was mostly reached in the calibration 
processes of the strong (^max > lOOmJybeam"^) sources, the noise level in 
images of which are dominated by residual sidelobes, whereas the latter limit 
was reached in the case of NGC 3079, which has 5max < 60mJybeam~^ at all 
observed frequencies, causing the residual sidelobes to be at or below thermal 
noise. 

Polarization measurements require the calibration of the D-terms. A few 
scans on two polarization calibrators were carried out during all observing runs 
except for the 1999/2000 observations to determine them. After imaging in 
Difmap, the calibrated data were re-imported to AlPS. An image was formed 
using the same Clean loop gain, number of iterations and Clean windows as 
in Difmap. The implementation of the Leppanen (1995) calibration algorithm 
in AIPS (see §2.2.3) requires the source model to be divided into less than ten 
components, i.e., one has to define windows in each of which the Clean model 
components are thought to represent one physical source component with uni- 
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form polarization. A plot of the Clean components in a three-dimensional plot 
can help doing this, revealing which of them are strong and which represent 
only noise (Fig. 4.1). The new merged components are used by the AIPS task 
LPCAL to derive the D-term solutions which are stored in a calibration table. 
This table was copied to the target source data files, where they were used to 
form Stokes I, Q and U images, from which the polarized intensity and polar- 
ization angle maps were calculated. 




Relative RA/mas 



Figure 4.1: Plotting Clean components to reveal their relative strengths. Left panel: 5.0 GHz contour 
map of the polarization calibrator 1159+2914 observed in September 2002. Crosses indicate the locations 
of Clean components. Right Panel: The Clean components from the same image plotted to reveal their 
relative strengths. This representation reveals that almost the entire flux is concentrated in the central peak, 
and, e.g., that the group of components in the blue ellipses can be neglected in the Z)-term calibration. As a 
consequence, the calibration was done using a single component model. 

The quality of the calibration can be estimated by comparing the D-term 
solutions from the two polarization calibrators; two examples are shown in 
Fig. 4.2. The D-terms derived from the September 2002 observations of 
NGC 3079 (with the exception of the 2.3 GHz D-terms, see below) and from 
the July 2002 observations of Cygnus A were very good and consistent among 
the two calibrators. The D-terms derived from the July 2002 observations of 
NGC 4261, Hydra A and Centaurus A were a little less well constrained, but 
still of good quality. To obtain the most sensitive polarization images, how- 
ever, images were made from this epoch without the stations at Mauna Kea, 
Hancock and Saint Croix, whose D-terms showed large scatter. These images 
did not differ significantly from those that were made from all data, and the 
polarization calibration was considered to be final. We have estimated the de- 
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4.2.1 NL D-terms at 1.7 GHz on September 22, 2002. 
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4.2.2 NL D-terms at 15.4 GHz on July 5, 2002. 



Figure 4.2: North Liberty D-terms at 1.7 GHz and 15.4GHz derived from the September 22, 2002 and 
July 5, 2002 observations, respectively. Each colour represents the measurements of a D-term in a specific 
IF, derived from observations of two calibrators. The scatter is representative of the polarization cahbration 
of the aiTay on these days (i.e., there are antennas with both smaller and larger scatter). The D-terms change 
systematically between the IPs at 1.7 GHz, whereas there is no such trend at 15.4 GHz. This reflects the 
different fractional bandwidths, 0.47 % at 1.7 GHz and 0.05 % at 15.4GHz. 
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tection limit of polarized emission as follows: two polarization images of a 
calibrator were made from the same set of calibrated (m,v) data and using the 
same imaging parameters (weighting, number of iterations of Clean, etc), but 
in one case using the Z)-terms derived from the calibrator itself and in the other 
case using the Z)-terms derived from the other calibrator. The polarized inten- 
sity images were divided by the associated total intensity images to obtain two 
images of fractional polarization of the calibrator. Finally, the fractional po- 
larization images were differenced, and the average difference was measured 
over the source region. The average differences in fractional polarization were 
found to be 0.3 % for the observations of NGC 3079 at 1.7 GHz and 5.0 GHz; 
0.5 % for NGC 1052; 1 % for NGC 4261, Hydra A and Centaurus A; and 0.3 % 
for Cygnus A. 

A special problem occurred during the D-term calibration of the NGC 3079 
2.3 GHz observations. The Z)-term amplitudes were unusually large, around 
5% in both RCP and LCP, and clustered around -160° (RCP) and -100° 
(LCP) degrees, whilst the 1.7 GHz and 5.0 GHz were of the order of 1 % to 
2 %. Surprisingly, the D-terms were similar at all antennas, suggesting a con- 
nection to the antenna structure. A feature of the 2.3 GHz system at the VLBA 
is the presence of a dichroic beamsplitter in the optical path for simultaneous 
observations at 2.3 GHz and 8.4 GHz, the so-called geodetic S/X-band mode. 
The beamsplitter is a permanently installed piece of circuit board with an im- 
printed copper pattern that efficiently reflects emission at 8.4 GHz away from 
the feed onto a mirror over the 8.4 GHz feed horn. 

Applying the D-terms to the calibrator data yielded fractional polarizations 
of 0.8% and 1.7%, respectively. Applying the D-terms to the NGC 3079 
data, however, yielded a fractional polarization of about 5 %, whilst it was 
unpolarized at 1.7 GHz and 5.0 GHz. A fractional polarization of the same 
order as the D-terms is highly unlikely. The beamsplitter in the 2.3 GHz 
optical path might explain the unusual similarities of the D-terms and their 
particularly high ampUtudes. But it does not explain why the systematic effect 
has not been calibrated on NGC 3079. We have explored numerous variations 
to the calibration strategy, e.g., changing the partitioning of the calibrator 
sources into sub-models, repeating the calibration without exporting the data 
to Difmap and different data weightings, but without success. We therefore do 
not present a 2.3 GHz polarization image of NGC 3079. 

The final images are shown in Figs. 4.3 to 4.19. Figs. 4.9-4.19 have su- 
perimposed polarization ticks, except for the NGC 3079 2.3 GHz image in 



66 



4.3. The Images 



Fig. 4.10, and display the fractional polarization as colours. The lowest con- 
tours have been set to three times the image rms noise, and the polarization 
ticks have been cut off below five times the polarization image rms, except 
for the Cygnus A image, where the polarization has been cut off below seven 
times the polarization image rms to emphasise the detection. 

4.3 The Images 

4.3.1 NGC3079 

The first three epochs of NGC 3079 shown in Figs. 4.3-4.8 were observed using 
identical setups and have comparable quality. They show three components at 
5.0 GHz and two components at 15.4 GHz. These epochs have about twice the 
resolution of the fourth epoch, because they included Effelsberg at 5.0 GHz 
and 15.4 GHz. 

• 5.0 GHz Our images display three compact components and an elon- 
gated feature between the southern and the western feature. Based 
on this morphology and adopting the notation from Irwin & Seaquist 
(1988), we identify the components as A (south-east), B (west) and C 
(in between A and B). The compact component north of A was not de- 
tected by earlier 5 GHz observations, and we name it E. Component A 
is the strongest component at 5.0 GHz, causing fringe-fitting to shift it to 
the image centre. It is compact and shows only a little extension to the 
north-west, in direction of B and the jet-Uke feature C on the A-B axis. 
Component B is resolved, which is clearly visible in the second and third 
epoch, where it is resolved into three sub-components. Component E is 
also resolved, and its major axis has rotated from north-south in the first 
epoch (Fig. 4.3) to east-west in the third epoch (Fig. 4.7). The jet-like 
component C is continuous in the first and second epoch and resolved 
into sub-components in the third epoch. 

• 15.4 GHz Our images show two compact components. Based on the 
separation and relative P.A. of the two matching that of A and B at 
5.0 GHz, we identify the component towards the east with A and the 
component towards the west with B. Our identification is further sup- 
ported by 22 GHz observations by Trotter et al. (1998) who detected A 
and B with similar separation and relative P.A., and by 8.4 GHz observa- 
tions by Kondratko (2003), who noticed that E is fading towards higher 
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Source 


Epoch 


Frequency 


^maj 




P.A. 


Peak 


rm.s 


Dyn. 


P-Peak 


Prms 


Dyn. 






GHz 




ni&s 




mjybeam ' 


range 


mJybeam ' 


range 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


NGC3079 


Nov 20, 1999 


5.0 


1.92 


1.07 


-30 


20 


0.08 


270 


1 












15.4 


3.10 


1.66 


36 


45 


0.34 


140 


1 








NGC3079 


Mar 6, 2000 


5.0 


1.20 


0.85 


-29 


17.1 


0.06 


270 


1 












15.4 


2.01 


1.33 


-56 


51 


0.35 


150 


1 








NGC3079 


Nov 30, 2000 


5.0 


1.40 


0.87 


-23 


17.8 


0.06 


290 


1 












15.4 


2.25 


1.69 


-24 


51 


0.27 


210 


1 








NGC3079 


Sep 22, 2002 


1.7 


5.12 


3.49 


-22 


9.1 


0.06 


140 


1 


0.29" 


0.04 


7.3 






2.3 


6.02 


4.34 


12 


19.8 


0.13 


150 


1 












5.0 


3.86 


3.11 


-40 


22 


0.11 


200 


1 


0.43° 


0.05 


8.6 


NGC 1052 


Dec 28, 2001 


13.4 


1.49 


0.58 


-4 


344 


0.36 


960 


1 


1.36" 


0.15 


9.1 






15.4 


1.33 


0.50 


-5 


245 


0.32 


770 


1 


1.39° 


0.16 


8.7 


NGC4261 


Jul 5, 2002 




1.07 


0.60 


-6 


129 


0.18 


720 


1 


0.83 


0.09 


9.2 


Hydra A 


Jul 5, 2002 




1.58 


0.56 


3 


144 


0.29 


500 


1 


1.81° 


0.20 


9.1 


Centaurus A 


Jul 5, 2002 


13.4 


2.09 


0.47 





1080 


10.87 


100 


1 


7.98 


0.71 


11.2 






15.4 


2.15 


0.42 


1 


907 


7.46 


120 


1 


9.91 


1.05 


9.4 


Cygnus A 


Jul 7, 2002 


15.4 


0.59 


0.34 


-21 


316 


0.23 


1400:1 


1.39 


0.11 


12.6 



Table 4.2: Image parameters, a: the peak in polarized intensity is outside the image shown. 
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frequencies and probably becoming invisible at 15 GHz and that A and 
B were still strong at this frequency. Thus, our identification appears to 
be correct, although we detect neither E nor C at 15.4 GHz. A and B are 
compact in all epochs, and component B is the strongest at 15.4 GHz, 
so it is in the image centre. Fringe-fitting flags data that have too low 
SNR, and the remaining data of the first two epochs (Figs. 4.4 and 4.6) 
were distributed over a pronounced elliptical (m, v) coverage, causing an 
almost 90° change in position angle of the synthesized beam. This also 
caused the change in position angle of the major axes of A and B. No 
detections were made on baselines to Effelsberg, Mauna Kea, Hancock 
and Yl during the first epoch, to Effelsberg, Mauna Kea and St. Croix 
during the second epoch and to Brewster, Effelsberg, Hancock, Mauna 
Kea and St. Croix in the third epoch. These are predominantly long 
(> lOOMA) baseUnes, so that the resolution at 15.4 GHz during these 
epochs is lower than at 5.0 GHz. 

The flux densities of A and B changed on a significant level, especially 
between epoch two and three (Figs. 4.6 and 4.8), where A brightened 
remarkably. The ratio of the peak flux densities, 5'v^/5'v,B,is0.19±0.01 
in epoch one, 0. 1 1 ± 0.01 in epoch two and 0.28 ± 0.02 in epoch three. 
The changes have a significance of 9 a and 8 a, respectively, and cannot 
be due to amplitude calibration errors because they scale the flux density 
of both components by the same factor, leaving the ratio constant. 

In the fourth epoch of NGC 3079, Effelsberg did not observe at 2.3 GHz 
and 5.0 GHz, and so the beam size was significantly larger than in the first 
three epochs. 

• 1.7 GHz (Fig. 4.9) Two components are visible, a more compact one 
to the north-west and an extended one to the south-east. The south- 
eastern component, we named F, is the strongest at 1.7 GHz but is barely 
visible at 5.0 GHz. Due to their similar separation and relative position 
angle, E and F are easily misidentified with A and B (as unfortunately 
happened to Sawada-Satoh et al. 2000), but phase-referencing allowed 
us to identify the north-western component with component E seen at 
5.0 GHz during the same epoch. All components are unpolarized, and 
neither A, nor B nor C were detected at 1.7 GHz. 

• 2.3 GHz (Fig. 4.10) Three components are visible. Component E is the 
strongest, while F has become much fainter compared to the detection 
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at 1.7 GHz. Component A, although weak, is rehably detected, whereas 
component B is not seen. 

• 5.0 GHz (Fig. 4.11) The well-known triple structure appears to be al- 
most unchanged since the previous epoch in 2000. Component A is 
compact, component B is slightly extended in the north-south direction, 
and E is extended east-west. Due to the lower resolution, C appears to 
be connected to A and is not resolved into sub-components. F is very 
weak at 5.0 GHz. All components are unpolarized. 

Component D seen by Trotter et al. (1998) was not detected in any of our 

observations. 

The observations from 2002 and the 15.4 GHz observation from Novem- 
ber 2000 were used to compose a four-colour image of NGC 3079, in which 
red corresponds to 1.7 GHz, yellow to 2.3 GHz as, green to 5.0 GHz and blue 
to 15.4 GHz (Fig. 4.12). We prefer this somewhat arbitrary transfer function 
of frequencies to colours because a linear transfer function would map the fre- 
quencies of 1.7 GHz to 5.0 GHz into those 35% of the visible wavelengths 
(380 nm to 780 nm) that appear as red. This is due to the large gap between 
5.0 GHz and 15.4 GHz. The colour representation that we chose emphasises 
the spectral differences of the components. 

Additional Resuits from NGC 3079 

The observations presented in this thesis allow a more detailed study of 
NGC 3079 than for the other five objects. To parameterize the observations, 
integrated flux densities and component positions have been measured (Ta- 
ble 4.3. The images in this chapter display the "best" images with full resolu- 
tion obtained from each particular dataset. To minimise resolution effects and 
peak shifts due to changing imaging parameters, integrated flux densities and 
positions have been measured from uniformly weighted images tapered to a 
common resolution of 4 mas x 3 mas (not shown). 

The integrated flux densities have then been measured by summing over 
the pixels of the source region. This method is also sensitive to very weak and 
possibly extended emission which is below the lowest contour in the images 
shown. If a well-defined component was visible, positions have been measured 
by fitting a two-dimensional parabola to the brightest and the eight surround- 
ing pixels from the images tapered to 4 mas x 3 mas resolution. The fits also 
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yielded peak flux densities. This procedure was found superior to fitting Gaus- 
sians because the sources were not always well described by a Gaussian. 

In the cases of component F at 5.0 GHz and component B at 2.3 GHz, 
no component is visible in the contour plots (i.e., the peak flux density was 
below 3 a), but significant flux density was measured over the regions where 
they are visible at other frequencies. In the case of F at 5.0 GHz, the emission 
is weak and probably also resolved due to the lack of short (m,v) spacings. 
We therefore give the flux densities in Table 4.3, but do not give coordinates, 
because our procedure to measure peak flux densities did not yield meaningful 
results in these cases. We estimate the flux density errors to be of the order of 
5 %, in agreement with the VLBA's observational status summary^ and with 
what is commonly adopted in the literature for VLBA observations. Exploring 
five images from a reasonable imaging parameter space revealed a position 
uncertainty of 0.15 mas (1 a). 

The data presented in Table 4.3 can be used to derive spectral indices 
and component positions. Only the fourth epoch was phase-referenced and 
provides absolute coordinates. For the proper motion analysis, we therefore 
give the relative separations among the three components A, B and E, because 
they are well-defined and visible during all epochs. Spectral indices and rel- 
ative positions are given in Tables 4.4 and 4.5. Assuming a 1 a error of 5 % 
for the integrated flux densities yields spectral index errors of 0.09 for the 
5.0 GHz/15.4 GHz and 1.7 GHz/5.0 GHz pairs, 0.33 for the 1 .7 GHz/2.3 GHz 
pair and 0.13 for the 2.3 GHz/5.0 GHz pair. Each position has an error of 
0.15 mas, so that separations have errors of \/2 x 0.15 mas = 0.21 mas. 

4.3.2 NGC1052 

The 13.4 GHz and 15.4 GHz images of NGC1052 are shown in Figs. 4.13 
and 4.14. The excellent agreement between these two images confirms that 
the data calibration and the imaging process were correct. Furthermore, our 
15.4 GHz image nicely fits into the time series of images shown by Vermeulen 
et al. (2003). Compared to their last epoch observed in 2001.21 (ours being 
observed at 2001.99), some significant evolution is visible in the two inner and 
most compact components, whereas the outer extended components appear to 
be unchanged. 



http : / /www . aoc . nrao . edu/ vlba/obstatus/obs sum/ obs sum. html 
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• 13.4GHz (Fig. 4.13) The image shows the well-known pronounced 
double-sided jet structure, dominated by a central component with an 
extension to the west. The eastern jet is stronger than the western jet and 
widens at 5 mas distance from the strongest feature in the image cen- 
tre. It then gently narrows to 2 mas width before it fades in an extended, 
low surface-brightness feature. The western jet displays a bright, com- 
pact component with a clear extension in the direction of the jet axis 
and is connected by a very weak bridge of emission to the main com- 
ponent. There is a sharp transition to a weak, continuous feature that 
bridges 3 mas to the next, extended jet feature. Further down the west- 
em jet is a weak, isolated component. NGC 1052 is entirely unpolarized 
at 13.4 GHz. 

• 15.4 GHz (Fig. 4.14) This image differs only Httle from the 13.4 GHz 
image. Due to the slightly higher resolution, the main component is 
more clearly resolved into two components, and, probably due to a nega- 
tive spectral index, the coimection to the eastern jet has become weaker. 
The weak object at the end of the eastern jet has dimmed to the one- 
contour level and is disconnected from the jet. The connection between 
the two brightest components was not detected, now leaving a gap be- 
tween them. Compared to the 13.4 GHz image, the changes in the west- 
em jet are marginal. NGC 1052 is entirely unpolarized at 15.4 GHz. 

4.3.3 NGC 4261 

The 15.4 GHz image in Fig. 4.15 displays a double-sided, continuous struc- 
ture. The central feature gently fades into the western jet that extends out 
to 9 mas from the core. Compared to this transition, there is a rather sharp 
drop-off between the central peak and the eastem jet due to the circumnuclear 
disc. The polarization ticks 2 mas south-west of the core have a maximum 
of 0.68mJybeam^^ at a position where the total intensity is 14.4mJybeam^^ 
(4.7 % fractional polarization). The polarization was found to change consider- 
ably when the imaging parameters (weighting, tapering) were changed. They 
are probably due to residual deconvolution errors. We therefore conclude that 
NGC 4261 is entirely unpolarized. 

NGC 4261 has previously not been observed at 15.4 GHz, and we therefore 
compare our image to the 8.4 GHz and 22.2 GHz images by Jones & Wehrle 
(1997) and Jones et al. (2000). NGC 4261 is undergoing considerable spectral 
evolution between these two frequencies, and our image nicely fits inbetween. 
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At 8.4 GHz, the absorbed gap in front of the counter-jet is clearly visible. At 
15.4 GHz, the absorber is manifest by somewhat denser isophotes < 1 mas east 
of the core, and at 22.2 GHz, there is almost no indication of the absorber. The 
lengths of the jet and counter-jet are shorter at 15.4 GHz than at 8.4 GHz, and 
shorter at 22.2 GHz than at 15.4 GHz. This is due to the negative spectral 
index, the longer integration time used at 8.4 GHz and the higher sensitivity at 
this frequency. 

4.3.4 Hydra A 

The 15.4 GHz image in Fig. 4.16 shows a relatively small, double-sided core- 
jet structure, the transitions from the core into the jets being equally smooth on 
both sides. Aside from the jets, a weak component was detected 3 mas south- 
west of the core which is connected to the jets via a weak bridge of emission. 
Hydra A is entirely unpolarized. 

We see no significant structural difference between our image and the 
15.4 GHz image shown by Taylor (1996). The peak flux density in our im- 
age, (144ib7)mJybeam~\ is higher than the integrated flux density reported 
by Taylor (1996), (127 ± 6) mJy, but almost in agreement within the errors. 

4.3.5 Centaurus A 

The 13.4 GHz and 15.4 GHz images are shown in Figs. 4.17 and 4.18. The 
iterative self-calibration and cleaning was found to give multiple solutions due 
to the poor (m,v) coverage at the source declination of —43°. The similarity 
of the final images, however, indicates that the calibration of both frequencies 
has converged on the most likely source structure. 

• 13.4 GHz The image in Fig. 4.17 shows a core and a single-sided jet 
structure extending over almost 10 mas to the north-east. The flux den- 
sity drops off to 25 % of the peak flux density within the first 1 .5 mas, 
reaches a plateau and fades at 8 mas core distance. Large noise patches 
surrounding the source and a high image rms indicate significant side- 
lobes. Centaurus A has a formally significant degree of polarization of 
4.9 %, but there are also formally significant regions of polarization out- 
side the source area which obviously are due to deconvolution errors. 
These lie outside the region shown in Fig. 4.17, but similar effects can 
be seen in the 15.4 GHz image in Fig. 4.18. We therefore consider Cen- 
taurus A to be unpolarized at 13.4 GHz. 
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• 15.4 GHz The image in Fig. 4.18 is similar to the 13.4 GHz image, 
the only difference being a barely visible bend in the jet base close to 
the core. A detection of this bend has been claimed by Fujisawa et al. 
(2000) from VLBl observations at 4.8 GHz, but their (w, v) coverage was 
extremely sparse. Tingay et al. (2001) have monitored the source dur- 
ing eight epochs at frequencies of up to 22 GHz and did not detect such 
a bend in any of their epochs, and we interpret the bend claimed by 
Fujisawa et al. (2000) as an imaging artefact. We consider the polariza- 
tion seen in this image as being due to deconvolution errors, Uke in the 
13.4 GHz image. 

Our images have considerably higher resolution (2 mas x 0.5 mas) than the 
8.4 GHz images (13 mas x 3 mas) and the 22.2 GHz images (5 mas x Imas) 
presented by Tingay et al. (2001). Hence, most of the jet structure and also the 
weak, extended counter-jet feature seen by these authors is resolved out in our 
images. 

4.3.6 Cygnus A 

The 15.4 GHz image in Fig. 4. 19 shows a continuous double-sided jet structure 
spanning more than 12 mas. The core appears to be elongated and slightly 
bean-shaped. The eastern jet comprises a compact component of 1 mas extend 
and a weak, extended component 2 mas to 4 mas east of the core. The western 
jet has a brightness temperature similar to that of the compact component in 
the eastern jet, but extends further out. At a position 3.5 mas west of the core, 
the flux density sharply drops to less than 1 % of the peak flux density, and 
a diffuse, irregular structure of almost 5 mas length marks the continuation of 
the jet. Cygnus A is the only source presented here that shows significant, 
believable polarization. The central region shows 1.4mJybeam^' of polarized 
emission on top of a Stokes Ipeak of 316mJybeam^\ i.e., the region is 0.44 % 
polarized. The D-term calibration in this data set was very good, and the image 
has the highest dynamic range of all observations presented here (1400:1). 
Furthermore, the data calibration was repeated independently by a colleague 
who found very similar results (U. Bach, priv. comm.). We therefore conclude 
that Cygnus A is polarized. 

Our image displays basically the same structure as has been found by Bach 
et al. (2003), although a comparison is difficult due to the higher resolution in 
our image that was caused by the participation of Effelsberg. The faint and 
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extended jet and counter-jet structures in our image are therefore resolved and 
weaker than in the image by Bach et al. (2003). 
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full resolution 






tapered 




Epoch 


Preq. 




c 

Jmax 




5 


5max 


X 


RA 


Dec 












Component A 






Nov 20, 1999 


5.0 


27 


21 


0.78 


28 


25 


0.89 








15.4 


9.59 


6.30 


0.66 


9.03 


6.21 


0.69 






Mar 6 2000 


5.0 


26 


17 


0.65 


26 


24 


0.92 








15.4 


14.8 


6.81 


0.46 


6.75 


7.02 


1.04 






Nov "^0 2000 


5.0 


23 


17.8 


0.77 


24 


23 


0.96 








15.4 


19.0 


15.1 


0.79 


18.4 


17.5 


0.95 






Sen 22 2002 


1.7 




















2.3 


1.88 


1.24 


0.66 


2.45 


1.13 


0.46 


57.80373 


47.2436 




5.0 


25 


22 


0.88 


26 


22 


0.85 


57 80368 


47.2430 












Component B 






Nov 20, 1999 


5.0 


14.9 


6.92 


0.46 


16.9 


13.9 


0.82 








15.4 


51 


45 


0.88 


48 


46 


0.96 






Mar fi 2000 


5.0 


19.4 


5.86 


0.30 


17.9 


14.6 


0.82 








15.4 


59 


51 


0.86 


59 


54 


0.92 






Nov ^0 2000 


5.0 


19.2 


8.05 


0.42 


18.7 


14.6 


0.78 








15.4 


68 


51 


0.75 


66 


62 


0.94 






Sen 22 2002 


1.7 




















2.3 








0.64 












5.0 


22 


14.4 


0.66 


22 


14.8 


0.67 


57 801 12 

J 1 .OVA lii 


47.2594 












Component E 






Nov 20, 1999 


5.0 


4.14 


1.61 


0.39 


4.12 


3.45 


0.84 








15.4 














- 




Mar 6, 2000 


5.0 


3.28 


1.04 


0.32 


2.99 


3.22 


1.08 




- 




15.4 














- 




Nov 30, 2000 


5.0 


7.18 


1.42 


0.20 


5.14 


3.96 


0.77 




- 




15.4 


















Sen 22 2002 


1.7 


5.43 


4.23 


0.78 


3.75 


2.94 


0.78 


57 80429 


47.2633 




2.3 


21 


19.8 


0.94 


18.0 


19.8 


1.10 


57 80420 


47.2605 




5.0 


7.78 


5.41 


0.70 


111 


5.57 


0.72 


57 80413 


47.2603 












Component F 






Nov 20, 1999 


5.0 
15.4 








4.12 










Mar 6, 2000 


5.0 
15.4 








3.80 










Nov 30, 2000 


5.0 








2.88 












15.4 


















Sep 22, 2002 


1.7 


15.3 


9.71 


0.64 


17.1 


9.30 


0.54 


57.80632 


47.2521 




2.3 


7.18 


2.11 


0.29 


9.31 


2.18 


0.23 


57.80596 


47.2530 




5.0 


6.89 






8.32 











Table 4.3: NGC 3079 component data. "Tapered" flux den.sities have been measured from images tapered 
to the same resolution and convolved with a beam of 4 mas x 3 mas. "Full resolution" flux densities have 
been measured from images using ah available data. In September 2002, the tapered images differed only 
marginally from the full resolution images. Integrated flux densities, 5, are given in mJy, peak flux den- 
sities, 5niax, are given in mJybeam"', the ratio x = Smim/S is dimensionless, coordinates are seconds of 
right ascension relative to 10:01:00 and arcseconds of declination relative to 55:40:00. Only the the last 
epoch used phase-referencing and hence provides coordinates measured with respect to the phase calibrator 
J0957-H5522 at RA 09:57:38.1849709 and Dec -h55:22:57.769241. When a flux density, but no position is 
given, the emission is very extended and possibly also faint, sometimes not being visible in the images. The 
relative separations of the components measured during all four epochs are listed in Table 4.5. 
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Epoch 


SI 


A 


B 




E 


F 


1999-11-20 


„5.0 
"15.4 


-0.99 


0.92 


< 


-1.24 




2000-03-06 


„5.0 
"l5.4 


-1.19 


1.06 


< 


-0.93 




2000-11-30 


5.0 
"l5.4 


-0.25 


1.13 


< 


-1.64 




2002-09-22 


ai-'' 
"2.3 


>8.52 


>4.14 




5.11 


-1.64 




„2.3 
"5.0 


3.02 


4.52 




-1.07 


-0.28 




"5.0 


>4.56 


>4.41 




0.67 


-0.67 



Table 4.4: Spectral indices of NGC 3079 components. Limits have been calculated using three times the 
image rms. Assuming ale error of 5 % for the integrated flux densities yields spectral index errors of 0.09 
for the 5.0 GHz/15.4 GHz and 1.7 GHz/5.0 GHz pairs, 0.33 for the 1.7 GHz/2.3 GHz pair and 0.13 for the 
2.3 GHz/5.0 GHz pair. Although our observations vifere designed to avoid resolution effects, we caimot rule 
out that the spectral index of F is affected by its extend and by some of the 2.3 GHz flux density missing 
due to the lack of short baseUnes. 



Epoch 


Freq. 




A-B 






B-E 






A~E 








ARA 


ADec 





ARA 


ADec 


e 


ARA 


ADec 


© 


1999-11-20 


5.0 
15.4 


21.1 
21.1 


15.5 
16.5 


26.2 
26.8 


26.2 


1.5 


26.3 


5.1 


17.0 


17.7 


2000-03-06 


5.0 
15.4 


21.3 
21.2 


15.6 
16.3 


26.4 
26.8 


26.1 


1.4 


26.2 


4.8 


17.0 


17.7 


2000-11-30 


5.0 
15.4 


21.3 
21.3 


15.8 
16.2 


26.5 
26.8 


25.9 


1.2 


25.9 


4.6 


17.0 


17.6 




1.7 




















2002-09-22 


2.3 






















5.0 


21.7 


16.4 


27.2 


25.5 


0.9 


25.5 


3.8 


17.3 


ll.l 



Table 4.5: Relative positions of NGC 3079 components. All separations are in mas and have errors of 

0.21 mas 



Epochs 


Freq. 


A 
A© 


-B 

V 


B 
A© 


-E 

V 


A 
A© 


-E 

V 






mas 


c 


mas 


c 


mas 


c 


1999-11-20/ 


5.0 


0.2 


0.16 


0.1 


0.08 


0.0 


0.0 


2000-03-06 


15.4 


0.0 


0.0 










2000-03-06/ 


5.0 


0.1 


0.03 


0.2 


0.07 


0.1 


0.03 


2000-11-30 


15.4 


0.0 


0.0 










2000-11-30/ 


5.0 


0.7 


0.09 


0.4 


0.05 


0.1 


0.01 


2002-09-22 

















Table 4.6: Proper motions derived from separations in Table 4.5. A© has errors of 0.21 mas, and hence v 
has errors of 0.17c, 0.07c and 0.03c in the first, second and third row, respectively. 
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Figure 4.3: NGC 3079 5.0 GHz VLBA+EB+Yl image from November 20, 1999. Peak is 20mJybeam-^ 
contours are at 0.23mjy beam^' x 2^, natural weighting, beam is 1.92mas x 1.07 mas in P.A. —30°. 
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Figure 4.4: NGC 3079 15.0GHz VLBA+EB+Yl image from November 20, 1999. Peak is 
45mJybeam~\ contours are at 1.02mJybeam~^ x 2^^, natural weighting, beam is 3.10mas x 1.66mas 
inP.A. 36°. 



78 



4.3. The Images 




40 30 20 10 -10 -20 

Relative RA / mas 



Figure 4.5: NGC3()79 5.0GHz VLBA-l-EB-l-Yl image from March 6, 2000. Peak is 17.1 mJybeam-\ 
contours are at 0.19mjybeam^' x 2"^, natural weighting, beam is 1.20mas x 0.85 mas in P.A. —29°. 
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Figure 4.6: NGC3079 15.0GHz VLBA+EB+Yl image from March 6, 2000. Peak is 51mJybeam-\ 
contours are at 1.05mjybeam^' x 2^, natural weighting, beam is 2.01mas x 1.33mas in P.A. —56°. 
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Figure 4.7: NGC3079 5.0GHz VLBA-l-EB-l-Yl image from November 30, 2000. Peak is 
17.8mjybeam^' , contours are at O.lSmJybeam^' x 2^, natural weighting, beam is 1.40mas x 0.87 mas in 
P.A. -23°. 
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Figure 4.8: NGC3079 15.0 GHz VLBA-l-EB-l-Yl image from November 30, 2000. Peak is 
Slmjybeam^', contours are at 0.82mjybeam"' x 2^, natural weighting, beam is 2.25 mas x 1.69mas 
inP.A. -24°. 
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Figure 4.9: NGC3079 1.7GHz VLBA+EB image from September 2002. Peak 9.70mJybeam"', 
contours 0.19mjybeam"' x 2^, nat. weighting, beam 5.12mas x 3.49mas in P.A. —22°, pol. ticks 
Imjy beam"' mas"' where > 5a = 0.22mJy beam"' . Axes relative to RA 10 01 57.7906 Dec 
55 40 47.1788 (J2000), colours show fractional polarization. 
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Figure 4.10: NGC3079 2.3 GHz VLBA image from September 2002. Peak is 19.8mJy beam"' , contours 
are at 0.39mjybeam"' x 2^, natural weighting, beam is 6.02mas x 4.34mas in P.A. 12°. Axes relative to 
RA 10 01 57.7906 Dec 55 40 47.1788 (J2000). 
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Figure 4.12: Four-colour image of NGC 3079, composed of the observations from 2002 and tlie 15.4 GHz 
observations from November 2000. 1.7 GHz is represented as red, 2.3 GHz as yellow, 5.0GHz as green and 
15.4GHz as blue. 
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Figure 4.13: NGC1052 13.4GHz VLBA image. Peak is 344mJybeam"\ contours are at 
l.OSmJybeam"' x 2^, uniform weighting, beam is 1.49mas x 0.58mas in P.A. —4°, pel. ticks are 
SmJybeam"' mas"' where > 5(T = 0.74mJybeam"', colours show fractional polarization. 
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Figure 4.14: NGC1052 15.4 GHz VLBA image. Peak is 245 mJy beam" ' , contours are at 
0.94 mJy beam" ' x 2^, uniform weighting, beam is 1.33mas x 0.50mas in P.A. —5°, pol. ticks are 
5 mjybeam"' mas"' where > 5 C7 = O.SOmJybeam"' , colours show fractional polarization. 
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Figure 4.15: NGC4261 15.4GHz VLBA image. Peak is 129mJybeam"', contours are at 
0.54mjy beam^' x 2^, uniform weighting, beam is l.OVmas x 0.60mas in P.A. —6°, polarization ticks 
are 5mjybeam^' mas^' where > 5cT = 0.45 mjy beam" ' , colours show fractional polarization. 
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Figure 4.16: Hydra A 15.4GHz VLBA image. Peak is 144mjybeam"' , contours are at 
0.87mJybeam"' x 2^, uniform weigiiting, beam is 1.58mas x 0.56mas in P.A. 3°. The length of the po- 
larization ticks is 5mjybeam^' mas^' and start at the polarization 5 (T level of l.OOmJybeam^' . Colours 
show fractional polarization. 
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Figure 4.17: Centaurus A 13.4GHz VLBA image. Peak is l.OSJybeam"' , contours are at 
32.61 mjybeam^' x 2^^, uniform weighting, beam is 2.09mas x 0.47 mas in P.A. 0°, polarization ticks 
are 20mjybeam"' mas"' where > 5 ff = 3.55mjy beam"' , colours show fractional polarization. 
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Figure 4.18: Centaurus A 15.4GHz VLBA image. Peak is 907 mJy beam" ' , contours ai'e at 
22.38mjy beam"' x 2^^, uniform weighting, beam is 2.15mas x 0.42mas in P.A. 1°, polarization ticks 
are 20mjybeam"' mas"' where > 5 (T = 5.25mjybeam"' , colours show fractional polarization. 
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Figure 4.19: Cygnus A 15.4GHz VLBA+EB+Yl image. Peak is 316mjybeam"', contours are at 
0.69mjybeam^' x 2^, uniform weighting, beam is 0.59mas x 0.34mas in P.A. —21°, polarization ticks 
are SmJybeam^' mas"' where > 7 fj = 0.77mjybeam"' , colours show fractional polarization. 
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Chapter 5 



Discussion 



The discussion of the observational results presented in the previous chapters 
is divided into two sections: one treating the polarimetry results and one treat- 
ing the multi-epoch observations of NGC 3079, considering spectral indices, 
proper motions and general properties of nuclear radio components in Seyfert 
galaxies. 

5.1 Polarimetry Results 
5.1.1 Introduction 

The jet in M87 was the first extragalactic jet in which polarized emission 
was found (Baade 1956), and was immediately identified as synchrotron 
emission. With synchrotron emission being established as the process by 
which AGN jets radiate in the radio regime, polarization measurements 
allowed investigations of magnetic field orientations and strengths, shock 
properties, and optical depths. 

These investigations are based on the degree, orientation, location in the 
source and temporal evolution of polarized emission. A lack of polarized emis- 
sion, however, is a fundamental observational difference between the sources 
observed for this thesis and most other AGN, and it indicates a process which 
is not commonly found. It would be an unnecessary complication to postulate 
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a different intrinsic emission mechanism in polarized and unpolarized sources, 
given their similarities in so many other respects. We therefore beUeve that 
the emission process in the source is also synchrotron emission and therefore 
is intrinsically significantly polarized, and that the observed lack of polarized 
emission therefore is due to a depolarization mechanism or a blending effect. 
The mechanism(s) reprocess(es) the radiation either inside or outside the emit- 
ting regions, and we refer to them as "intrinsic" and "extrinsic" depolarization 
mechanisms, respectively. Because blending effects would also occur inside 
the sources, they are treated as an intrinsic depolarization mechanism. 

5.1 .2 Significance of Lack of Polarization 

We first demonstrate the significance of our findings (significant polarized 
emission in only one out of six sources) by comparing our results to the sam- 
ples by Kellermann et al. (1998) and Zensus et al. (2002). They have observed 
171 compact sources at 15 GHz to investigate their jet structure and kinemat- 
ics, and to derive statistical properties. The continuation of the project, known 
as the MOJAVE prgramme, has also polarimetric information on a statistically 
complete subsample (Lister 2003). This sample comprises all sources with 
declination 6 > -20°, galactic latitude \b\ > 2.5° and total 15.4 GHz VLBA 
flux density of more than 1.5 Jy. The sources are core-dominated, beamed 
quasars because one of its goals was to test Monte-Carlo simulations of beam- 
ing in a flux-limited sample of jets. Preliminary results are available on the 
MOJAVE prgramme homepage'. As of October 2003, 129 sources with po- 
larization measurements were available (Lister, priv. comm.), with a detection 
limit of 0.3 % of fractional polarization. Given this limit, the sample has 10 
unpolarized sources. In our sample, not considering Centaurus A due to the 
low image fidehty and NGC 3079 because it is clearly not a radio galaxy, three 
out of four sources are unpolarized. In a difference-of-two-proportions test 
with Yates correction for continuity, we find a probability of 0.03 % that the 
two samples have been drawn from the same parent distribution. If we require 
a detection limit of more than 0.6 % for the MOJAVE sample sources, equal to 
the averaged detection limit of our three 15 GHz observing runs in December 
2001 and July 2002, 19 of the 129 sources in the MOJAVE sample are unpo- 
larized. This subsample has a probability of 1.2 % of being drawn from the 
same parent distribution as our sources, and so the polarization properties of 
the two samples are significantly different. 

'http : / /www . physics . pur due . edu/~mlister/MO JAVE/ 
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Furthermore, the fraction of double-sided jets (sources with structure on 
more than one side of the brightest component) in our sample is four out of six. 
If we exclude NGC 3079 because it is not a radio galaxy, then the fraction is 
four out of five, and if Centaurus A is treated as a double-sided source, because 
a counter-jet is not prominent but some evidence has been found (Tingay et al. 
2001), then the fraction is five out of five. In the MOJAVE sample, the number 
of double-sided sources is six out of 129. A difference-of-two-proportions test 
yields a probability of less than 0.01 % that the two samples have been drawn 
from the same parent distribution (the probability remains less than 0.01 % 
when Centaurus A is treated as a single-sided source). 

Our sample selection criteria (free-free absorbers in front of a jet, sources 
bright enough for VLBl polarimetry at 15.4 GHz and closer than 200 Mpc), 
therefore selected sources with significantly different polarization character- 
istics compared to those satisfying the MOJAVE selection criteria (basically, 
sources brighter than 1.5 Jy at 15.4 GHz). However, the statements about the 
fraction of polarized sources and the single- and double-sidedness are not 
strong because our sample is small. Furthermore, the MOJAVE programme 
is flux-density limited at 15.4 GHz and so is biased towards highly beamed 
quasars. 

5.1.3 Intrinsic Causes 
Tangled Internal Magnetic Fields 

If the emitting regions are optically thin, the radiation that one receives is the 
superposition of radiation from along the line of sight through the source. In 
this case, for a homogeneous magnetic field and a population of relativistic 
electrons with a power-law distribution of the form N{E)dE = NqEq'' dE, the 
theoretical maximum of polarization is 

m(p) = 100 X -^i- (5.1) 

(Le Roux 1961), where m is the percentage of polarization. For p = 2, 
m{p) = 70%, and these high degrees of polarization are in fact found in ra- 
dio lobes (e.g., in Cygnus A, Carilli et al. 1989) and also in BL Lac objects 
(e.g., Gabuzda et al. 1994). In these cases, the magnetic fields must be very 
homogeneous on angular scales similar to or larger than the resolution of the 
instrument. In general, magnetic fields are not homogeneous, but are com- 
posed of a uniform component, Bq, and a random component, B,^, which varies 
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on a typical length scale, thus making "cells", in which the magnetic field is 
homogeneous. Provided that varies on scales much less than the source di- 
ameter, the intrinsic degree of polarization in each cell, mi{p), is averaged over 
the source to the observed degree of polarization, mo{p), given by 



g2 

mo{p)=mi{p) X " 

0^ (5.2) 
Energy in uniform field 



Energy in total field 

(Burn 1966). For p = 2, the fractional polarization in each cell is 70 %. To 
depolarize this below our detection threshold of 0.3 %, the magnetic field en- 
ergy would need to be dominated by the random component by a factor of 
more than 200. The jets would have to be turbulent, and very little ordering 
of the magnetic fields by the overall outward motion of the jet flow would be 
permitted. 

The fractional polarization decreases as i"^'''^, where i is the number of 
independent magnetic cells across the emitting region (Jones & Odell 1977). 
To depolarize an intrinsic fractional polarization of 70 % to 0.3 % requires 
54000 cells, and to depolarize an intrinsic fractional polarization of 10 % still 
requires 1 100 cells. Thus, the jets would need to be interspersed throughout 
with randomly oriented magnetic fields, a situation which does not agree with 
two observational facts. First, AGN jets are shaped by magnetic fields from 
sub-pc to kpc scales, requiring at least some ordering of the fields in all parts 
of the jets. Second, jets evolve along their axes, which has been observed 
in many radio jets on parsec scales. Shocks in jets compress the magnetic 
fields and so leave behind regions of ordered magnetic fields (Laing 1980). 
Evolving jets have been observed in NGC 1052 (Vermeulen et al. 2003) and 
Cygnus A (Krichbaum et al. 1998), and it is therefore unlikely that these jets 
are dominated by random magnetic fields. Furthermore, turbulent magnetic 
fields require a source of energy to sustain the disorder because ordered fields 
have lower energy. 



In the case of an optically thick source, the maximum intrinsic degree of 
polarization is 10 % to 12 % (Jones & Odell 1977), and one receives emission 
only from the surface. The same considerations as for the optically thin case 
apply, and to depolarize the radiation below the detection threshold requires 
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more than 1 100 cells across the observing beam, so that in this scenario, the 
magnetic fields inside the jets must be ordered to confine the jet, but must be 
turbulent in the radiating surface. The idea of such an outer boundary layer 
(also called shear layer) which is physically different from the jet interior 
has been raised for the first time by Owen et al. (1989) to explain their VLA 
observations of the kpc-scale jet of M 87. In the case of M 87, however, the 
boundary layer is optically thin, as indicated by up to 50 % of polarization 
(Perlman et al. 1999), and the magnetic field in the boundary layer is highly 
ordered, indicated by similar P.A. of the polarized emission over regions 
> lOOpc in size. On smaller spatial scales, evidence for boundary layers has 
been found in some BL Lac objects and FRl radio galaxies by, e.g., Gabuzda 
(2002) and Attridge et al. (1999). FRl radio galaxies and BL Lacs are thought 
to be similar objects, the difference being due to different inclination angles. 
Chiaberge et al. (2000) find that FRl radio galaxies are more luminous by 
factors of 10 to 10^ than predicted from spectral energy distributions of 
BL Lacs seen edge-on. They explain the discrepancy with a stratified jet 
structure, in which a highly relativistic spine is surrounded by slower (but 
still relativistic) layers. When viewed from smaller angles, emission from the 
beamed spine dominates the jet emission, and the jet morphology is that of a 
BL Lac. When viewed from larger angles, the beaming cone points away from 
earth, but emission from the unbeamed boundary layer remains, preventing 
the full 10 to lO'^-fold decrease of luminosity, and the morphology is that 
of an FRl. Thus, there is theoretical and observational evidence for a re- 
gion around radio jets that has different physical parameters than the jet inside. 

However, there are considerable objections to the depolarization being due 
to jet boundary layers in the sources presented here. First, even though the 
physical properties of jet boundary layers are very much different from the 
embedded spine, that does not necessarily mean that the magnetic fields are 
turbulent. Observations of the distribution of polarized emission have yielded 
evidence for the existence of boundary layers, not against it (like the above- 
mentioned observations by Perlman et al. 1999 and Attridge et al. 1999), and 
so the boundary layers must be ordered. Second, numeric models of jets by 
Aloy et al. (2000) show that the boundary layer magnetic field is well aligned 
with the jet direction, and third, Laing (1980) showed that compression of a 
volume with randomly oriented magnetic fields onto a plane almost always 
causes the radiation to be polarized by more than 10 %. Thus, only if the ini- 
tial boundary-layer random field is uncompressed and the jet does not evolve 
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and stretch the boundary layer magnetic fields, is there the possibility that the 
radiation from the boundary layer is unpolarized. Given the works mentioned 
above, we deem this situation very unlikely. 

Internal Faraday Rotation 

In the case of internal Faraday rotation, polarized emission from various depths 
along the line of sight through the source is Faraday-rotated by the source it- 
self, the degree of rotation depending on the depth of the emitting region. In- 
ternal Faraday rotation is only significant at the transition from frequencies 
where the source is optically thick to where it becomes optically thin, because 
self-absorption goes as ~ v^"', whereas Faraday rotation goes as v^^, shal- 
lower than the self-absorption. Towards low frequencies, the t = 1 surface 
(the depth in the source from which most emission is received) approaches the 
surface of the source, and only few material is left that participates in Faraday 
rotation. Towards high frequencies, the source becomes transparent, and all 
material contributes to the Faraday rotation, but the rotation measure quickly 
drops as and the effect is small. Furthermore, internal Faraday rotation 
requires a significant fraction of "cold" electrons in the jet with low Lorentz 
factors (7niin ~ 1 to 10), because the relativistic mass growth of the electrons at 
higher Lorentz factors reduces their response to the electric fields. As electrons 
in jets of radio galaxies are expected to be highly relativistic (/min ~ 100), we 
do not consider this explanation further. However, there is the possibility that 
the situation in NGC 3079 is different. Modelling the NLR line emission to 
be powered by energy and momentum transfer from radio jets, Bicknell et al. 
(1998) find that electrons in Seyfert radio components may be only mildly 
relativistic. 

Given the careful tuning of conditions that is required to explain the lack 
of polarized emission with intrinsic causes, and given how unlikely they are 
to be fulfilled, we prefer to reject those explanations, and the depolarization 
must have a different origin. 

Both intrinsic and extrinsic depolarization mechanisms impose an upper 
limit on the typical length scale on which the magnetic field changes direction. 
Assuming that this scale is similar for all sources, the best limit is given by 
NGC 1052 at a distance of 19.4 Mpc, since it is the closest of the double-sided 
jet sources in which a solid non-detection was made. In NGC 1052, our highest 
resolution image at 15.4 GHz resolved 0.5 mas, corresponding to 0.047 pc in 
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the source. Within a volume of 0.047 pc on a side, i = 54000 independent 

cells are required to depolarize an intrinsic degree of 70 % polarization to the 
observed upper limit of 0.3 %. Thus, the typical length scale of a cell is < 
54000 of the observing beam, or 0.0012 pc. In the case of the source 
being optically thick, the typical cell size must still be < 1/10 of the beam, or 
0.0047 pc. 



5.1.4 External Causes 
Bandwidth Depolarization 

At each wavelength, the electric vector position angle of the waves after pass- 
ing through the medium is different, and integrating over the observing band- 
width yields an apparent net polarization that is lower than the initial state. 
The depolarization across the observing bandwidth can be calculated as fol- 
lows. Consider an ensemble of waves with wavelengths A in the observed 
bandwidth AA, so that Ao < A < -I- AA. Let us assume that all waves have 
the same electric vector position angle. If the waves pass through a magnetized 
plasma, the electric vectors E{X) are rotated by Faraday rotation. 

Assuming that a homogeneously magnetized slab of plasma has a rotation 
measure RM, the change of EVPA with wavelength is given by 



e{Xi)=RMxXf + do (5.3) 
^A0(Ai,Ao)=/?M(Af-Ao2), 

where do is the EVPA of the waves before passing through the medium and 
( A ) is the wavelength-dependent EVPA after passing through the medium. 
With respect to the reference wavelength Aq, the electric vectors £(A) are ro- 
tated by A0 (Ai , Aq): 

£(A)=£(Ao)e'^^(^'^) 

^ ,,2 ,2^ (5-4) 

= £(Ao)e'«^(^-^). 

Integrating all vectors in the observing band yields 
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Al — Ao JXi 

= r E{A^)e"'''^''-'^UX (5.5) 

Ai — Ao JXi 

This integral can be solved numerically (Fig. 5.1), and one finds that for 
the radiation to be depolarized to 50% of initial value requires RM = 2.2 x 
10^ in a 32 MHz wide band at 15.0 GHz. The largest values of RM that have 
been measured in quasar cores are a few thousands (e.g., Zavala & Taylor 
2003). Assuming a rotation measure of RM = lO^radm^^ causes 0.00012% 
depolarization. Bandwidth depolarization is unlikely to be important because 
unprecedently high and uniform (to rotate the electric vector position angle 
in one direction only) rotation measures would be needed to depolarize the 
sources, requiring very uniform conditions in the Faraday screen. 




Figure 5.1 : Bandwidth depolarization in the 15 GHz band. 
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Beam Depolarization 

If the rotation measure in a foreground Faraday screen changes on angular 
scales much smaller than the observing beam, regions with similar degrees 
of polarization but opposite signs average out and the observed degree of po- 
larization is decreased. Thus, spatially highly variable RMs can in principle 
depolarize the source. Estimates of RM and the length scale on which the 
magnetic field needs to be tangled or the electron density changes follow simi- 
lar approaches as in the previous sections. The electric vectors must be rotated 
by 180° on average to remove any residual polarized emission, requiring the 
rms of the variations in the Faraday screen, Orm, to be 



e{l) xRM 

^RM = d{X)/X' 

;rrad „ , _j 

^^^= (c/15.4GHz)2 ='''Q"^"^" 

or larger at 15.4 GHz. If the RM changes are caused only by density fluctua- 
tions in the medium and L constant), then the rms of the fluctuations is of 
the order of 10"* cm^-'. For NGC 3079, the constraints on the Faraday rotation 
are weaker because the highest frequency at which no polarized emission has 
been observed is 5.0 GHz, and so the rms of RM must be 872radm~-^ or more. 
The estimate of the rotation measure cell size was calculated in §5. 1 .3, yielding 
upper Umits of 0.0012 pc to 0.0047 pc for NGC 1052, the closest double-sided 
source. 



Free-Free Absorption 

Conditions that cause such high rotation measures can also produce strong 
free-free absorption, which is seen only towards comparatively small parts of 
the jets, due to which these sources were selected, while most of the jet parts 
appear to be unabsorbed (with the exception of NGC 1052, where free-free 
absorption has been detected towards larger parts of the jets (Kameno et al. 
2001; Vermeulen et al. 2003). The rotation measure depends linearly on the 
electron density, whereas free-free absorption depends on the electron density 
squared: 



RMoc / n.dl 



(5.7) 
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Tfo< J nidi. (5.8) 

Thus Ty decreases more rapidly than RM with decreasing electron density. 
Furthermore, a turbulent magnetized plasma can depolarize 

the radio emission without causing too much free-free absorption if the 
absorber is either extended (so the density is low) or hot, or both. 

We can test the feasibility of such a gas by estimating the free-free absorp- 
tion and Faraday rotation caused by conditions typically inferred for circum- 
nuclear regions in AGN. The free-free optical depth (Eq. 1.12) depends on the 
temperature, T, the observing frequency, V, and the density of positive and 
negative charges, A'^^. and A^_, respectively. We assume that any circumnuclear 
gas is dominated by hydrogen, and that therefore, = N- = rie,, where 
is the electron density. We further assume T = lO^K, = 3 x lO^cm"^ and 
L = 0.5 pc. The temperature of gas exposed to ionizing UV radiation is ther- 
mostatically regulated to a temperature of 10"* K. This is because the most ef- 
ficient cooling is by low excitation levels (1.9 eV to 3.3 eV) in metals (Spitzer 
1978). The probability of electrons in a Maxwellian distribution having the 
average energy of 2.6 eV increases by a factor of 32, when the temperature of 
the gas increases from 5000 K to 15000 K, and the gas therefore cools much 
more efficiently when a temperature of more than 10000 K is reached. The 
electron density and the path length are typical values for the well-studied ab- 
sorbed gaps in NGC 1052 (wg > 3.1 x 10^ cm^^ and L = O.Vpc, Kameno et al. 
2001, ne = lO^cm'3 and L = 0.5pc, Vermeulen et al. 2003) and NGC 4261 
(ne > 3 X 10"* cm"^ and L = 0.3 pc, Jones et al. 2001), and Eq. 1.12 yields 
Ty = 0.95 at 15.4 GHz. The same conditions would give rise to a rotation mea- 
sure of/?M= 1.2x 10^'radm^,ifS|| = 0.1 mG (comparable to what Lobanov 
(1998) derived for Cygnus A). 

Whilst the abovementioned conditions are a good fit for the well-defined 
free-free absorbers in front of the jets, the conditions along the unabsorbed but 
still depolarized parts of the jets must differ to produce an optical depth close 
to zero and maintain the rotation measure still reasonably high. To satisfy 
RM = lO^radm"-^, the electron density needs to be a factor 120 lower than in 
the example above, and the optical depth would then be Ty = 0.00. 

We can now summarize the absorber properties as follows. To yield 
RM = lO'^radm"^ within 5pc of the AGN (so that the jets fit into the 
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absorber) with a tangled magnetic field of an average strength of = 0.1 mG 
requires an electron density of only 25 cm^^. An important consequence 
is that the depolarization cannot be due variable Faraday rotation caused 
by density fluctuations. The fluctuations would need to be of the order of 
lO^cm"^ as shown above, and assuming a Gaussian distribution of these 
fluctuations, the average density then is 5 x 10^ cm~^, yielding too high optical 
depths. 



5.1 .5 Origin of the Faraday Screen 

Let us assume that a uniform absorber is present within several parsecs of the 
AGN. The absorber must be ionized on each line of sight to parts of the jets 
which were found to be unpolarized, i.e., the ionized regions in the absorber 
must be larger than the jet lengths. The radius of a sphere which can be com- 
pletely ionized by a source of ionizing photons with A < 91.1 nm is called 
the Stromgren sphere radius. It is given by, e.g., Osterbrock (1989) (in units 
convenient for AGN astrophysics): 

Ruu = m.6T^'^N^I^n-^'\ (5.9) 

Te is the electron temperature in K, the number of UV photons with 
X < 91.1 nm emitted per second, and is the electron density in cm~^. We 
estimate the number of ionizing photons as follows. The spectral energy dis- 
tribution of most AGN can be approximated by a power-law over the entire 
spectrum. The flux density at a certain frequency is then given by 

5v=5o(^)" (5.10) 
The number of photons required to produce this flux density is 

^v = ^ = r^f-V = 4;v«-i (5.11) 

hv hv \VoJ V^h 
Integrating over a given frequency range yields 
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For example, NGC4261 has 5o = 5(8.2 x lO^^Hz) = 3.1 x 
10-28wm-2Hz-i (de Vaucouleurs et al. 1991), and 5(5.3 x lOi'^Hz) = 
2.5 X lO^^^Wm^^H^-i (Fabbiano et al. 1992). The optical-to-7 spectral 
index therefore is a = —1.81. Integrating from the ionization limit of 
13.6eV = 3.3 x lO^^Hz to lO^^Hz yields 21000m-2s-i. At the distance of 
NGC4261 (35.8 Mpc, Nolthenius 1993), this corresponds to a total production 
rate of ionizing photons of 3.2 x 10^^ s^^ /?hii then follows to be 260 pc, 
when We = 25cm~^ and = 10^ K are assumed. Thus, the Stromgren sphere 
radius exceeds the required minimum of ~ lOpc imposed by the longest 
unpolarized jets in Cygnus A and Hydra A by two orders of magnitude. A 
Faraday screen which is not seen as a free-free absorber therefore can be 
ionized by the AGN. Note that UV radiation is easily absorbed in galaxies and 
that the intrinsic UV fluxes which contribute the bulk of the ionizing photons 
might be considerably higher. Similar estimates for the other sources are listed 
in Table 5.1. 

We conclude that the non-detection of polarized emission is due to a fore- 
ground screen in which non-uniform Faraday rotation causes beam depolar- 
ization. The low electron density in the screen imposed by constraints on the 
optical depth requires the changes in Faraday rotation to be caused by turbulent 
magnetic fields and not by density fluctuations. Such a low-density screens can 
be completely ionized by the AGN radiation. 



Source 


Vl 






Ref. 




V2 


5(V2) 




Ref. 


'•s 




Hz 


Wm-^Hz-' 






Hz 


Wm-^Hz-l 




pc 


NGC 3079 


8.2 X 10'** 


4.3 X 10- 


-28 


1 


5.3 


X 10" 


1.1 X 10- 


-ii 


2 


144 


NGC 1052 


8.2 X 10'** 


3.4 X Kr 


-28 


1 


5.3 


X 10'^ 


1.0 X ur 


-33 


2 


155 


NGC 4261 


8.2 X lO''* 


3.1 X 10" 


-28 


1 


5.3 


X lO" 


2.5 X 10- 


-33 


2 


264 


HydA 


9.6 X 10* 


6.5 X 10- 


-25 


3 


3.3 


X 10" 


1.4 X 10- 


-31 


4 


373 


Cen A 


8.4 X W 


3.9 X 10- 


-24 


5 


6.0 


xlO's 


1.3 X 10- 


-32 


6 


83 


CygA 


1.4 X 10' 


1.6 X 10- 


-23 


7 


1.8 


X 


4.0 X 10- 


-29 


8 


436 



Table 5.1 : Stromgren sphere radii for the observed sources. References: 1 - de Vaucouleurs et al. (1991), 
2 - Fabbiano et al. (1992), 3 - Kiihr et al. (1981), 4 - Brinkmann et al. (1994), 5 - Jones & McAdam (1992), 
6 - Steinle et al. (1998), 7 - Kellermann et al. (1969), 8 - Lebofsky (1981) 



5.1.6 Possible Absorbers 

In this section, we investigate the absorption and depolarization to be expected 
from various regions of ionized gas surrounding AGN. 



5.1. Polarimetry Results 



101 



Broad line region clouds Forbidden lines have almost never been ob- 
served in BLR clouds, and so they must have densities of at least lO^cm"^, 
because the forbidden transitions are coUisionally de-excited before they can 
radiate (e.g., Osterbrock 1989). Adopting = lO^cm"^, Sy = 0.1 mG and 
L = 0.5 pc, the rotation measure is 4 x lO^radm"^, sufficient to depolarize 
a background radio source, but the optical depth of the absorber would be 
Ty = 1.1 X 10^ and one would not see the background source. Also, BLR 
clouds are found in only the innermost tens of light-days of the AGN and have 
filling factors of the order of a few percent, whereas the depolarization in the 
jets extends at least out to 1 pc and is continuous. The depolarization therefore 
cannot be due to BLR clouds. 

Narrow line region clouds NLR clouds have much lower densities than 
BLR clouds, typically between lO'^cm"^ and lO^cm"^. They can extend to 
scales of kpc and have volume filling factors of 10~^ to 10~^. The free-free 
optical depths at 15.4 GHz following from these parameters cover the range 
of 0.01 to 100, and at 5.0 GHz cover 0.11 to 1100. In NGC 1052 free-free 
absorption has been detected towards the radio jets (Kameno et al. 2001; Ver- 
meulen et al. 2003). Adopting T = 1 for the highest observed turnover fre- 
quency of 43 GHz, Vermeulen et al. (2003) derive = 10^ cm^^ if the gas is 
uniformly distributed over 0.5 pc. The optical depth decreases with turnover 
frequency towards larger distances from the core. If a magnetic field of 10 /iG 
is adopted, comparable to interstellar magnetic fields, we find a rotation mea- 
sure of 4 X lO^radm^^. We conclude that NLR gas cannot be excluded as the 
cause of depolarization in NGC 1052. The situation is similar in NGC 3079. 
Our observations yield good evidence for a foreground free-free absorber (see 
§5.2.1) with optical depths of order unity at 5.0 GHz and hence may be due to 
NLR gas. The expected rotation measure is of the order of a few thousand and 
so exceeds the 870radm~^ required to depolarize NGC 3079 at 5.0 GHz. 

No free-free absorption has been detected along the jets in the other four 
sources, and so their NLR gas is unlikely to be the cause of depolarization. 

H II Regions and Interstellar Gas The line of sight to the radio jets in- 
tersects H II regions and interstellar gas in the sources. The H II clouds are 
ionized by UV radiation from stars, and so may contribute to the depolariza- 
tion of background radio emission. Although the interstellar medium, with 
densities of 0.1 to 10 atoms per cubic centimetre, is mostly neutral, it has a 
spatially variably ionized component. Both H II regions and interstellar gas 
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therefore contribute to depolarization. However, one would expect the rotation 
measures to be similar to those on lines of sight through our Galaxy. Galac- 
tic rotation measures have been well determined using pulsar observations, 
(e.g., Weisberg et al. 2004, Mitra et al. 2003 and Han et al. 1999), and these 
authors find rotation measures of mostly < lOOradm"-^, even for galactic lati- 
tudes below 10° and pulsar distances of several kiloparsecs. Only few pulsars 
exceed lOOradm^^, but no pulsar shows > lOOOradm^^. Only the highest ro- 
tation measures therefore are of the order of magnitude required to depolarize 
NGC 3079. The argument is strengthened considering that NGC 3079 shows 
more star forming activity than our galaxy and so the interstellar gas has a 
higher ionized fraction. We therefore cannot rule out that the radio emission 
from NGC 3079 is depolarized by H II regions and interstellar gas. The inter- 
stellar medium is unlikely to be important in the other five objects which have 
elliptical host galaxies. First, ellipticals generally show fewer star forming ac- 
tivity than spirals and second, the interstellar gas in these objects has a lower 
density. 



Spherical accretion onto AGN The origin of the material which settles 
around AGN in accretion discs is largely unknown. We present here a simple 
model with a spherical infall of gas towards the AGN (Bondi accretion). The 
gas is ionized by the UV radiation from the accretion disc and therefore causes 
Faraday rotation. 

The mass distribution in the sphere can be approximated with a rotation- 
ally symmetric density profile, and the expected rotation measure can then be 
obtained by integrating over the density profile. 

Let us assume that gas with a density p (in kgm~^) is freely falling radi- 
ally inwards in the direction of the AGN, so that its density is a function of 
the radius R. The free-fall velocity is given by vg = a(2GM,//?)^/-^, in which 
M, is the black hole mass in kg, R is the radius in m, the dimensionless vis- 
cosity parameter k is commonly set to 0.05 and G is the universal gravitation 
constant. The density profile is then given by 

^ 4nR^vs 47rJc(2GM.)i/2 

where M is the accretion rate in kgs^^ . We estimate the magnetic field in the 
freely falling gas by means of equipartition: 
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1 1 M/J-3/2 



4;r 



where denotes the permeability of free space. 

For the rotation measure, we need the particle density, which can be de- 
rived by dividing the density p (R) by 2{xmu, where mu is the proton mass and 
the factor 2/i with /i = 0.622 accounts for the fraction of helium. 

p(R) M 
* 2ixmn Snmn7tK{2GM,y/^ 
Inserting both B and into Eq. 1.10 now yields 



RM = 2.63 X 10"^^ 



/ 



M/?-3/2 

Sixmn7iK(2GM.y/^ 

1/2 



X ( iE^foM!^^ MV2m1/4/?-5/4^^ 

2.63 X lO-l^/i,)/^ M3/2 / _7/4 „7/4\ 
" Jcl/2(2G)l/4l6^mH;r3/2 ^ r-"!" -^max J 

= 4..9x,0-*;(«-'^-.i/^). 

We exemplify the calculation with data from NGC4261. The black hole 
mass was estimated by Ferrarese et al. (1996) to be 4.9 x 10** Mq = 9.8 x 
10^^ kg, and assuming that the accretion rate is 1 % of the Eddington accretion 
rate (see below) one gets O.OlMQyr"^ = 6.9 x 10-^*'kgs~^ Assuming that the 
infall starts at lOpc = 3.09 x 10^^ m (the results do not change much when 
^max > lOpc) and stops at 0.2pc = 6.17 x 10^^ m (where the disc is seen in 
NGC4261), one gets RM = 3.1 x lO^^radm^^ According to Eq. 5.15, the 
electron densities at the inner and outer edges are ne(0-2pc) = 3.0 x 10^ m~^ 
and ne(lOpc) = 8.5 x 10^ m~^, corresponding to 3.0 x 10^ cm~^ and 85cm~^, 
respectively. 

We also need to verify that the optical depth is not too high. Equation 1.12, 
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converted from units of cm ^ for and pc for L to m ^ and m, respectively, 
takes on the form 

= 2.67x lO-^^r-i-^^v-^ i JnldL (5.17) 

Inserting Eq. 5.15 into this equation, substituting the path length, L, by the 
radius, R, and integrating over R from /?nim to R^ax yields 



= 2.67 X lO"^*^ r-^-^^ v"^-^ 



MR-^- 




\S^lmu7tK{2GM.y/^ 

2.67 X lO-^'^ y-^-^M^ 
7Z X / ■'i//? (5.18) 



1 2.67 X 10-3*' r-i-35v-2-iM2 

2 ^ (8AimH7rK-)2 2GM. 



V max "min/ 



= 5.85xl033r-1.35^-2.1^^(^-^2_^-2^) 

Converting to convenient units, the equation becomes 

1.23 X l0i6r-i-35 X (/?^„-/?-L) • (5.19) 

In this equation, T is in K, v is in GHz, M and M, are in M0 yr~^ and Mq, 
respectively, and R is in pc. For NGC4261, the free-free optical depth then 
is 0.00094, so that radiation at 15.4 GHz can pass through almost unaffected 
by free-free absorption. Thus, a spherical accretion halo around the AGN can 
depolarize the synchrotron radiation without being seen as an absorber. As this 
is a promising approach, we calculate RM and T for each galaxy in the sample. 
Accretion rates are difficult to determine, but they can be estimated using the 
black hole masses as follows. The luminosity at which the force generated 
by radiation pressure balances the gravitational force is called the Eddington 
luminosity. It can be derived by equating the force due to radiation pressure 
and the gravitational force 
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^Edd<7T _ GM,m 

AKcGM.m "^^-^^^ 

A-Edd = , 

(7t 

where Lead is the luminosity in Js^^, c is the speed of light in ms~^ G is 
the gravitational constant, M, is the black hole mass in kg, m is the mass of 
a particle in kg (here, the proton mass) and is the Thompson scattering 
cross-section in m^. 

The Eddington accretion rate can then be estimated as the mass equivalent 
of the Eddington luminosity to place an upper limit on the accretion rate. 

M=^. (5.21) 

The fraction rj describes the actual accretion rate in terms of the Edding- 
ton accretion rate, but unfortunately is only poorly constrained by models. 
Brunner et al. (1997) have modelled lUE and ROSAT data from observations 
of 3 1 radio-quiet quasars as thermal emission from a thin accretion disc. The 
median black hole mass in their sample is 6.1 x 10^ (1.2 x 10^^ kg), and 
they derive accretion rates of the order of 10 % to 30 % of the Eddington 
accretion rate, with a median of 12%. O'Dowd et al. (2002) studied the 
relation between AGN power and host galaxy using 40 BL Lac objects and 22 
radio-loud quasars. They find black hole masses between 10^ M© and lO^^M© 
and derive accretion rates of (2 x lO^'^ to l)xLEdd- We use accretion rates 
of 1 % of the Eddington accretion rate because this is the average accretion 
rate found by O'Dowd et al. (2002) for the BL Lac objects. BL Lacs are 
considered "misaligned" FRl radio galaxies, the class to which NGC4261, 
Hydra A and probably also Centaurus A belong. NGC3079 and NGC 1052 
are classified as Seyferts or LlNERs, and it is not clear what the accretion rates 
in these types of objects are. NGC 1052, however, displays a double-sided 
radio jet which marks considerable nuclear activity. Unfortunately, the black 
hole mass is only poorly constrained but the suggested range overlaps with 
the black hole masses suggested fo FR 1 galaxies, and for now, we assume the 
same accretion rate. Cygnus A is a FR 11 galaxy and so is probably a powerful 
quasar viewed edge-on. O'Dowd et al. (2002) find higher accretion rates 
for quasars than for BL Lacs, and so the accretion rate and hence rotation 
measure we find in Cygnus A are probably only lower limits. The model 
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details for each source are listed in Table 5.2. The rotation measures are 
higher than the required 8290 radm^^ in each of the sources, and is higher 
than the required 872radm^^ (due to the lower highest frequency of 5.0 GHz 
at wliich no polarized emission was found) in NGC 3079. The optical depths 
in all sources are small enough to be neglected. We conclude that the spherical 
accretion model, if correct to an order of magnitude, can significantly reduce 
the observed degree of polarization. 

Two arguments need to be considered. First, the rotation measures 
we derive critically depend on the value of Rmia- We have adopted the 
scales on which the well-defined free-free absorbed gaps were observed in 

NGC 1052 and NGC 4261, indicating where a different source geometry sets 
in. The depolarization screens in NGC 4261, Hydra A and Cygnus A pro- 
vide enough Faraday rotation as long as Rmin < 3 pc, but the rotation measures 
in NGC 3079, NGC 1052 and Centaurus A then drop to lOOradm'^ and the 
spherical accretion model fails. We are, however, convinced that adopting 
^min < Ipc is justified, and the model then predicts enough Faraday rotation. 

Second, the accretion rates in NGC 3079 and NGC 1052 might be signif- 
icantly lower than the 1 % of the Eddington accretion rate. The bolometric 
luminosities of these two AGN are at least two orders of magnitude lower 
than in FRl radio galaxies, suggesting that the accretion rate also is lower. 
With M = 0.0001 X LEdd» the model predicts rotation measures of less than 
50radm^^ for both objects. We conclude that our model cannot explain the 
depolarization in NGC 3079 and NGC 1052. 



Source 


M. 


i-Edd 


M 


RM 


Tff 






IO^Mq 


lO^^w 




radm~^ 






NGC 3079 


2" 


2.5 


4.44 


3.2 X 10^ 


3.9 X 10- 


-6 


NGC 1052 


10 to 1000* 


13 to 1300 


22 to 2200 


2.4 X W* to 7.6 X 10* 


1.9 X 10" 


-5±1 


NGC 4261 


490 ±100'' 


620 ±130 


1090 


(3.1 ±0.4) X 10'' 


9.5 X 10" 


-4 


Hyd A 


490'' 


620 


1090 


3.1 X lO* 


9.5 X 10" 


-4 


Cen A 


10" 


13 


22 


2.4 X lO'' 


1.9 X 10" 


-5 


CygA 


nof 


220 


377 


8.3 X lO'* 


3.3 X 10" 


-4 



Table 5.2: Rotation measures for the sample sources. The accretion rate has been assumed to be 1 % of 
the Eddington accretion rate, but might be significantly lower in NGC 3079 and NGC 1052. References: a - 
Kondratko (2003), b - Guainazzi et al. (2000), c - Ferrarese et al. (1996), d - Woo & Urry (2002), e - Israel 
(1998), /-Conway (1999) 
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King model In an alternative model, we assume that the gas distribution 
follows the model by (King 1972). The model was invented to describe the 
density in stellar and galaxy clusters. If we assume that stars are present even 
in the very vicinity of the AGN, the gas kinematics are dominated by the stars 
due to their comparatively large mass. However, the black hole begins to 
dominate the gas kinematics at radii 1 < r < lOpc (depending on its mass), 
and as we are observing radio emission on scales of the order of 1 pc, we 
actually should not apply the King (1972) model. We nevertheless present the 
calculation to estimate the order of magnitude. 

The gas density is given by the King model as 



Here, no and rc are normalization constants. We require ?ie = 10^ m ^ = 
lO^cm-^ at r = 6.17 X lO^^m = 0.2pc, and so no = 2.83 x 10^ m'^ and = 
6.17X lO^^m. 

Inverting Eq. 5.15 yields the density p 




(5.22) 



(5.23) 



which we use to derive the equipartition magnetic field: 




(5.24) 




Inserting this and the electron density into Eq. 1.10 yields 
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RM = 2.64 10"^^ X 



no 



1 + 



-3/2 



l+(- 



-3/4 



= ax [ 
Jr 

= a I 7?niax iPl 



dr 



1 + 



1+1- 



-9/4 



19 3 Rl 



2' 4' 2' 



-9/4 



dr 



19 3 Rt 



2' 4' 2' 



(5.25) 

where a = 2.64 10-i^(4/XoM'"h«o)'''^ = 2-872 x 10"!^ and 2Fi[k, I; m; n] is 
the Hypergeometric2Fl function 



2Fi[k, I; m; n] = Y,ik)i{l)i/{m)in' 

!=0 



in which brackets denote 



(5.26) 



{x)o = 1 

{x)i =x{x+l){x + 2)...{x + i) 



(5.27) 



The integral can be evaluated numerically and yields a rotation measure 
of 1.9radm^^, much too low to depolarize the radiation. Even choosing no 
to be one order of magnitude higher than in the previous case yields a rotation 

—9 3 /2 

measure of only 59radm , because Eq. 5.25 goes with «q . We conclude that 
the King (1972) model underestimates the Faraday rotation by several orders 
of magnitude. 
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Isothermal gas distribution Let us assume that the circumnuclear gas 
is isothermal and has a rotationally symmetric power-law density profile with 
index — i. Its pressure can then be written as 

p{r)=por'kT (5.28) 

The gas pressure has to balance the gravitational potential of the black 
hole: 



Por dr 

_GM^J_ , ^3 29) 

GM _ kT 

The two sides of the equation can be equal only at a single point because 
of the different dependencies on the radius. An isothermal gas is therefore not 
able to exist in the gravitational field of a black hole. 



5.2 Additional Results from NGC 3079 

This section discusses more aspects of the multi-epoch, multi-frequency ob- 
servations of NGC 3079. As an orientation, we have sketched NGC 3079 in 
Fig. 5.2 together with the frequencies at which the various components have 
been detected. 



5.2.1 Spectra 

The flux densities of components A, B, E and F measured in November 2000 
and September 2002 are shown in Fig. 5.3, displaying three unusual features. 
First, there are unusually inverted spectra, e.g., a^f^{A) = 3.02 and ajs (£") = 
5.11. Second, there is the fairly steep spectrum of F, a^Q^F) = —1.64. Third, 
there is a striking correlation between the peak frequencies and the component 
positions in the source, namely that that the turnover frequencies are increasing 
with decreasing right ascension. This is illustrated in Fig. 5.3. We now address 
these three points in turn. 



110 



5.2. Additional Results from NGC 3079 



E B 
1.7-5.0 GHz 5.0-15.4 GHz 




A 

2.3-15 GHz 



Figure 5.2: Sketch of all components seen in NGC 3079, showing the frequencies at which each compo- 
nent was detected. 



Unusually Inverted spectra First we demonstrate that the emission from 
NGC 3079 is synchrotron emission. The brightness temperatures Tb (the phys- 
ical temperature a blackbody would need to have to produce the observed flux 
density) can be calculated by inverting the Rayleigh- Jeans approximation of 
Planck's law, yielding 

(5.30) 

Here, Tg is the brightness temperature in K, c is the speed of light in ms^\ 
By is the flux density per unit sohd angle, D., in Wm~^Hz~^ sr~\ V is the 
observing frequency in Hz and k is Boltzmann's constant in J/K. 

Q. can be computed as the area of a circle intersected by a cone with 
opening angle j3 (the FWHM diameter of the observing beam), yielding 
n = 4;rsin'^(j3/4), or, using sin(x) x for x << 1 and inttoducing units of 
mas, Q. « 1.846 x 10~^^ x j3^. Converting By and v to Jy/beam and GHz, 
respectively, yields 

TB = l.76xl0^^m^Ks-^r^^^. (5.31) 

The high brightness temperatures that we calculated are shown in Table 5.3 
and provide good evidence that the emission process of the radio components 
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Epoch 


Freq. 


■^max 


7b 




GHz 


mjybeam"' 


10* K 




Component A 




Nov 20, 1999 


5.0 


25 


140 




15.4 


6.21 


3.6 


Mar 6, 2000 


5.0 


24 


130 




15.4 


7.02 


4.0 


Nov 30, 2000 


5.0 


23 


130 




15.4 


17.5 


10 


Sep 22, 2002 


1.7 








2.3 


1.13 


29 




5.0 


22 


120 




Component B 




Nov 20, 1999 


5.0 


13.9 


76 




15.4 


46 


26 


Mar 6, 2000 


5.0 


14.6 


80 




15.4 


54 


31 


Nov 30, 2000 


5.0 


14.6 


80 




15.4 


62 


36 


Sep 22, 2002 


1.7 








2.3 








5.0 


14.8 


81 




Component E 




Nov 20, 1999 


5.0 


3.45 


19 




15.4 






Mar 6, 2000 


5.0 


3.22 


18 




15.4 






Nov 30, 2000 


5.0 


3.96 


22 




15.4 






Sep 22, 2002 


1.7 


2.94 


140 




2.3 


19.8 


510 




5.0 


5.57 


30 




Component F 




Nov 20 1 999 


5.0 
15.4 






Mar 6, 2000 


5.0 
15.4 






Nov 30, 2000 


5.0 








15.4 






Sep 22, 2002 


1.7 


9.30 


440 




2.3 


2.18 


56 




5.0 







Table 5.3: Brighmess temperatures of NGC 3079 components 



is synchrotron emission. The maximum spectral index of synchrotron radia- 
tion, however, is +2.5 (Rybicki & Lightman 1979, eq. 6.54), if the energies of 
the electrons in the jet are power-law distributed. Higher spectral indices can 
only be due to external media, and we therefore are confident that the compo- 
nents with a > +2.5 are obscured by a foreground absorber. We identify the 
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absorption mechanism as free-free absorption with 

5v = 5v,o e-< = Syfi v"^ e-'"""''" (5.32) 

where 5v and 5v,o are the observed and intrinsic flux densities of a component, 
respectively, Oq is the intrinsic spectral index, is the optical depth as given in 
Eq. 1. 12 and is the optical depth without the frequency term. The frequency 
dependence of causes the flux density to decrease exponentially towards 
lower frequencies, so that the spectral index can be arbitrarily high. Eq. 5.32 
has three unknowns: Sy^, Oq and T^. We have detected both A and E at three 
frequencies and can therefore solve for the three unknowns (we do not have 
three simultaneous measurements for A, and will therefore take the 5.0 GHz 
and 15.4 GHz data from November 2000 and the 2.3 GHz measurement from 
September 2002). The results are given in Table 5.4. The emission measure, 
k, is given by 

k = nil, (5.33) 

where L is the line-of-sight path length through the medium with electron 
density n^, assuming that T = 10^ K and N+ = N- = n^. For A and E, we 
find emission measures of 6.4 x lO^pccm^^ and 5.2 x lO^pccm^^. For an 
estimate of the line-of-sight column densities, n^L, we assume that the ab- 
sorber has the same depth as the FWHM of the synthesized beam. At the dis- 
tance of NGC 3079, 1 mas corresponds to 0.073 pc, and a circular beam with 
the same area as the elliptical beam used for the flux density measurements, 
4mas X 3 mas, has a FWHM of 3.46 mas, or 0.25 pc. The column densities then 
follow to be 1.23 x lO^^cm"^ for A and 1.11 x lO^^cm"^ for This is similar 
to the absorbers that have been found in, e.g., NGC 1052 (3.7 x 10^^ cm'^, Ka- 
meno et al. 2001), NGC 4261 (5 x lO^^cm'^, Jones et al. 2001) and Hydra A 
(3.7 X 10^^ cm^^, Taylor 1996) based on the same technique, and so indicates 
that the absorber in front of A and E has a similar width and similar physical 
conditions. 



Component 


Sv,o 


Oq 


noHz 


k n^L 




mjy 






pc cm^* cm^^ 


A 


180 


-0.8 


21 


6.4x10' 1.23x10^^ 


E 


10000 


-4.1 


17 


5.2x10' 1.11x10^2 



Table 5.4: Intrinsic properties of components A and E. Oo and are dimensionless quantities. 
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Component A is reasonably well modelled with an intrinsically optically 
thin (cUq = —0.8) synchrotron source behind a moderately thick (TfQjj^ = 21, 
'^Ighz ~ 0-72) free-free absorber The fit results from component E (tf^jj^ = 
17, t|qu2 = 0.58), however, raise an uncomfortable problem: Our model pre- 
dicts the intrinsic spectral index to be —4.1, much less than commonly ob- 
served in synchrotron sources. This result is robust with respect to changes 
of the flux densities at 1.7 GHz and 2.3 GHz by as much as factors of two, 
and so amplitude calibration errors and measurement errors from the image 
would have only Uttle effect. Such extremely steep spectra are not impossible 
to explain. One possibility is that the electrons in the source age and produce 
an exponential cutoff at high frequencies, but this process requires typically 
10^ yr in an undisturbed environment without further energy injection and so 
is unlikely to happen within parsecs of an AGN. Another possibiUty is that 
the electrons are simply not power-law distributed, but have a relativistic ther- 
mal distribution. Such models have been considered by, e.g., Jones & Hardee 
(1979) and Beckert & Duschl (1997), and they usually reveal an exponential 
cutoff at high frequencies. Both explanations are unsatisfactory because they 
require the electron distribution cutoff to he very close to the frequency at 
which T = 1, but we see no better solution. 

The limits on aj-j and aj j of component B might also be due to free- 
free absorption, but as B was detected only at 5.0 GHz and 15 GHz, we cannot 
solve Eq. 5.32. We can, however, draw a conclusion from its extent. B and 
E are resolved in the 5.0 GHz images from November 1999, March 2000 and 
November 2000. If the inverted spectrum of B was due to synchrotron self- 
absorption, one would expect it to be compact, so its extend indicates that it 
is unUkely to be synchrotron self-absorbed and the inverted spectrum is more 
Ukely to be due to free-free absorption. 

Steep spectra Component F has a spectral index of —1.64. Its turnover 
frequency Ues below the frequency range covered by our observations, so we 
cannot make a strong statement about physical parameters. But as all other 
components have spectral indices > 4.00 in parts of their spectra, we can as- 
sume that this is also the case for F, and consequently, that an east-west density 
or path length gradient in the free-free absorber has moved the peak frequency 
to frequencies < 1.7 GHz. 

Correlation of peak frequencies and riglit ascension The spectra 
shown in Fig. 5.3 have been arranged in order of decreasing right ascensions 



114 5 .2. Additional Results from NGC 3079 



of the components, and one sees a systematic shift in the turnover frequency 
towards higher frequencies from top to bottom. This arrangement suggests a 
connection between the peak frequencies and the component position in the 
source. If this effect is not accidental (which has a 12.5 % probability), it 
might be due to the sub-structure of a free-free absorber in front of the radio 
source. Free-free absorption has been suggested to be involved in A, E, and 
probably also B, so perhaps a single absorber covers all components, whose 
optical depth increases towards the west. As an example, it could be that the 
absorber has the shape of a wedge, the thin part of which is in front of F and 
the thicker part of which is in front of B. 
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Figure 5.3: Spectra of components in NGC 3079 with the spectral indices labelled. The turnover frequency 
of each component increases with decreasing right ascension, suggestive of a wedge-shaped absorber in 
front of the AGN. All spectral indices were calculated from the September 2002 data, except for aj5'4 for 
component A and B, which were calculated from the November 2000 data to avoid variability effects. 



5.2.2 Proper Motion 

So far, the location of the AGN in NGC 3079 has not been accurately 
determined. Based on VLBI maser observations. Trotter et al. (1998) have 
suggested that the core lies on the A — B axis, where that axis intersects the 
line of masers, and this identification has been supported by Sawada-Satoh 
et al. (2000) and, recently, by Kondratko (2003). 
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Figure 5.4: Relative positions among components A, B and E in NGC 3079 at 5.0 GHz for all epochs with 
arrows showing the direction of relative motion. 
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The relative positions of components A, B and E measured at 5.0 GHz are 
listed in Table 4.5 and are plotted in Fig. 5.4. In each plot, one component 
serves as a stationary reference point, and the arrows indicate in which direc- 
tion the other two components are moving. The errors are s/l x 0.15 mas = 
0.21 mas, because the error of the reference point was added to the other errors. 

A stationary If A is chosen as a stationary reference point, then B is seen 
to move radially away from A, and E is moving almost perpendicular to the 
A — E axis. The motion of E would not be expected in the standard AGN 
model, where radio components move radially away from the AGN. If E was 
ejected by the AGN located at or near A, then its path must have a large bend 
that first makes it travel away from A and then directs it in the same direction 
in which B is travelUng. We do not consider this scenario further. 

B stationary If B is the reference, then A is seen to move radially away 
from B, but E is moving towards the A — B axis. This situation, like the pre- 
vious, does not agree very well with the standard model, in which one would 
not expect components to move in a direction away from the AGN. Never- 
theless, this scenario is of particular importance. Sawada-Satoh et al. (2000), 
comparing their maser and continuum observations to those of Trotter et al. 
(1998), found B to be stationary with respect to the brightest maser features. 
In this case, plotting the separations relative to B yields the same speeds and 
directions as plotting them relative to the masers and hence to the true (radio- 
invisible) AGN centre of mass location, and E would indeed be moving in the 
direction of the nucleus. Furthermore, Kondratko (2003) recently found good 
evidence for a rotating disc centred on the A — S axis, confirming those results. 

E stationary If E is the reference, then A and B move in almost perpen- 
dicular directions. This scenario has one attractive advantage, namely that the 
back-extrapolation of the trajectories of A and B intersect in the past. This al- 
lows them to be ejected from the same location, south-west of B. E could be 
an object that is kinematically disconnected from the nucleus and is not part of 
the AGN ejecta. An obvious guess is to identify it as a supernova, supported by 
the fact that E has not been observed by Trotter et al. (1998) in 1992, but was 
detected by Sawada-Satoh et al. (2000) in 1996, although they misidentified it 
with B. Thus, E must have appeared between 1992 and 1996. At 5.0 GHz, it 
has an average flux density of 4.89 mJy, corresponding to 1.3 x lO^^WHz^^ 
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at a distance of 15 Mpc. This is comparable to the brightest supernova in, e.g., 
M82, 41.9 + 58, having 4.6 x IQi^WHz-i at 5 GHz (Kronberg et al. 1985), 
although most of the supernovae reported in that publication have ten times 
lower luminosity. However, E lacks the typical decrease in flux density that is 
seen as the supernova shells expand. Although Sawada-Satoh et al. (2000) find 
aflux density of (5.7 it 1.7)mJy at 1.4GHz, compared to (3.06±0.15)mJy in 
our 1 .7 GHz observations, the integrated flux densities were (4.15±0.21)mJy, 
(3.08±0.15)mJy, (5.20±0.26)mJy and (7.1 1±0. 36) mJy in our 5.0 GHz ob- 
servations. This sequence of measurements does not suggest a gradual decline 
of flux density with time like has been observed in M82's 41.9 + 58 during 
1971-1981, fading with a rate of between 9%yr-i and 24%yr-i (Kronberg 
et al. 1985). 

Core invisible It is also possible that none of the radio emitters is coinci- 
dent with the dynamical centre. This situation is common in radio galaxies 
and quasars where one usually sees a T = 1 surface in the jet, and not a com- 
ponent at the location of the black hole, but as a Seyfert 2 galaxy, NGC 3079 
is certainly an entirely different object. If the dynamical centre is invisible in 
the radio, the components may be shocks in an AGN outflow. 

IHelical paths Another interesting interpretation is that all components 
are moving on paths that follow a heUx with a fairly wide opening angle. If 
the axis of the helix is inclined at a small angle to the line of sight, inwards 
motion of components would be easily explained. 

None of the five scenarios sketched above yields compelling evidence for 
the location of the core. However, we think that the model in which B is sta- 
tionary (it might be equivalent to the model in which the core is invisible) and 
the helical paths model need to be considered in the future. Maser observa- 
tions are usually regarded as one of the "best possible" pieces of evidence for 
the location of black holes in AGN, but the kinematic centres of circumnu- 
clear discs do not necessarily need to coincide with these cores. This has been 
found in, e.g., NGC 4261 (Ferrarese et al. 1996, see §3.3), although in that 
case, the multiple centres are interpreted as being due to past merging events, 
and spiral galaxies like NGC 3079 are generally not considered as descendants 
of galaxy mergers. We plan to monitor NGC 3079 during the next two years 
using phase-referencing with good frequency coverage between 1 .7 GHz and 
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22 GHz, allowing us to follow absolute proper motions to resolve the ambigu- 
ous findings and solve for the core position. 

5.2.3 Remarks on Individual Components 

Component A is the most compact component in NGC 3079. The ratio of 
peak flux density to integrated flux density, x, has an average of 0.85 when 
measured from tapered images and of 0.69 when measured from the full 
resolution images from November 1999, March 2000 and November 2000 
(Table 4.3). Hence, it is only slightly resolved. 

Component B is clearly resolved at 5.0 GHz in the first three epochs 
(x = 0.39) and compact at 15 GHz (x = 0.83). This is probably not a res- 
olution effect because the 5.0 GHz and 15 GHz images have similar resolu- 
tion (the 15 GHz resolution even being sUghtly smaller), but a combination 
of steep spectral index and lower sensitivity can in principle account for the 
effect. The extension of B to the south in the 5.0 GHz images has four to 
five contours, the lowest contour being drawn at three times the image rms. 
The extension therefore has a surface brightness of 8 to 16 times 3 o, which 
is of the order of 1.6mJybeam"^ to 3.2mJybeam^^ Assuming a spectral 
index of —0.7 as is common for optically thin synchrotron radiation then pre- 
dicts a surface brightness at 15.4 GHz of 0.7mJybeam^^ to 1.5mJybeam^^ 
This is visible as an extension of the lowest contour in the 15 GHz images 
which are drawn at around ImJybeam^^ especially in the November 2000 
epoch, where this contour is at 0.82mJybeam^^. Furthermore, Kondratko 
(2003) finds the same extension of B at 22 GHz to have a surface brightness of 
(0.42 to 0.84)mJybeam~\ consistent with a similar extrapolation of the sur- 
face brightness found by us at 5.0 GHz. Thus, the northern and compact part 
of B has an inverted spectrum between 5.0 GHz and 15.4 GHz, whereas the 
extension to the south fades towards higher frequencies. If these two differ- 
ent spectral shapes on such small linear scales are interpreted as being due to 
free-free absorption, the situation does not agree well with the idea of a sin- 
gle, wedge-shaped absorber in front of the radio components as suggested in 
§5.2. 1. But there are too few observational constraints to conclusively tell what 
causes the spectral gradient in B. 

Component E is also resolved, with an average peak-to-integrated ratio 
of 0.30 in the full-resolution images. Furthermore, the position angle of its 
largest extent appears to change between the first three epochs. We have fitted 
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elliptical Gaussians to the (m,v) data of all 5.0 GHz epochs, using data with a 
(m,v) radius of < 100 MA during the first three epochs and using data with a 
(m, v) radius of < SOMA of the fourth epoch. The limits were chosen to reduce 
the ellipticity of the (w, v) coverages while keeping as many data as possible. 
The results of the fits are listed in Table 5.5. The P. A. changes between the 
first three epochs are 28.7° and 55.8°, whereas the P.A. change between the 
November 2000 and September 2002 epochs is difficult to determine due to 
the nn ambiguity of the P.A. However, the rates at which the P.A. changed is 
97.2° yr-i between November 1999 and March 2000 and 76.3° yr"! between 
March 2000 and November 2000. Extrapolating over the 1.811 yr period be- 
tween the third and fourth epoch yields P.A.s of 231.8° (= 51.8° when the PA.s 
are reflected onto the range 0° to 180°) for a rate of 97.2° yr^^ and 269.8° 
(= 89.8°) for a rate of 76.3° yr \ or an average P.A. of 70.8°. This is very 
close to the observed P.A. of 65.0°. 

However, we deem this result not very reliable. First, the P.A.s changed 
considerably when the maximum (m,v) radius in imaging was changed and 
second, the P.A.s predicted for the September 2002 observation, 51.8° and 
89.8°, span 38.0°, and hence there is a probability of 100 x (38.0°/180°) = 
21 % of the measured PA. falling into that range if it was random. The P.A. 
changes of E are not the result of a changing (m,v) coverage, as the restoring 
beam after tapering the images is independent of the P. A. of E. 

Combining this result with the observed spectra and proper motion, we 
find the least exotic explanation for the radio emission of component £ is an 
evolving shock in a large jet-like outflow. 
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deg 


Nov 1999 


2.69 


2.20 


-18.3 


2.50 


1.42 


9.2 


Mar 2000 


1.85 


1.68 


-29.6 


2.09 


1.17 


37.9 


Nov 2000 


2.35 


1.86 


-24.8 


2.32 


1.57 


93.7 


Sep 2002 


4.16 


3.81 


-50.6 


2.28 


2.07 


65.0 



Table 5.5: Position angle changes of component E. 



5.3 Statistics 



We have carried out a literature review to compare the findings in NGC 3079 
to other Seyfert galaxies and AGN to address the following questions. 
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• The spectral indices of the radio components in NGC 3079 were found to 
be unusually inverted. Is this peculiar to NGC 3079 or do other Seyferts 
have inverted spectra as well? The result might indicate whether free- 
free absorption is common in Seyferts or not. 

• NGC 3079, like some other Seyferts, shows a kpc-scale outflow, and 
if pc-scale components tend to be ejected in the same direction, this 
could help in the identification of components in NGC 3079. Is there a 
general correlation between the host galaxy rotation axes in Seyferts and 
the pc-scale emission? This might indicate that Seyfert radio ejecta are 
undergoing significant interactions with interstellar gas, or precession 
in the accretion discs, whereas powerful objects frequently are aligned 
from the innermost regions to the outer radio lobes. 

• How do the VLBI-measured spectral indices in Seyferts generally com- 
pare to those measured with smaller arrays? 

5.3.1 Compilation of the Sample 

The range of radio properties displayed by Seyfert galaxies has been explored 
extensively based on WSRT, VLA or ATCA observations at 1 .4 GHz, 5 GHz or 
8.4 GHz at arcsec resolution by de Bruyn & Wilson (1978), Ulvestad & Wilson 
(1984, 1989), Unger et al. (1987), Giuricin et al. (1990), Rush et al. (1996), 
Nagar & Wilson (1999), Morganti et al. (1999), Ulvestad & Ho (2001a), and 
Vir Lai (2001). They studied the range of radio luminosities, radio source sizes, 
morphologies, alignment between jet and host galaxy, and specttal indices, 
and looked for dependence of properties on Seyfert type. With the exception 
of Ulvestad & Ho (2001a), who found 50% of low-luminosity Seyferts to 
have a flat spectrum core, they all found that most Seyfert galaxies have steep, 
optically-thin synchrotron radio spectra, with 5 % to 35 % (depending on the 
sample) showing a flat-spectrum core, much less than the rate of occurrence 
of flat-spectrum cores in powerful radio galaxies. No difference was found 
between the spectral index distributions of the two Seyfert types. 

However, the 0.3" to several arcsec beams used in those studies do not al- 
low details of nuclear structure to be separated or the distribution of jet bending 
angles to be constructed. Further, they would probably blend together multi- 
ple radio components within the nucleus and could hide the presence of more 
absorbed cores. 
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Many further VLBI observations of Seyfert galaxies are available in the lit- 
erature but no statistical summary from those observations has been presented. 
Thus, we made a Uterature search for VLBI observations of Seyferts to com- 
pare to our observations of NGC 3079. The resulting list in Table 5.6 contains, 
we believe, all VLBI imaging observations of AGN classified as Seyfert 1, 2 
or 3 by Veron-Cetty & Veron (2001) (irrespective of any sub-classification), 
published as of August 2003 with a linear resolution of < 10 pc, although it 
is not impossible that we have inadvertently omitted some works. This sam- 
ple consists mostly of the nearest and best-known Seyferts. They were often 
selected for observing because they show bright radio nuclei at arcsecond reso- 
lution, and so the sample is not complete in any sense. Large, uniform samples 
of Seyfert galaxies are being observed with the VLBA by Murray et al. (1999) 
and Schmitt et al. (2002), but the results are not yet available. Proper mo- 
tions in Seyferts are being presented by Roy et al. (in prep.) and we confine 
ourselves here to spectral indices and jet bending. 

Spectral indices have been taken from the literature if they were measured 
simultaneously between L-band and C-band, i.e., between 1.4 GHz to 1.7 GHz 
and 4.9 GHz to 5.0 GHz. To compare pc and kpc-scale structure, we have 
searched for interferometric radio observations with a linear resolution of > 
lOOpc, made with, e.g., the VLA or MERLIN. The angles have been measured 
from the images by us, from the identified core or the brightest component to 
the next bright component. The complete list is shown in Table 5.6, those 
objects used for determining spectral indices are listed in Table 5.7, and the 
angles are shown in Table 5.10. We have plotted the data in two histograms 
shown in Fig. 5.5. 

5.3.2 Spectral Properties of Seyferts 

We selected from Table 5.6 a subsample consisting of the Seyfert galaxies for 
which dual-frequency VLBI observations have been made to study the compo- 
nent spectral properties at high resolution. This subsample is listed in Table 5.7 
and includes the spectra of A, B and E in NGC 3079 . 

Spectral Index Distributions 

The distribution of ^ from Table 5.7 is shown in Fig. 5.5 (top). The spectral 
indices lie between —1.8 and +4.6 with a median of —0.40 (using the Kaplan- 
Meier estimator constructed by ASURV rev 1 .2, Lavalley et al. 1992 to account 
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for lower limits and discarding upper limits). The distribution shows a tail to- 
wards positive spectral index due to four components in NGC 1068, NGC 2639 
and NGC 3079 that have strong low-frequency absorption. 

Separating the VLBl measurements by Seyfert type, we found the median 
ttg '* for the Sy Is was —0.75 and for the Sy 2s was —0.23, which were not 
significantly different under the survival analysis two-sample tests (P = 6 % to 
12 % that the spectral-index distributions came from the same parent distribu- 
tion) or under a Kohnogorov-Smimov test (P = 6 % to come from the same 
parent distribution), treating only lower hmits as detections. 

Despite the statistical insignificance of a difference between Seyfert types, 
we notice that the most inverted spectra all occur in objects that have evi- 
dence for an absorber. NGC 3079 has good evidence for free-free absorption, 
NGC 1068 and Mrk348 are type 2 Seyferts with type 1 spectra in polarized 
light and NGC 2639 shows strong X-ray absorption (Wilson et al. 1998). The 
inner edge of such an absorber is expected to be ionized by UV emission from 
the AGN which would produce free-free absorption and hence produce in- 
verted spectra preferentially in the type 2s. All four objects with inverted radio 
spectra also host H2O masers, which is also consistent with an orientation ef- 
fect if, for example, the radio core is viewed through a molecular disc or torus, 
which produces H2O maser emission. 

With spectral indices of a^-'^{A) > 4.56, al'^{B) > 4.41 and a^-'^{E) = 
0.67, NGC 3079 has contributed the two most inverted spectra seen in any 
Seyfert on parsec scales. However, many more radio components in Seyferts 
have flat or inverted spectra, and NGC 3079 seems to be an extreme object, 
possibly because it is an edge-on spiral galaxy and so there is more absorbing 
gas in the line of sight than in face-on Seyferts. 

Spectral Indices Measured with Smaller Arrays 

The tail of extreme positive spectral indices that we see in our sample is not 
seen at arcsec resolution in the works mentioned in §5.3.1, probably due to 
component blending. We made a statistical comparison of our VLBI spectral 
indices to arcsec resolution spectral indices using the Kolmogorov-Smimov 
test and chi-squared homogeneity test by treating the lower and upper limits 
as detections, and by applying survival analysis two-sample tests. The results 
are summarized in Table 5.8. Our spectral index distribution was significantly 
different (i.e. level of significance was less than the critical value of 5 %) 
from all the comparison samples according to all tests performed, with the 
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exception of the Palomar Seyferts by Ho & Ulvestad (2001), which was found 
to be significantly different in only the chi-squared test. This excess of flat- 
spectrum sources in the Palomar Seyfert sample was already noticed by Ho & 
Ulvestad (2001). 

Fraction of Flat-Spectrum Components 

Of the 16 objects in the subsample, a majority of 13 (81 %) has at least one 
flat or inverted spectrum VLBI component with ttj-^ > —0.3 (denoted by an 
asterisk in Table 5.7). This fraction is 60 % larger than that measured at lower 
resolution by Ho & Ulvestad (2001) for the Palomar Seyferts and four or 
more times larger than the fraction measured at arcsec resolution for the other 
samples (Table 5.9). The difference is highly significant for most samples, 
having a probability of < 0.1 % of occurring at random. The larger fraction of 
flat- and inverted-spectrum components found with VLBI is not unexpected 
since increasing resolution reduces the component blending. 



Summary 

Flat- or absorbed- spectrum cores are present in 13 of 16 Seyfert galaxies. The 
most strongly absorbed components tend to be in objects that also show other 
evidence for absorption, such as water maser emission or large X-ray absorp- 
tion columns. In those cases, according to the Seyfert unification schemes, 
the flat-spectrum components would be viewed through an edge-on obscur- 
ing disc or torus, the ionized inner edge of which could be causing free-free 
absorption. If so, the situation in those objects is different from the Blandford- 
Konigl jet model. However the case for free-free absorption is unambiguous in 
only NGC 3079 for which spectral indices are > -1-2.5. The other 12 Seyferts 
have —0.3 < a < +2.5 and so spectral turnovers could be produced by either 
synchrotton self-absorption or free-free absorption. 

5.3.3 Jet Misalignment 

Kiloparsec-scale linear radio sttuctures in Seyfert galaxies are thought to be 
energized by jets or outflows originating at the AGN core, and so one might 
expect pc and kpc-scale radio structures to be aligned. Misalignments could be 
caused by changes in the ejection axis or by pressure gradients in the ISM. To 
look for such effects, we compiled kpc-scale radio observations of the VLBI 



124 



5.3. Statistics 



Seyfert sample, where available. The resulting sample of 21 objects in Ta- 
ble 5.10 is, we believe, a complete list of Seyferts with both pc-scale and kpc- 
scale observations published, although seven of these objects either have no 
clear extended structure on one of the scales or the structure was deemed to be 
affected by processes on scales > Ikpc, e.g., buoyancy forces that are unre- 
lated to changes in the jet ejection axis. In one object, NGC 1167, the double- 
sided kpc-scale emission is bent, so we quote both. This sample has a large 
overlap with the sample in Table 5.10 and, likewise, is not complete. Biases 
such as selection for bright radio nuclei are probably present, but whether they 
affect the distribution of jet bending angle is unclear. Doppler boosting would 
favour jets that point towards us, but the jet speeds so far measured in Seyferts 
are low and boosting effects should be minor. 

The distribution of misalignment angles from Table 5.10 is shown in 
Fig. 5.5 (bottom). Misalignment in Seyfert jets between pc and kpc-scales 
is common, with 5 out of 15 objects bending by 45° or more. The null hy- 
pothesis that the jet bending-angle distribution was drawn from a parent pop- 
ulation of aligned jets was rejected by a chi-squared test at the 2 % level of 
significance. The null hypothesis that the distribution was drawn from a uni- 
formly distributed parent population could not be rejected by a chi-squared 
test, which returned a 25 % level of significance. Thus, the misalignment be- 
tween the pc-scale and kpc-scale radio structure in NGC 3079 is common in 
Seyfert galaxies. 

The frequent misalignments seen in Seyfert galaxy jets could be evidence 
either for changes in the jet ejection angle and hence of the accretion disc, per- 
haps by a mechanism such as radiation-driven warping (Pringle 1997), or for 
bending due to pressure gradients in the ISM (collision or buoyancy) through 
which the jets propagate. 

We compared the extended pc-scale radio structures to the host galaxy rota- 
tion axes, assuming that the host galaxy is circular and the rotation axis vector 
is projected onto the galaxy's minor axis, which can be measured from optical 
images. For all spiral galaxies with pc-scale radio structure, we used the major 
axes from de Vaucouleurs et al. (1991), subtracted 90° and reflected the angles 
onto the range 0° to 90°. Four objects were not listed in de Vaucouleurs et al. 
(1991), so we used data from Schmitt & Kinney (2000) in two cases and we 
measured the major axes from Hubble Space Telescope archival data in the 
other two cases. We performed a Spearman rank correlation test and found 
a correlation coefficient of 0.07, corresponding to a 82 % significance and so 
there is no correlation between the directions of the pc-scale structures and 
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the host galaxy rotation axes. However, the sample was small (13 objects) and 
therefore the test is not strong. Thus, the lack of aUgnment in Seyferts between 
pc-scale and kpc-scale structure does not allow one to identify a pc-scale ejec- 
tion axis in NGC 3079. Precession of the ejection axis is common and may be 
reflected in the confusing structure in the core of NGC 3079. 

We compared the bending angles in Seyferts to those in radio-loud objects 
using a comparison sample of core-dominated radio-loud objects compiled 
from Table 6 in Pearson & Readhead (1988), and Table 1 in Conway & 
Murphy (1993), which have both VLBI and VLA or MERLIN observations at 
5 GHz. Those objects show a surprising bimodal distribution of jet bending 
angle, in which the pc-scale jets tend to align with or to be perpendicular to 
the kpc-scale structure. We compared our jet bending angle distribution to 
the distribution in radio-loud objects using a Kolmogorov-Smimov test and 
found no significant difference, the probabihty of being drawn from the same 
parent distribution being 14%. However, the Seyfert sample is small and 
the distribution for the Seyferts was shown in the previous paragraph to be 
consistent also with a uniform distribution. 



Summary 

The position angle differences between pc-scale and kpc-scale radio emission 
were found to be uniformly distributed between 0° and 90°. Such bends could 
be due to changes in the jet ejection axis or due to pressure gradients in the 
ISM. No correlation was found between the axis of pc-scale radio structure 
and the rotation axis of the host galaxy. 

We conclude that NGC 3079 shares most properties with other Seyferts. 
Its spectra are extreme, but three other objects also showed inverted spectra 

and the transition between these three and NGC 3079 does not appear to mark 
a fundamental difference, and misalignment between pc-scale and kpc-scale 
radio emission also is a common finding. 



Table 5.6: Seyferts observed with VLBI. Column 1 lists the names we used, column 2 the names used by Veron-Cetty & Veron (2001), column 3 the 
Seyfert types as described by Veron-Cetty & V6ron (2001) , column 4 the linear resolution in pc, column 5 the instrument and frequency used, and 
column 6 the reference to the publication. 



Name 


V6ron name 


Class 


Lin. Res. 


Instrument 


Ref 


3C 287.1 




SI 


5.43 


VLBA 5 GHz 


Fomalont et al. (2000) 


3C 390.3 




S1.5 


1.09 


global VLBI 5 GHz 


Alef et al. (1996) 


III Zw 2 




S1.2 


0.260 


VLBA 15/43 GHz 


Brunthaler et al. (2000) 


Ark 564 




- 


0.598 


global VLBI 5 GHz 


Vir Lai (2001) 


IC 5063 


PKS 2048-57 


Slh 


3.30 


LBA 2.3 GHz + LBA HI 


Oosterloo et al. (2000) 


MCG 8-11-11 




S1.5 


0.635 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 1 




S2 


6.18 


LVN 1.7 GHz 


Kukula et al. (1999) 


Mrk 1 




S2 


0.550 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 3 




Slh 


5.24 


EVN1.7GHZ 


Kukula etal. (1999) 


Mrk 78 




S2 


1.75 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 231 




Sl.O 


0.343 


VLBA 15 GHz 


Ulvestad et al. (1999b) 


Mrk 231 




Sl.O 


4.09 


VLBA 1.4/2.3/4.8/8.4/15/22 GHz 


Ulvestad et al. (1999a) 


Mrk 273 




S2 


7.32 


VLBA 1.4 GHz 


Carilli & Taylor (2000) 


Mrk 348 




Slh 


0.131 


VLBA 15 GHz 


Ulvestad et al. (1999b) 


Mrk 348 




Slh 


6.41 


EVN 1.4 GHz 


Neff &deBruyn(I983) 


Mrk 348 




Slh 


1.49 


VLBA 1.7/5.0/8/15/22 GHz 


Peck et al. (2003) 


Mrk 477 




Slh 


2.39 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 530 


NGC 7603 


S1.5 


1.45 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 766 




S1.5 


0.413 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 926 




S1.5 


0.917 


VLBA 8.4 GHz 


Mundell et al. (2000) 


Mrk 1210 




Slh 


0.392 


EVN 1.6/5 GHz VLBA 1.6/5 GHz 


Middelberg et al. (2004) 


Mrk 1218 




S1.8 


0.582 


global VLBI 5 GHz 


Vir Lai (2001) 


NGC 1052 




S3h 


0.015 


VLBA 5/8.4/22/43 GHz 


Kadler et al. (2002) 


NGC 1052 




S3h 


0.048 


VLBI 1.4-43.2 


Vermeulen et al. (2003) 


NGC 1052 




S3h 


0.038 


VLBA 2.3/8.4/15.4 GHZ 


Kameno et al. (2001) 


NGC 1068 




Slh 


0.221 


VLBA 1.7,5,15 GHz 


Roy et al. (1998) 


NGC 1068 




Slh 


2.94 


EVN 1.4 GHz 


Ulvestad et al. (1987) 


NGC 1167 




S3 


0.895 


global VLBI 5 GHz 


Giovannini et al. (2001) 
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Name V&on name Class Lin. Res. Instrument Ref. 



NGC 1167 


S3 


0.639 


global VLBI 1.6 GHz 


Giovannini et al. (1990) 


NGC 1275 


S1.5 


0.416 


VLBA 2.3/5/8.4/15.4/22/43 GHz 


Walker et al. (2000) 


NGC 2110 


Sli 


0.151 


VLBA 8.4 GHz 


Mundell et al. (2000) 


NGC 2110 


Sli 


0.966 


EVN 1.6/5 GHz 


Middelberg et al. (2004) 


NGC 2273 




0.342 


global VLBI 5 GHz 


VirLal (2001) 


NGC 2639 


S3 


0.233 


global VLBI 5 GHz 


Vir Lai (2001) 


NGC 2639 


S3 


0.162 


VLBA 1.6/5/15 GHz 


Wilson et al. (1998) 


NGC 3079 


S2 


0.044 


VLBA 5/8/22 GHz 


Trotter et al. (1998) 


NGC 3079 


S2 


0.054 


global VLBI, 5.0 GHz 


Irwin & Seaquist (1988) 


NGC 3079 


S2 


0.042 


global VLBI 1.4/8.4/15/22 


Sawada-Satoh et al. (2000) 


NGC 3147 


S2 


0.219 


VLBA 1.6/2.3/5/8.4 GHz 


Ulvestad & Ho (2001b) 


NGC 3227 


S1.5 


3.74 


MERLIN 1.6/5 GHz 


Mundell et al. (1995) 


NGC 4151 


S1.5 


0.116 


VLBA 1.6,5 GHz 


Ulvestad et al. (1998) 


NGC 4151 


S1.5 


1.29 


EVN 1.7 GHz 


Harrison et al. (1986) 


NGC 4151 


S1.5 


4.82 


MERLIN 5 GHz 


Pedlar et al. (1993) 


NGC 4168 


S1.9 


0.236 


VLBA 5 GHz 


Nagar et al. (2002) 


NGC 4203 


S3b 


0.084 


VLBA 1.6/2.3/5/8.4 GHz 


Ulvestad & Ho (2001b) 


NGC 4258 


S2 


0.029 


VLBA 22 GHz 


Hermstein et al. (1997) 


NGC 4258 


S2 


0.232 


VLBA 1.4/1.6 GHz 


Cecil et al. (2000) 


NGC 4395 


S1.8 


0.221 


VLBA 1.4 GHz 


Wrobel et al. (2001) 


NGC 4565 


S1.9 


0.207 


VLBA 5 GHz 


Falcke et al. (2000) 


NGC 4579 


S3b 


0.118 


VLBA 1.6/2.3/5/8.4 GHz 


Ulvestad & Ho (2001b) 


NGC 4579 


S3b 


0.246 


VLBA 5 GHz 


Falcke et al. (2000) 


NGC 5252 


S2 


0.671 


VLBA 8.4 GHz 


Mundell et al. (2000) 


NGC 5506 


Sli 


0.815 


EVN 1.6/5 GHz 


Middelberg et al. (2004) 


NGC 5548 


S1.5 


0.366 


VLBA 8.4 GHz 


Wrobel (2000) 


NGC 5793 


S2 


0.135 


VLBA 1.7,8.4,15,22 GHz 


Hagiwara et al. (2001) 


NGC 5793 


S2 


0.316 


VLBA+Y27 1 .6/5 GHz 


Hagiwara et al. (2000) 


NGC 5929 


S3 


0.290 


global VLBI 5 GHz 


VirLal (2001) 


NGC 5929 


S3 


6.44 


MERLIN 0.4/L6/5 


Su etal. (1996) 


NGC 7212 


Slh 


1.37 


global VLBI 5 GHz 


Vir Lai (2001) 


NGC 7469 


S1.5 


0.844 


global VLBI 5 GHz 


VirLal (2001) 


NGC 7674 


Slh 


2.80 


VLBA+Y27+Arecibo 1.4 GHz 


Momjian et al. (2003) 
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Name 


V&on name 


Class 


Lin. Res. 


Instrument 


Ref. 


NGC 7674 




Slh 


2.58 


EVN 1.6/5 GHz 


Middelberg et al. (2004) 


NGC 7674 




Slh 


2.18 


global VLBI 5 GHz 


Vir Lai (2001) 


NGC 7682 




S2 


0.996 


global VLBI 5 GHz 


Vir Lai (2001) 


T0109-383 


NGC 424 


Slh 


0.226 


VLBA 8.4 GHz 


Mundell et al. (2000) 
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Source 


Flat/Inv. 


Comp. 


at' 


Refs. 


Mrk231 


* 


N 


-0.86 


Ulvestad et al. (1999a) 






C 


+0.88 


Ulvestad et al. (1999a) 






s 


-1.54 


Ulvestad et al. (1999a) 


Mrk 348 


* 




+0.93 


Barvaini.s & Lonsdale (1998) 


Mrk 463E 




L 


-1.06 


Norris et al., unpublished 






R 


-0.62 


Norris et al., unpublished 






1 


-0.67 


Norris et al., unpublished 






2 


-0.77 


Norris et al., unpublished 






3 


-1.00 


Norris et al., unpubUshed 


Mrk 1210 




NW 


-1.26 


Middelberg et al. (2004) 






SE 


-0.78 


Middelberg et al. (2004) 


NGC 1068 


* 


NE 


-1.5 


Roy et al. (1998) 






C 


>2.2 


Roy et al. (1998) 






SI 


>0.6 


Roy et al. (1998) 






S2 


<-1.9 


Roy et al. (1998) 


NGC2110 


* 


c 


+0.12 


Middelberg et al. (2004) 


NGC 2639 


* 




+ 1.78 


Wilson et al. (1998) 


NGC 3079 


* 


E 


+0.67 


This thesis 






A 


> +4.56 


This thesis 






B 


>+4.41 


This thesis 


NGC 3147 


* 




+0.20 


Ulvestad & Ho (2001b) 


NGC 4151 


* 


E 


-0.5 


Ulvestad et al. (1998) 






D 


< -0.2 


Ulvestad et al. (1998) 






F 


< -1.3 


Ulvestad et al. (1998) 


NGC 4203 






+0.38 


Ulvestad & Ho (2001b) 


NGC 4579 


* 




+0.20 


Ulvestad & Ho (2001b) 


NGC 5506 


* 


BO 


+0.06 


Middelberg et al. (2004) 






Bl 


-0.73 


Middelberg et al. (2004) 






B2 


-0.30 


Middelberg et al. (2004) 


NGC 5793 


* 


CIC 


-0.70 


Hagiwara et al. (2000) 






CINE 


-1.00 


Hagiwara et al. (2000) 






C2C 


-0.46 


Hagiwara et al. (2000) 






C2W 


-0.72 


Hagiwara et al. (2000) 






C2E 


+0.13 


Hagiwara et al. (2000) 


NGC 7469 


* 


2 


-0.34 


Norris et al., unpubUshed 






3 


+0.09 


Norris et al., unpubUshed 






5 


-0.12 


Norris et al., impubUshed 


NGC 7674 






-1.78 


Middelberg et al. (2004) 



Table 5.7: Seyferts from Table 5.6 that have dual-frequency (1 .4 GHz and 5 GHz) VLBI observations with 
matching beam sizes, from which the spectral index is derived. Asterisks denote objects with at least one 
component with a flat or inverted spectrum. 
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Level of Significance 



Sample 








KS 


Logrank 


Gehan 


Peto-Peto 


Ungeretal. (1987) 


X-ray Seyferts 


0.3' 


fo 


0.7% 


0.3% 


8.3% 


8.3% 


Rush et al. (1996) 


CfA Seyferts 


0.1' 


7c 


0.1% 


0.0% 


0.5% 


0.5% 


Rush et al. (1996) 


12 nm Seyferts 


0.1' 


7c 


0.6% 


0.0% 


2.0% 


2.0% 


Morgantietal. (1999) 


Distance Limited 


0.1' 


7c 


0.3% 


0.0% 


1.2% 


1.2% 


Ho & Ulvestad (2001) 


Palomar Seyferts 


1.0' 


7c 


27% 


44% 


75% 


77% 


Ulvestad & Wilson (1984) 


Distance Limited 


0.1' 


7c 


0.3% 


0.0% 


0.7% 


0.7% 



Table 5.8: Comparison of the spectral index distribution of the VLBI Seyfert sample to various other 
Seyfert samples. The level of significance gives the probability that the spectral index distribution of the 
comparison sample and the VLBI seyfert sample in this thesis were drawn from the same parent population. 



Sample 




Na>-03 


A'total 


% 


P 


Middelberg et al. (2004) 


own obs + lit search 


17 


38 


45 




Unger et al. (1987) 


X-ray Seyferts 


1 


20 


5 


0.3% 


Rush et al. (1996) 


CfA Seyferts 


3 


46 


6.5 


<0.1% 


Rushetal. (1996) 


12 Seyferts 


4 


40 


10 


< 0.1 % 


Morgantiet al. (1999) 


Distance limited 


13 


60 


22 


1.8% 


Ho & Ulvestad (2001) 


Palomar Seyferts 


21 


45 


47 


50% 


Ulvestad & Wilson (1984) 


Distance limited 


3 


40 


7.5 


< 0.1 % 



Table 5.9: Fraction of flat- and inverted-spectrum components in various Seyfert samples. The column 
'Na>-0.3' gives the number of components with spectral indices a > —0.3 (i.e., components with flat or 
inverted spectra). The column Wtotai' gives the total number of components in each sample. The column 
% gives (A'a>-o.3/Motai) x 100. The column P gives the level of significance from a comparison of the 
fraction of flat- and inverted-spectrum components in the VLBI sample of this thesis (limits treated as 
detections), to that of each comparison sample at arcsec resolution, using the di£ference-of-two-proportions 
test with Yates correction for continuity (e.g. Glanz 1992). 



Table 5.10: Subset of sources from Table 5.6 for which both pc and kpc-scale structures have been observed. The misahgnment angle, A0, is calculated 
as 0(pc) — 0(kpc), and has then been reflected onto the range of 0° to 90°. For the seven objects at the end of the table, either the P.A. information of 
the kpc scale structure was not completely reUable or was unavailable. The position angles of the galaxy minor axes in the last column were drawn from 
de Vaucouleurs et al. (1991), except for those marked a taken from Schmitt & Kinney (2000) and those marked h measured by us from Hubble Space 
Telescope archival images (see §5.3.3). 



00 



Name 


Res. 


0(pc) 


(kpc) 


A© 


Array 


Refs. 


©(opt.) 




(pc) 


(deg) 


(deg) 


(deg) 






deg 



3C 390.3 


1.09 


323 







global VLBI 5 GHz 


AM etal. (1996) 






1210 




323 




VLA 5 GHz 


AM etal. (1996) 




Mrk78 


1.75 


40 




55 


global VLBI 5 GHz 


Vir Lai (2001) 






504 




275 




VLA 5 GHz 


Vir Lai (2001) 




Mrk231 


0.343 


93 




77 


VLBA 15 GHz 


Ulvestad et al. (1999b) 


100 




3270 




170 




VLA 1.5/8.4/15 GHz 


Ulvestad et al. (1999a) 




Mrk 1218 


0.582 


330 




13 


global VLBI 5 GHz 


Vir Lai (2001) 


130* 




294 




317 




VLA 5 GHz 


Vir Lai (2001) 




NGC 1052 


0.015 


65 




30 


VLBA 5/8.4/22/43 GHz 


Kadler et al. (2002) 






162 




275 




MERLIN 1.4 GHz 


Kadler et al. (2003) 




NGC 1068 


0.221 


11 




29 


VLBA 1.7,5,15 GHz 


Roy et al. (1998) 


160 




217 




40 




VLA 5 GHz 


Ho & Ulvestad (2001) 




NGC 1167 


0.895 


140 


313 


7 


global VLBI 5 GHz 


Giovannim et al. (2001) 






131 




66 


74 


VLA 5 GHz 


Sanghera et al. (1995) 




NGC 1275 


0.416 


353 




23 


VLBA 2.3 to 43 GHz 


Walker et al. (2000) 






382 




150 




VLA 1.4 GHz 


Pedlar etal. (1990) 




NGC 3079 


0.044 


125 




65 


VLBA 5/8/22 GHz 


Trotter etal. (1998) 


75 




310 




60 




VLA 1.4/4.9 


Duric et al. (1983) 




NGC 4151 


0.116 


257 




3 


VLBA 1.6,5 GHz 


Ulvestad et al. (1998) 


140 




104 




80 




VLA 5 GHz 


Ho & Ulvestad (2001) 




NGC 4258 


0.029 


3 




8 


VLBA 22 GHz 


Herrnstein et al. (1997) 


60 




125 




355 




VLA 1.5 GHz 


Cecil et al. (2000) 




NGC 7212 


1.37 


257 




52 


global VLBI 5 GHz 


Vir Lai (2001) 


133« 




609 




25 




VLA 5 GHz 


Vir Lai (2001) 




NGC 7469 


0.844 


28 




35 


global VLBI 5 GHz 


Vir Lai (2001) 


35 
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Name 


Res. 


0(pc) 


(kpc) 


AO 


Array 


Refs. 


(opt.) 




(pc) 


(deg) 


(deg) 


(deg) 






deg 




193 




353 




VLA 5 GHz 


Vir Lai (2001) 




NGC 7674 


2.80 


274 




22 


VLBA+Y27+AO 1.4 GHz 


Momjian et al. (2003) 


64" 




178 




296 




MERLIN 1.7 


Momjian et al. (2003) 




inzw2 


0.260 


303 




63 


VLBA 15/43 GHz 


Erunthaler et al. (2000) 






3460 




240 




VLA 1.5 GHz 


Branthaler, priv. comm. 




Ark 564 


0.598 


282 






global VLBl 5 GHz 


Vir Lai (2001) 


28 


IC 5063 


3.30 


66 




49 


LEA 2.3 GHz + LEA HI 


Oosterloo et al. (2000) 






242 




115 




ATCA 8.3/1.4GHZ 


Morgantiet al. (1998) 




NGC 2110 


0.151 


350 




16 


VLBA 8.4 GHz 


Mundell et al. (2000) 






240 




186 




VLA 1.4/5 GHz 


UIvestad& Wilson (1983) 




NGC 3147 


0.219 


317 






VLBA 1.6/2.3/5/8.4 GHz 


Ulvestad & Ho (2001b) 


65 


NGC 5506 


0.815 


73 




76 


EVN 1.6/5 GHz 


Middelberg et al. (2004) 


1 




359 




177 




VLA 4.9 GHz 


Colbert et al. (1996) 




NGC 7682 


0.996 


125 






global VLSI 5 GHz 


Vir Lai (2001) 


68* 
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Figure 5.5: Histograms of 
Seyfert galaxy literature data. 
Top: Histogram of spectral 
indices with a bin width of 
0.5. Bottom: Histogram of an- 
gles between pc and kpc-scale 
structure with a bin width of 
15°. 
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Chapter 6 



Fast Frequency Switching 



The invention of Very Long Baseline Interferometry was mainly driven 
by the need for higher angular resolution radio observations. Single-dish 
observations would never have approached the resolutions achieved by 
optical telescopes owing to the long wavelengths. However, once VLBI was 
available, it was straightforward to progressively increase the resolution by 
increasing the observing frequencies. In the recent few years, observations 
at 86 GHz have evolved from the experimental stage to a commonly used 
observing technique which achieves the highest angular resolutions possible 
today, and first detections at even higher frequencies have already been made 
(Krichbaum et al. 2002). The scientific benefits are undoubted: even in the 
closest objects, only millimeter- VLBI can resolve the linear scales on which 
extragalactic radio jets are launched, a process yet not understood. 86 GHz 
VLBI observations also contribute to investigations of internal jet structure, 
composition and coUimation (e.g., Doeleman et al. 2001, Gomez et al. 1999). 

Unfortunately, higher frequency observations involve a number of se- 
rious problems: the sources are usually weaker because the emissivity of 
optically thin synchrotron sources drops as v^°-^, the efficiencies of typical 
centimetre radio antennas drop to about 15 % because the surface rms is 
« A/10, the receiver performances become disproportionately worse because 
the amplifiers contribute more noise and the atmospheric coherence that 
limits the integration time decreases as 1/v. This last problem, however. 
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6.1. Introduction 



is not a hardware property and would be solved if the atmospheric phase 
fluctuations could be calibrated. In fact, the so-called phase-referencing first 
demonstrated by Alef (1988) is now an established calibration technique in 
centimetre VLB I to lower the detection thresholds to the sub-mJy level and 
to accurately determine the positions of sources in VLBI observations. In 
phase-referencing, a sttong calibrator source is frequently observed (every 
few minutes, depending on observing frequency) to calibrate the visibility 
phases of the target source integrations, i.e., the telescopes cycle between 
the target source and the phase calibrator source. In millimetre VLBI the 
technique is not commonly used, although a successful proof of concept exists 
(Porcas & Rioja 2002). Problems arise from the need of a suitable, strong 
phase calibrator in the vicinity of target source and from a combination of 
short atmospheric coherence and telescope slewing times. 

In this chapter, I describe a new phase-referencing technique for VLBI 
observations at high frequencies, fast frequency switching. 



6.1 Introduction 

The main source of phase noise in VLBI observations at frequencies higher 
than 5 GHz is turbulence in the troposphere causing refractive inhomo- 
geneities. These path-length variations are non-dispersive, and one can 
self-calibrate the visibility phases at one frequency and use the solutions to 
calibrate visibilities at another frequency after multiplying the phases by the 
frequency ratio, r, provided that the lag between the two measurements does 
not exceed half the atmospheric coherence time (Fig. 6.1). This is possible 
with the VLB A because frequency changes need only a few seconds. After 
multiplying the phase solutions by the frequency ratio and applying them 
to the target-frequency phases, there remains a constant phase offset, AO, 
between the signal paths at the two frequencies, which must be calibrated. It 
can be monitored with frequent observations of achromatic, strong calibrators, 
and must be subtracted from the high frequency visibility phase. 

Following the notation from Eqs. 2.4ff, the observed visibiUty phases using 
this calibration scheme can be described as 
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Mh) = Uh) + &2) + 0p^fe) + 0^t(^2) + €M + C'n(f2) (6.1) 

^,(f3) = 0r(?3) + C;(f3) + Csit3) + <^>^t(f3) + «/>atm(?3) + ^(^3)- 

Here, the indices Vr and Vt denote terms of the reference and target fre- 
quency, respectively, and (pj- and are the true visibiUties at the reference and 
target frequency, respectively. Interpolation of the reference frequency terms 
to the times of the target frequency scans yields 

Mt2) = $r{t2) + &2) + ^p^(?2) + &{t2) + &(?2) + &t2). (6.2) 

Self-calibration at the reference frequency is then used to obtain a source 
model, and self-calibration of the visibiUties after dividing by this model yields 
the required phase corrections. These correction terms are interpolated to the 
times where the target frequency was observed, and scaled by the frequency 
ratio, r. In practice, the phase solutions derived in this process comprise con- 
tributions from all five sources of error, and so do the scaled phase corrections. 
The difference between the target frequency visibilities and the interpolated 
reference frequency visibilities then is 

In this equation, is zero because a model of the source has been sub- 
tracted. The instrumental phase offset is constant and can be determined from 
calibrator observations, and so the term (0-^^'^, — r^j^'J is known and can be re- 
moved. The antenna position errors also cancel out because this type of error 
causes a shift of the source on the sky, and to compensate for it requires a phase 
correction which scales with frequency. E.g., assume an antenna position error 
of 1 cm along the line of sight on a 1000 km baseline. The error shifts the fringe 
pattern on the sky by (0.01 m/lOOOkm) = 10^^ rad = 2.06 mas. At a frequency 
of 15 GHz (A = 2.00cm), the fringe spacing on the sky is A/Z) = 4. 13 mas, and 
hence a phase correction of 179° is required to account for the error. At a fre- 
quency of 45 GHz (A = 0.67 cm), three times higher than the other frequency, 
the fringe spacing is 1.38 mas, and a phase correction of 537° is required to 
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compensate for the error, three times more than at the lower frequency, and 
so (<^'^t ~ '"&) = 0- Similar to the antenna position error, the source posi- 
tion error scales with frequency when measured in units of turns of phase, 
but a frequency-dependent core shift, as predicted by conical jet models (e.g., 
Lobanov 1998), modulates {(f)^^ — r^pg^) on each baseUne with a sinusoid. 
The period of the sinusoid is 23*^56'", its amplitude depends on the magnitude 
of the shift and its phase depends on the direction of the shift relative to the 
baseline direction. It is zero only in the abscence of a core shift, and therefore 
cannot be neglected. The atmospheric phase errors also scale with frequency, 
and hence ((/>^ — r^^) = 0, but the ionospheric phase noise is difficult to de- 
termine; something which we, unfortunately, did not discover beforehand the 
observations presented here. It can probably be determined using interspersed, 
wide-band scans at a low frequency, e.g., in the 1.4 GHz band, where the ef- 
fect is strong. Similar to Eq. 2.7, the remaining terms describe the difference 
between the target frequency phase and the interpolated reference frequency 
phase as the target frequency visibihty phase plus the position offset and an 
interpolation term, (pim: 

K - h = ^i + (<os - '•^p^'os) + fet- (6.4) 

The result is the high-frequency visibilities phase-referenced to the source 
at low frequency, and so the technique can measure the position shift of AGN 
cores with frequency and can prolong coherence. The coherence time can 
in principle be prolonged indefinitely, however, in our project, unmodelled 
ionospheric path length changes limited the coherence to half an hour in the 
worst case. Thus, in the following sections, we use fast frequency switching 
to calibrate the short-term phase fluctuations and then use self-calibration at 
the target frequency with a half-hour solution interval to remove the long- 
term phase drifts. The extension from the 30 s atmospheric coherence time 
to a coherent integration time of 30 min yields a large sensitivity improvement. 

A similar observing strategy has been developed by Kassim et al. (1993) 
for the VLA, who used scaled phase solutions from 330 MHz to calibrate 
simultaneously observed 74 MHz data. In this case, the dominant source of 
phase errors was ionospheric path length changes. Using their new calibration 
technique, they were able to increase the coherent integration time from 
< 1 min to > 10 min and to make 74 MHz images of several radio sources. 
Fast frequency switching is also being considered as a standard calibration 
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mode for ALMA in the future (D'Addario 2003). 



14.375 vs. 43.125 GHz Visibility Phases 

Upper panel: Hancock, lower panel: Kitt Peak 
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Figure 6.1: Demonstration of the scalability of phase solutions. 43 GHz fringe-fitted phase solutions 
(solid lines) on 3C 273 compared to 15 GHz fringe-fitted phase solutions multipUed by the frequency ratio 
(dashed lines) from Hancock (upper panel) and Kitt Peak (lower panel). Ten minutes of data are shown, 
and the phases foUow each other very well. 



6.2 Observations 

We observed a pilot project (BR073) on January 5, 2002, to attempt to measure 
the diameter of the AGN in M81 (M81*) at 86 GHz. We obtained only very 
weak detections of M81*, but we learnt a number of important points about 
designing fast-frequency switching experiments (Middelberg et al. 2002). We 
incorporated those considerable improvements in the observing strategy for 
our next project (BM175C), which we describe here. 

We observed NGC4261 as a fast frequency switching target and 3C273 
and 3C279 occasionally for the inter-frequency offset and to test the tech- 
nique on strong sources on May 5, 2003 (project BM175C). A previous run 
(BM175B) had to be re-observed because of high fringe rates due to weather 
and a number of equipment failures. BM175C was observed in a period of 
superb weather, using 256 Mbps to record a bandwidth of 64 MHz with 2-bit 
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Table 6.1 : Summary of project BM175 observed on May 5, 2003. The sources observed are in column 
(1), the scan duration is in column (2) and the frequency pairs in GHz are in colirain (3). An a indicates 
fringe finder scans during which no frequency switching was done. 
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sampling. We recorded LCP only, the data were divided into 8 IFs with 8 MHz 
bandwidth, each of which were subdivided into 64 spectral channels 125 kHz 
wide. The correlator integration time was one second to allow monitoring of 
the phases with high time resolution. All antennas performed well, except for 
Fort Davis, where a receiver problem caused complete loss of 43 GHz data. A 
summary of this observing run is given in Table 6.1. 

Several considerations influenced the experiment design: 

• Frequency Choice The target frequency should be an integer multiple 
of the reference frequency to avoid having to unwrap phase wraps. For 
example, if the frequency ratio, r, is the non-integer value of 2.5, and 
the reference frequency phase wraps from 359° to 0°, then the scaled 
target frequency phase will jump from 897.5° (= 177.5°) to 0°, intro- 
ducing a phase jump of 177.5° into the calibration phase for a 1° phase 
change at Vref. In contrast, choosing r = 2.0, when the reference fre- 
quency phase wraps from 359° to 0°, the scaled target frequency phase 
changes from 718° (= 358°) to 0°, corresponding to a change of 2° for 
a 1° phase change at Vref. Hence, for integer values of r, phase wraps at 
the reference frequency introduce integer multiples of 360° at the target 
frequency and so have no effect. 

We chose a reference frequency of 14.375 GHz since the third and sixth 
harmonics at 43.125 GHz and 86.25 GHz lie within the VLBA receiver 
bands. For convenience, we will refer to these frequencies as "15 GHz", 
"43 GHz" and "86 GHz", respectively. In this document, two consecu- 
tive integrations at Vret and Vt will be called a "cycle", each integration 
of which is called a "half-cycle", and a sequence of cycles on the same 
source is called a "scan". Long (several minutes), continuous integra- 
tions on a single source at one frequency, e.g. fringe-finder observations, 
will also be called "scans". 

• Integration Times We chose a cycle time of 50 s, of which 22 s were 
spent at the reference frequency of 15 GHz and the remaining 28 s were 
spent at the target frequency, either 43 GHz or 86 GHz. An average time 
of 7 s per half-cycle was lost in moving the subreflector between the feed 
horns, resulting in net integration times of 15 s at Vref and 21 s at Vt. The 
integration times are a compromise between source brightness, antenna 
sensitivity and expected weather conditions. This setup yielded a 5 (7 
detection limit of 89 mJy in 15 s at 15 GHz for the VLBA. 
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• Calibrator Scans Calibrators must be observed frequently to monitor 

the phase offset A<I> between the two frequencies. An important con- 
straint is that the calibrators must be achromatic, i.e., they must not have 
their own frequency-dependent core shifts. We included adjacent scans 
on two different calibrators to measure AO before and after target source 
observations. Five of the calibrator scans were twice as long to provide 
more time for tests with strong signals. 



6.3 Data Reduction 

6.3.1 Standard Steps 

Data reduction was carried out in AIPS. A calibration table entry was gen- 
erated every 4.8 s to provide high temporal resolution. The amplitudes were 
calibrated using Tsys and gain measurements provided by the VLBA's auto- 
mated calibration transfer, and amplitude corrections for errors in the sampler 
thresholds were performed using autocorrelation data. Phase corrections for 
parallactic angles were applied and a simple bandpass correction was derived 
at each frequency from one of the fringe finder scans. 

The VLBA's pulse calibration system did not deliver usable data to cali- 
brate the phase offsets between the IFs because the frequencies were changed 
too quickly. The pulse calibration system has a default integration time of 
10 s, plus one second for readout. If any part of that integration time is during 
a part of the scan that is flagged by the online system, then the integration is 
not accepted. Thus, one looses the first two integrations because the online 
system conservatively flags 10 s to lis after a frequency change, and there is 
an overlap between the flagged time and the second integration. During the 
whole experiment, only ~ 10 useful pulse-cal measurements were recorded 
per station at 15 GHz (half-cycle time 22 s), and ~ 100 at 43 GHz (half-cycle 
time 28 s), compared to the number of half-cycles at these frequencies of 
~ 630 and ~ 230, respectively. We therefore used self-calibration on the 
same fringe finder scans as were used for bandpass calibration to correct for 
instrumental delays and phase offsets. These offsets were found to be stable 
over the experiment. 
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6.3.2 Ionospheric Correction 

Although the frequencies used in the project are quite high, ionospheric effects 
can not be neglected and will prevent a successful phase transfer if uncorrected. 
A typical ionospheric delay at a frequency of Vi = 100 MHz is 0.1 /^s (10 
turns of phase), the exact number depending on the time of day. The delay 
scales with frequency as v^^, and so for a delay at 100 MHz of 0.1 jUs, the 
ionospheric delays at 14.375 GHz, 43.125 GHz and 86.250 GHz are 4.84 ps, 
0.54 ps and 0.13 ps, respectively, corresponding to 25°, 8.4° and 4.0° of phase. 
In contrast, the path length through the troposphere is non-dispersive and so 
the linear phase versus frequency scaling law used by fast frequency switching 
caimot correct phase changes that are induced by the ionosphere. Ionospheric 
phase changes have much longer time scales than tropospheric changes, and 
they can be calibrated before fringe-fitting when the electron content of the 
ionosphere along the line of sight is known. The AlPS task TECOR can use 
maps of ionospheric total electron content derived from GPS data to calculate 
phase and delay corrections. Unfortunately, the error in these maps can be 
quite high, up to 20 % when the TEC is as high as a few tens of TEC units 
(1 TECU=10^^electronsm^^), and up to 50 % or higher when the TEC is of 
the order of a few TEC units. We have used the TEC maps produced by the 
Center for Orbit Determination in Europe (CODE^) to calibrate the effects of 
the ionosphere. We found that these maps yielded better results than those 
from the Jet Propulsion Laboratory (JPL), i.e., the residual phase errors after 
scaling were smaller. We do not know whether this finding is coincidental or 
the CODE maps generally give better results. 

6.3.3 Flagging 
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Table 6.2: Flagging times in seconds at each station and frequency applied to the beginning of each half- 
cycle with QUACK. 



From the start of a new half-cycle 1 s to 2 s are required to set up the 
electronics and 5 s to 8 s are required to position the subreflector. Data 

'http : / / www . aiub . unibe . ch/ ionosphere . html 
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during that time should be flagged. The flags generated by the VLBA onUne 
system turned out to be too conservative, flagging more data than needed 
and sometimes leaving only 10 s of data per integration at either Vref or 
Vt. We used amplitude vs. time plots of caUbrator scans to estimate the 
switching times for each antenna separately. We found that the times needed 
to move the subreflector are frequently asymmetric (e.g., moving it from Vt 
to Vref took more time than moving it back). Furthermore, they vary slowly 
with time, presumably due to elevation-dependent gravitational force on the 
subreflectors. Compromise flagging times for each antenna and frequency 
at the start of each half-cycle were applied using the AIPS task QUACK 
(Table 6.2). 

The effect of the flagging is illustrated in Fig. 6.2. Amplitudes need 2 s to 
3 s longer to reach their final values than do the phases. Thus, the visibility 
phases are not very sensitive to even large errors in subreflector rotational 
position. The repeatability of the phases from cycle to cycle shows that 
positioning the subreflector along the optical axis is repeatable to < 5° of 
phase at 15 GHz, which is much less than other sources of phase error in fast 
frequency switching. 

Occasionally, a subreflector does not move at all, resulting in a half-cycle 
of noise. These events have been flagged by calculating the duration of the 
entries in the Ust of flags provided by the VLBA and flagging those cycles that 
have flagging entries for longer than 30 s and therefore exceeded the half-cycle 
time. 

6.3.4 Fringe-Fitting and Phase Scaling 

We fringe-fitted the 15 GHz data using the AIPS task FRING, with an SNR 
threshold of 5 and delay and rate search windows of 20 ns and 50mHz, 
respectively. As NGC4261 has an extended, double-sided jet structure at 
15 GHz, we made a 15 GHz image that we used as a source model in a second 
run of FRING, so that the phase solutions did not contain structural phase 
contributions (Fig. 6.10). The solution interval was set to 1 min, yielding one 
phase, delay and rate solution per half cycle. The detection rate was ~ 90%. 
The 15 GHz solution (SN) table was written to a text file with TBOUT to do 
the phase scaling outside AIPS. We have written a Python program (FFSTG, 
the Fast Frequency Switching Table Generator) that processes a SN table in 
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the following way. 

First, a series of timestamps is generated from each pair of consecutive en- 
tries in the input table such that they coincide with the Vt half-cycles. Second, 
from each pair of consecutive Vref phase solutions, a solution is interpolated 
for the new timestamps consisting of a phase, a phase rate derived from the 
Vref phase solutions and the time interval between the Vref scans, and an inter- 
polated delay. Both the phase rate and the delay do not need to be scaled by r 
because the rate is stored in a frequency-independent format (in units of ss~^) 
and the delay is non-dispersive. Third, the interpolated solution is scaled by 
the frequency ratio, r, yielding Vt phase solutions for the times at which the 
source was observed at Vt. 

We did not use the phase rates derived by FRING because each of those 
was derived from one half-cycle of 22 s length, whereas the phase rates interpo- 
lated from two consecutive half-cycles as described above used two half-cycles 
separated by 50 s and therefore have much better SNR. A further advantage of 
deriving phase rates from pairs of phase solutions is that each Vref phase so- 
lution is used both in the determination of the phase rate to the preceding and 
the succeeding Vref solution, causing a smoothing of the phases with time and 
reducing the effect of outliers. Finally, the interpolated phase, phase rate and 
delay solutions are stored in an output table together with the Vt frequency 
ID. FFSTG provides an interface to Gnuplot to plot and inspect the input and 
output phases. One can select regions of interest and check the results of the 
scaling. 

The table is imported to AlPS using TBIN and can be used to update the 
most recent calibration table at the target frequency. 



6.4 Results 

6.4.1 43 GHz 

NGC 4261 was detected on most baselines at all times after scaling the 15 GHz 
solutions to 43 GHz, with correlated flux densities of 30mJy (800 MA) to 
160 mJy (SOMA). The 43 GHz half-cycle average visibilities on baseUnes 
to BR and LA are shown in Figs. 6.4 and 6.6. The short-term fluctuations 
introduced by the troposphere are almost perfectly calibrated, but residual 
phase drifts remain on longer timescales, especially at the beginning and at 
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the end of the experiment, when the sun was setting or rising. 



Structure Functions 

Structure functions of the 43 GHz visibility phases from two 25 min scans on 
NGC4261 are shown in Fig. 6.8. Solid black lines show raw data observed 
at night with delay calibration only, and solid red lines show the same data 
calibrated with scaled-up phase solutions from fringe-fitting at 15 GHz. The 
dashed lines show the same stages of calibration from a 25 min scan observed 
during the evening. 

The structure functions constructed from the calibrated phase series show 
a residual phase noise on the shortest timescales (60 s) of 50° during sunset 
and 33° at night (taking the median of the phase noise of all baseUnes). The 
phase noise increases toward longer timescales, probably due to errors in the 
ionospheric models. Evidence for the long-term phase noise being dominated 
by ionosphere is as follows. 

First, applying an ionospheric TEC correction significantly improved the 
coherence, as one can see by comparing Figs. 6.6 (lower panel) and 6.9. Sec- 
ond, the residual phase errors decreased at night. After sunset at 19:56 local 
time in Los Alamos (=02:56 UT), representative for the south-western anten- 
nas, the phase errors rapidly dropped. 

We used FRING on the 43 GHz data calibrated with fast frequency 
switching with a solution interval of 30 min, to self-calibrate the data and 
so to remove the residual long-term phase drifts. The structure functions 
after applying this correction are shown as blue lines in Fig. 6.8; the 
long-term phase noise is reduced as expected, and the phase noise on 60 s 
timescale is sUghtly lower, with a median of 44° during sunset and 3 1 ° at night. 



Expected Phase Noise 

The expected phase noise in the visibilities calibrated with fast frequency 
switching consists of three parts: (1) thermal phase noise at the reference fre- 
quency scaled by the frequency ratio, (2) thermal phase noise at the target 
frequency and (3) tropospheric phase changes the two integrations. One can 
estimate the expected rms phase noise from sensitivity calculations as follows. 
The 1 a thermal noise on a baseline is given by 
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1 SEFD 
A5 = — X , (6.5) 

TJs V2 X Av X T 

(Walker 1995), where T7s is a scaling factor, SEFD is the antenna's system 
equivalent flux density, Av is the bandwidth and T the integration time. The 
efficiency factor depends on accurately known effects like losses due to sig- 
nal quantisation, clock drifts, and others, but also on time-variable, unknown 
effects hke antenna pointing errors and moisture on the feeds and reflecting 
surfaces. It therefore has to be determined empirically, and for the VLBA usu- 
ally is assumed to be 0.5 for 1-bit sampling. The observations presented here 
were done using 2-bit sampling, increasing the sampUng efficiency from 0.637 
to 0.881, or by 38 %, and T7s then is 0.69. 

According to the VLBA Observational Status Summary^, the SEFD 
for a VLBA dish is 550 Jy at 15 GHz, 1436 Jy at 43 GHz and 5170 Jy at 
86 GHz. In our observations, T was 15 s at 15 GHz and 21s at 43 GHz 
and 86 GHz, and Av was 64 MHz. One therefore expects thermal noise 
levels of 17.8 mJy, 39.2 mJy and 141 mJy at 15 GHz, 43 GHz and 86 GHz, 
respectively. At 15 GHz, NGC4261 has a correlated flux density, 5, between 
60mJy and 300mJy, depending on baseline length, and we do the following 
calculations assuming either 80mJy, representative for long baselines, or 
200 mJy, representative for short baselines (the results of which we give in 
brackets). 



(1) On a single baseline, the expected SNR of a detection at 15 GHz is 
S/AS = 80mJy/17.8mJy =4.49 (11.2) when averaging over the band. Fringe- 
fitting, however, uses all baselines to one particular antenna to derive a phase 
correction, and so the SNR is increased by \/N, where A'^ is the number of 
baselines. At any given time, N ^S, so that the SNR of a detection at 15 GHz 
increases to 12.7 (31.7), corresponding to a phase error of 1 rad/SNR = 4.51° 
(1.81°). These errors are scaled by the frequency ratio to 13.5° (5.43°) at 
43 GHz and 27.1° (10.9°) at 86 GHz. 

(2) Assuming that NGC 4261 has a compact flux density of 100 mJy at both 
frequencies, the thermal noise contributions are 21.4° and 54.7°, respectively. 
Adding those in quadrature to the scaled-up thermal rms phase noise from 

^http : / / www . aoc . nrao . edu/vlba/ obstatus/ obssum/ obs sum. html 
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15 GHz yields 25° (22.1°) at 43 GHz and 56.3° (55.0°) at 86 GHz. 

(3) We estimated the tropospheric phase noise within the switching cycle 
time using structure functions of 3C273. On 3C273, the thermal noise con- 
tributions at 15 GHz and 43 GHz are 0.04° and 0.4°, respectively, so that the 
visibihty phases are essentially free of thermal noise and any phase changes 
during and between the half-cycles are due to changes in the troposphere. We 
found the median rms phase noise after self-calibration with a 30 min solution 
interval to remove the residual long-term phase drift to be 13.3° at 43 GHz, or 
26.6° at 86 GHz. 

Adding those three constituents in quadrature yields 28.3° (25.8°) at 
43 GHz and is in good agreement with the measured rms phase noise of 31°. 



Coherence 

Another way of illustrating the performance of the corrections is with coher- 
ence diagrams (Fig. 6.11). For the plots shown, NGC4261 43 GHz data from 
all basehnes to three stations, BR, KP and LA, have been used. All data were 
calibrated using fast frequency switching (no subsequent fringe-fitting), and 
were averaged over the band and in time to obtain one visibility per baseline 
and half-cycle. Assuming unit amplitude, the vector sum of these visibilities 
has been computed over progressively longer time intervals and the length of 
the resulting vector has been normalized to one by dividing by the number 
of visibilities. All of these sums have then been averaged for each baseline 
separately to yield one coherence measurement per baseline and coherent 
averaging time interval. The left panels show data from a 25 min scan centred 
on 01:06 UT (early evening) and the right panels show data from a 25 min 
scan centred on 03:40 UT (at night). Bold lines show the average over all 
baselines and the error bars show the 1 a scatter of the baseline average. The 
coherence achieved at night with fast frequency switching is significantly 
better than that achieved during the early evening. 

Before making an image from the 43 GHz visibilities calibrated with fast 
frequency switching, the residual phase offsets and phase rates were removed 
using FRING with a solution interval of 30 min so that one solution per scan 
was obtained. The resulting dirty image (Fig. 6.12) has a peak flux density 
of 79mJybeam^' and an rms noise of 4.4 mJy beam" \ yielding a dynamic 
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range of 18:1. After several cycles of phase self-calibration with a solution 
interval of 30 s and one cycle of amplitude self-calibration with a solution in- 
terval of 12h, the peak flux density was 95mJybeam^' and the rms noise was 
0.78 mJy beam" ^, so the dynamic range improved to 122:1 (Fig. 6.13). The 
theoretical rms noise at 43 GHz was expected to be 0.46mJybeam~^. Unfor- 
tunately, the ionospheric phase errors and the need to run FRING to caUbrate 
them did not allow a core shift measurement to be made. 
The image dynamic range can be estimated using 



DNR = 




(6.6) 



(Perley 1986), where M is the number of independent measurements, is 
the number of antennas and (^nns is the rms of the phase noise. The number 
of independent measurements is given by the number of 43 GHz half-cycles, 
corresponding to the basis on which 0rms was calculated. With seven 43 GHz 
scans, 30 half-cycles per scan, and multiplying by the number of baselines that 
gave good data (A^ = 8, so 28 baselines), M = 5880. The average of 0rms during 
the evening and at night was (44° -|-31°)/2 = 38°, and the expected dynamic 
range follows to 16:1. This is in excellent agreement with what we measured 
in the first 43 GHz image (18:1). Note that this agreement is independent of 
the half-cycle times. Doubling the cycle time decreases (^rms by \/2, but M is 
decreased by two at the same time, and the result remains unchanged. 

6.4.2 86 GHz 

Following the same data reduction path as for the 43 GHz data, we obtained 
good detections of NGC4261 at 86 GHz on baselines among the four stations 
FD, KP, LA, and PT. Uncorrected and corrected visibility phases are plotted 
in Figs. 6. 14, structure functions and coherence plots from a 25 min scan are 
shown in Figs. 6.15 and 6.16 and images are shown in Fig. 6.17. The peak 
flux density in the upper panel of Fig. 6.17 is 59.3 mJy. This is, we believe, 
the first detection of NGC 4261 with VLBl at 86 GHz. It is also probably the 
weakest continuum object ever detected with VLBl at this frequency, because 
it is less than the 85mJy in the 86 GHz image of 1308 + 328 by Porcas & 
Rioja (2002), using conventional phase-referencing to a caUbrator 14.3' away. 
With only delay calibration applied, the median rms phase noise of the best 
25 min scan is 104°, after applying the scaled 15 GHz phase solutions is 70° 
and after self-calibration with a 30 min solution interval is 80°. The increase 
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in rms phase noise after removal of phase rates is due to the small number 
of stations in the sample. The expected thermal phase noise is the quadrature 
sum of the scaled-up 15 GHz noise and the 86 GHz noise, the former of which 
we estimated to be 27.1° on long and 10.9° on short baselines, and the latter, 
assuming a correlated flux density of 100 mJy, is 54.7°. The quadrature sum 
is 61.0° on long and 55.7° on short baselines, and adding the tropospheric 
phase errors of 26.6° is 66.6° and 61.7°, in agreement with the measured noise 
levels. The rms noise in the final 86 GHz image is 7.4mJybeam ' , compared 
to an expected thermal noise of 3.7mJybeam~^. The dynamic range in the 
first image was 7:1, and assuming seven antennas and an average rms phase 
noise of 64° gives an expected dynamic range of 9:1. Thus, the results of fast 
frequency switching agree very well with the theoretical expectations. 

6.5 Summary 

• Fast-frequency switching can be used to calibrate tropospheric phase 
fluctuations if the switching cycle time is shorter than the atmospheric 
coherence time. The accuracy is limited by tropospheric phase changes 
between and during the reference frequency integrations. 

• Insufficient knowledge about the ionosphere's total electron content 
(TEC) has prevented calibration of the inter-band phase offset and hence 
making a pure phase-referenced image without using selfcal was not 
possible. This also prevented the detection of a core shift. Current global 
TEC models derived from GPS data have errors that are too large to suf- 
ficiently calibrate the ionospheric component of the phase changes. A 
possible solution could be to insert frequent (every lOmin) scans at L- 
band with widely separated IF frequencies to derive the effective TEC 
on the line of sight to the target source. 

• The VLBA pulse calibration system might be able to monitor instru- 
mental phase changes if it is modified to cope with the unusually short 
integration times. 

• The instrumental phase offset seen on the LA-PT baseline was stable to 
<1 rad over 10 h, and it seems to be sufficient to determine the offset a 
few times throughout the experiment. Measuring it every 30 min like we 
attempted therefore is not necessary. 
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• Truly simultaneous observations at two bands would reduce the residual 
phase errors even further because tropospheric phase changes within the 
switching cycle time would be entirely removed. 

The primary use of fast frequency switching is to detect sources that are 
too weak for self-calibration at the target frequency within the atmospheric 
coherence time, but can be reliably detected at a lower frequency. The 5 a 
detection limit of the VLBA at 43 GHz within 120 s, with 64 MHz bandwidth 
and 2-bit sampling, and neglecting coherence loss due to tropospheric phase 
changes, is 81.9 mJy (i.e., a thermal noise level of 16.4 mJy on a single base- 
line). The 1 a noise level using fast frequency switching with half-cycle times 
of 22 s at 15 GHz and 28 s at 43 GHz (yielding net integration times of 15 s and 
21 s, respectively), after one cycle is 66.2 mJy, of which 39.1 mJy is thermal 
noise in the raw 43 GHz visibility, 53.4 mJy is thermal noise in the phase 
solutions after scaling from 15 GHz, and neglecting noise due to tropospheric 
phase changes. The 1 o noise level after 120 s of fast frequency switching (i.e., 
2.4 cycles) is 42.7 mJy. The noise level of fast frequency switching reaches 
that of a conventional 120 s integration at 43 GHz after 815 s (13.6 min), i.e. 
16.4 mJy. Any longer integration with fast frequency switching then yields 
detection thresholds that are out of reach with conventional methods. 

Fast frequency switching, however, is limited by the source strength at the 
reference frequency, and 100 mJy at 15 GHz is, from our experience, a reason- 
able minimum required for a successful observation. However, many sources 
exist that meet this criterion. The 2 cm survey by Kellermann et al. (1998), for 
example, detected 130 sources with peak flux densities of 100 mJy or more; all 
of them can be observed at 86 GHz using fast frequency switching. Another es- 
timate of the number of new sources made accessible to 86 GHz VLBI by fast 
frequency switching can be made by means of the logA'^ — log 5 relation. We 
assume that the sources observable with conventional 86 GHz VLBI are dis- 
tributed homogeneously in the universe, i.e., lowering the detection limit does 
not yield the detection of a new population, so the number of sources above a 
given flux density limit, S, increases as S~^'^. We further make the conserva- 
tive assumption that these sources have a = between 15 GHz and 86 GHz, 
because they are compact. The current 5 a detection limit of the VLBA at 
86 GHz is 600 mJy (thermal noise of a 30 s integration with 128 MHz band- 
width on a single baseline). In fringe-fitting, all N antennas are used to derive 
a phase solution, so the noise level is reduced by V^, and with N = 1, assum- 
ing one VLBA station capable of 86 GHz observations is missing at any time. 
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the 5 (7 detection limit is 230 mJy. Fast frequency switching allows one to ob- 
serve any compact source with S^ecHz > lOOmJy, a factor of 2.3 fainter than 
with conventional methods, and one expects eight times more sources brighter 
than 100 mJy than there are brighter than 230 mJy, following the logA'^ — log 5 
relation. A little less conservative estimate would be to allow for spectral 
indices as steep as a = ^0.3, so that sources with 5i5ghz = lOOmJy have 
SgeoHz > 60mJy. The number of observable sources would then increase by a 
factor of 29 of the number accessible to conventional techniques. 

The absolute number of presently observable sources is much more dif- 
ficult to determine due to selection effects in the various radio source cata- 
logues. However, we can make a simple estimate using the recent 86 GHz 
VLBl survey by Lobanov et al. (2002). They aimed at observing more than 
100 sources collected from the literature that had an expected compact flux 
density of SgecHz > 300 mJy. Applying the fast frequency switching sensitiv- 
ity estimates then yields between 800 and 2900 sources observable with the 
VLBA at 86 GHz. 

Furthermore, the proposed capability of fast frequency switching to detect 
core shifts in AGN can supplement jet physics with observational constraints 
which are otherwise difficult to obtain. 

6.6 A Recipe 

This section summaries the experiment design and the data reduction for those 
who want to use fast frequency switching for their observations. 

6.6.1 Experiment Design 

• Frequencies The choice of target and reference frequencies needs to be 
made with careful attention to the following issues. The calibration of 
Vt critically depends on the ability to self-calibrate the Vref data. Ob- 
servers should check whether the Vref provides sufficient baseline sen- 
sitivity within the available Vref integration time to yield reliable detec- 
tions. Second, the frequency ratio r needs to be integer. As this can be 
true for only an infinitely small bandwidth, one should select the fre- 
quencies such that the centres of the Vt and Vref bands have an integer 
ratio (this condition also prefers 2-bit sampling over 1-bit sampling be- 
cause of the smaller bandwidth at the same recorded bitrate). If the band- 
width is 64 MHz and the centre frequencies are Vref = 14.375 GHz and 
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Vt = 43.125 GHz or Vt = 86.250 GHz, the frequency ratios at the lower 
edge of the bands (where the deviation from an integer ratio reaches a 
maximum) are 

_ 43.125GHz- 32MHz _ 
14.375 GHz 

86.250GHz -32MHz 

r = TT^rrTTT^TT = 5.9978. 

14.375 GHz 

Thus, for each phase turn at 15 GHz, one accumulates a phase error 
of 0.79° at 86 GHz, and after 72 turns the error would be one radian. 
However, averaging reduces this error because it has opposite sign at the 
other edge of the band. Furthermore, the phase does normally not wind 
constantly in one direction throughout the observation. 

• Ionospheric corrections Given the results from the data presented in 
this memo, we think that regular (every lOmin), short (1 min) scans at 
L-band with widely separated IFs should be observed during the obser- 
vations to derive the ionospheric dispersive delay. This should allow one 
to separate the ionospheric and tropospheric phase errors and to obtain 
coherence over several hours even during dusk and dawn. We have not 
tested the procedure but have estimated the effect on a data set from a 
different project (BM170), where the flux density on the longest base- 
lines was 50 mJy. We obtained one phase solution per IF every 30 s using 
fringe-fitting, and the delays were in good agreement with our expecta- 
tions, although the effect was small due to the small bandwidth that was 
covered by the IPs. 

• Switching cycles The switching cycle times depend very much on the 
weather conditions. As they are unknown, one should apply for dy- 
namic scheduling so that a certain level of "good" weather can be as- 
sumed. During the observations presented here, the coherence times 
were mostly 5 min to 10 min at 43 GHz, occasionally decreasing to 
2 min. We chose a cycle time of 50 s, short enough to follow the phase 
at Vref most of the time. 



Calibrator scans One should insert regular scans on strong, nearby cal- 
ibrators to be able to check the phase transfer from Vref to Vt during 
postprocessing and to measure the inter-frequency instrumental offset. 
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The calibrators must be strong enough for self-calibration at both Vt and 
Vref within the respective integration times. At least 5 to 7 of those cal- 
ibrator observations should be lOmin to ISmin long to estimate the at- 
mospheric coherence times and to help debug the calibration. There is 
the danger of transferring a possible core shift to the target. This can be 
avoided if at least two, but preferably more, calibrators are used to check 
the calibration. 

• Correlation The correlator integration time should be set to 1 s to allow 
good tracking of the visibility phases with time. 

6.6.2 Data Reduction 

• Loading data A sufficiently small interval between the calibration table 
entries should be chosen, 5 s to 8 s seemed appropriate to us. 

• Flagging Flagging of the data that were recorded during a frequency 
change is required. The times that need to be flagged at the beginning 
of a half-cycle with QUACK can be estimated using plots of visibility 
ampUtude with time on caUbrator scans. 

• Standard steps The following steps do not differ from standard VLBI 
data reduction: amplitude calibration using Tsys and gain measurements 
and calibration of sampler threshold errors, removal of parallactic angles 
from the visibility phases, bandpass cahbration and manual calibration 
of the inter-IF phase offsets using self-calibration on a strong source. 

• Ionospheric errors The ionospheric delay needs to be caUbrated care- 
fully. GPS -based models applied with TECOR were not accurate 
enough. If the source has been observed at L-band, the dispersive delay 
should be estimated and applied to all frequencies. We have no experi- 
ence with this procedure and cannot give advice on details. 

• Self-calibration Self-calibration should be carried out at Vref and will 
yield one phase, delay and rate solution per half-cycle if the solution 
interval exceeds the half-cycle time. 

• Phase scaUng The solution (SN) table from self-caUbration is written 
to a text file using TBOUT. The table is read and processed outside 
AlPS with FFSTG, the fast frequency switching table generator, and 
read into AlPS with TBIN. The scaled-up phase solutions are applied 
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to the most recent Vt calibration table. The data can now be averaged 
in time and frequency to increase the SNR and for scientific interpre- 
tation. Contact Enno Middelberg (emiddelb@mpifr-bonn.mpg.de) for 
help installing Python and using FFSTG. 
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Figure 6.2: Illustration of the 
flagging scheme. The period 
03:08-03:13 UT shows 3C279 
data, the period 03:13-03:23 

UT shows 3C 273 data. Occa- 
sional outliers are due to vary- 
ing subreflector rotation times. 
Top: 1 5 GHz visibility am- 
plitudes on live baselines to 
the reference antenna, LA. No 
flagging has been applied. Bot- 
tom: The same visibility am- 
plitudes after flagging. 
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Figure 6.3: Phases of the 
same visibilities shown in 
Fig. 6.2. Top: Phases before 
iiagging. Bottom: Phases after 
flagging. 



200 -BR-LA'(1 -5J ' ' - 

100 - * ^ * ^ I * * * * * * 

-100 - « ;^ — -«^>»,.«*»».«V.»_ 

\ H h \ H * I ' 

200 -l<P-LA(4-5) ^ ^ ^ ^ ^ ^ ' I 

100-} * * f * * * : ^ I* * * t 
_ ^ * V V *- V : . * ^ ; . 

-100 X J' . 't. ^ ^i. ^ ^ ^ 3. 

200 -LA-Nl'(5-7) ^ ^ ^ ^ ' ^ 

i 100-* * * * ^ * r * * V - * - ?■ 

? - J- - ^ r , . .* . 

■■"W - , * * , ^ • - 

200 -LA-Ov'(5-8) ^ ^ ^ ^ ^ | ~ 

100 - f » ^ **" ^ * * t 

** * * * * * * t .Z 

•lOO - 1 * t , * 

200 -LA-Pt'(5-9) ^ ^ ^ f I 

I""-* .* * * - *- - > * . *- i« , - ' 

0-,^ , ^ ^ f f . , t * , * , - 

-100 - i * * t I * * * * *. 

lI t ij 1 \ t \ ^ L ^ I 

1/0310 12 14 16 18 20 22 

TIME (HOURS) 

200^BR-LA'(1-5) ' ' ' ' ' ' - 

100 * 

- - ^ . 

.100- ••»«.---♦--. 

200 KP-LA(4-5) ^ '' ^ ^ ^ ^ ' ' 

100 * 

- - - * » - • 

-100- 

200 -LA-Nl'(5-7) ^ ^ ^ ^ ^ ' ^ 

Sloo- * 

lo---'****" 
a 

■100 - 

200 -LA-Ov'(5-8) ^ ^ ^ ^ ^ ' - 

100 - - - ' * " ** 
0- * 
-100 - " * 

200 -LA-Pt'(§-9) ^ ^ ^ ^ ^ ' - 
100- 

, . ' - ' 

-100 - 

^ \ \ \ \ \ I 

1/0310 12 14 16 18 20 22 
TIME (HOURS) 



158 



6.7. Figures 



100 



100 




•*•- 

-i '-r — 45 — 


++^ 


7 ' 


^ 




100 


100 




KP ( 1 ■ 4J tV* 

•tn * 
H '-^ — ^ — 


^ i 
'J't 








— 4^ 


100 


100 




l-A ( 1 ' 5 ) 


it 




^ * 1 


:v - 


100 


100 




JjL ( 1 - 7 ) 








— ^ 


100 


100 


-Bp 


i)V(1,#) J* 


*t 


ft 


+ f + 
+ ^ 


— ^ 

+++ 

— ?f- 


100 



100 




i>T(1-'9(*) *.'■• 








— ^ 


100 


100 


??. 


ic ( 1 ;.i*o ) 






1" 1 






1/00 1/02 


1/04 1/06 
TIME (HOURS) 


1/08 


1/10 


100 


100 


BR 


'HN(ftiS) if: 










100 


100 




KP(1-4) 

V 


\ 

— — 


% 




— 
— 


100 



100 




LA(1*-5) >^ 

♦} 

+4 * 




^» 


% 


^ 


100 


100 


in 












100 


100 


_BR 


ov(i-J) Jf* 

* 


> 




% 

** 


H- 

1^ 


100 


100 


_BR 


J 


\ 








100 


100 




'sc(iy>o) J.^ 


*> 1 


+^ 1 


1 





TIME (HOURS) 



Figure 6.4: Raw and cal- 
ibrated visibility phases on 
baselines to Brewster. Top: 
NGC4261 raw 43 GHz visi- 
bility phases on baselines to 
Brewster with only delay cal- 
ibration applied. Bottom: 
NGC4261 calibrated 43 GHz 
visibility phases on baselines 
to Brewster Calibration used 
scaled-up phase solutions from 
fringe-fitting with a clean com- 
ponent model at 15 GHz. 
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Figure 6.5: Expanded sec- 
tions of Fig. 6.4 
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Figure 6.7: Expanded sec- 
tions of Fig. 6.6 
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Figure 6.8: Structure functions from 25 min of observation of NGC4261 at 43 GHz (the structure func- 
tions have been cut off after 1150 s due to a lack of points at longer lags). Solid black lines show raw data 
observed at night with delay calibration only, solid red lines show the same data calibrated with scaled-up 
phase solutions from fringe-fitting at 15 GHz, and solid blue lines show the same data without the residual 
phase rate. Dashed lines show the same stages of calibration, respectively, but from a scan observed during 
the late evening, where ionospheric effects increased the phase noise. The horizontal dashed black lines 
indicate the theoretical phase noise. The median residual phase noise of calibrated data, with residual phase 
rates, is 33° at night and 50° during sunset, and without residual phase rates is 31° at night and 44° during 
sunset. 
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Figure 6.9: The effect 
of ionospheric corrections. 
NGC4261 calibrated 43 GHz 
visibility phases on baselines 
to Los Alamos as in Fig. 6.6 
(lower panel), but no iono- 
spheric correction has been 
applied. The residual phase 
rates are much higher. 



100 


-100 

100 


-100 

100 


-100 

! 100 

h 
' -100 

100 


-100 

100 


-100 



Jitj - LA ( 3 ^5 ) 



BB-LA(1 - 3j 



KP-LA(4-5) 



1., 



TuaTnlTT^'SI ^ 

4 



LA-OV(5-,8) 



LA,-F^T(5-9) ^ 



1/00 1/02 



ifH - 



H ^ 



1/04 1/06 1/08 

TIME (HOURS) 




Figure 6.10: The effect 
of using a model in self- 
calibration. NGC4261 
calibrated 43 GHz visibility 
phases on baselines to Los 
Alamos. Calibration has been 
done with scaled-up phase 
solutions from fringe-fitting at 
15 GHz using a point source 
model (light blue) and a clean 
component model (green), 
showing the effect of source 
structure at 15 GHz in the 
calibration at 43 GHz. 
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Figure 6.1 1 : Coherence improvement at 43 GHz. NGC 4261 was observed in the early evening (centred 
on 01:06 UT, left column) and at night (centred on 03:40 UT, right column) at 43 GHz with interleaved 
15 GHz observations for calibration. 43 GHz data with (solid lines) and without (dashed Unes) phase- 
calibration from 15 GHz have been averaged in frequency and in time over each 43 GHz half-cycle. The 
normaUzed vector svmi has been computed on all baselines to a given station over progressively larger time 
intervals and assuming unit amplitude. Bold lines show the average over aU basehnes. 
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Residual LL map. Array: BFHKLMNOPS 
NGC4261 at 43.101 GHz 2003 May 05 




4 2 -2 
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Contours %: -16 16 32 64 



Figure 6.1 2: Naturally weighted, full-resolution image of NGC 4261 at 43 GHz, calibrated with scaled-up 
phase solutions from 15 GHz. FRING has been used to solve for one residual phase and rate solution per 
25 min scan before exporting the data to Drfmap. No further self-calibration has been applied. The image 
noise is 4.4 mjy beam"' and the dynamic range is 18:1. 
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Clean LL map. Array: BFHKLMNOPS 
NGC4261 at 43,101 GHz 200,3 May 05 
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Figure 6.13: Data and imaging parameters as in Fig. 6.12, but several cycles of phase self-calibration with 
a solution interval of 30 s and one cycle of amplitude self-calibration with a solution interval of 12 h have 
been applied. The image noise is 0.78 mjy beam" ' and the dynamic range is 122:1. 
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Figure 6.14: 86 GHz visibil- 
ity phases. Top: NGC4261 
raw 86 GHz visibility phases 
on baselines among Fort 
Davis, Kitt Peak, Los Alamos, 
Owens Valley and Pie Town 
with only delay calibration 
applied. Bottom: 86 GHz 
visibility phases on baselines 
among Fort Davis, Kitt Peak, 
Los Alamos, Owens Valley 
and Pie Town. Calibration 
has been done with scaled-up 
phase solutions from fringe- 
fitting at 15 GHz using a clean 
component model. Good 
detections were made during 
almost every 25 min scan 
observed at night between 
2:00 h UT and 7:00 h UT, 
when the ionosphere was 
stable. 
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Figure 6.15: Structure functions of a 25 min scan at 86 GHz. Colours show the same stages of calibration 
as is Fig. 6.8. 
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Figure 6.16: The coherence improvement at 86 GHz with fast frequency switching on some of the base- 
lines shown in Fig. 6.14. Errors have not been computed due to the reduced number of data points. 
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6.7. Figures 



Residual LL map. Array: BFHKLMNOPS 
NGC4261 at 86.226 GHz 2003 May 06 
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Clean LL map. Array: BFHKtMNOPS 
NGC4261 at 86.226 GHz 2003 May 06 
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Figure 6.17: 86GHz images 
ofNGC4261. Top: Naturally 
weighted, full-resolution im- 
age of NGC4261 at 86 GHz, 
calibrated with scaled-up 
phase solutions from 15 GHz. 
FRING has been used to 
solve for one residual phase 
and rate solution per 25 min 
scan before exporting the 
data to Difmap. No further 
self-caUbration has been 
applied. The image noise 
is 8.4mJybeam^' and the 
dynamic range is 7:1. Bottom: 
Data and imaging parameters 
as above, but several cycles 
of phase self-calibration have 
been applied. The image noise 
is 7.4mjybeam^' and the 
dynamic range is 12:1. 
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Conclusions 



In an attempt to measure the magnetic field strength in the sub-pc scale 
environment of AGN through Faraday rotation, we have used the VLBA at 
15.4 GHz to look for polarized emission in NGC 3079, NGC 1052, NGC 4261, 
Hydra A, Centaurus A and Cygnus A. We selected these nearby (D < 200 Mpc) 
objects for high hnear resolution and good evidence for free-free absorption 
in a foreground medium. If the medium is interspersed with a magnetic 
field, Faraday rotation is expected to occur. Surprisingly, we did not detect 
any polarized emission in five out of six sources, and we found one source 
(Cygnus A) to be polarized at a level of only (0.44 zt 0.3) %. The results for 
one source (Centaurus A) were inconclusive due to poor (m,v) coverage. The 
non-detection of polarized emission is not unexpected in regions where a 
dense free-free absorber is seen, but is unexpected for regions that apparently 
are unabsorbed. We have supplemented our observations of NGC 3079 
with multi-epoch, multi-frequency VLB I observations to study spectral prop- 
erties and proper motions. From our study, we draw the following conclusions. 

1 . The non-detection of polarized emission is significantly different from 
VLBI surveys targeted at highly beamed quasars. We find that our selection 
criteria biased the sample towards unbeamed, double-sided radio jets in which 
circumnuclear absorbers are seen edge-on. 

2. We expect that the radio emission from the sources is synchrotron 
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emission which is intrinsically polarized, and conclude that the radiation 
is depolarized in a foreground screen. The rotation measure in the screen 
must be variable on scales much smaller than the synthesized beam size so 
that differential Faraday rotation causes a net decrease of polarized emission 
when averaged across the beam area. In our sources, such screens are 
detected through their free-free absorption only in very small regions across 
parts of the radio jets where they are interpreted as free-free absorption in 
circumnuclear tori, but their opacity rapidly decreases away from that region 
and they are not detected along the jets. This imposes constraints on the 
physical parameters in the foreground screen. The screen must be dense 
enough to depolarize at the relatively high frequency of 15.4 GHz and at the 
same time, must be transparent to radio emission along most of the jet. We 
find that the electron density in the screen must be of the order of 25cm~^ 
for magnetic fields of the order of 0. 1 mG and for path lengths of a few parsecs. 

3. Considering possible absorbers in front of the jets in which the Faraday 
depolarization occurs, we present a model in which ionized material is 
accreted in a spherical flow from scales of > lOpc to sub-pc scales (Bondi 
accretion). Assuming equipartition magnetic fields, the model predicts the 
necessary Faraday rotation and, at the same time, is transparent to radio 
emission at 15.4 GHz. The accretion rates are probably signifcant lower in 
the Seyferts NGC 3079 and NGC 1052 than in the other four objects, and the 
rotation measures predicted by the model are too low. However, we cannot 
rule out that the depolarization in these two objects is due to a combination of 
the ionized fraction of the interstellar gas, H II regions and NLR clouds. 

4. Our observations of NGC 3079 have yielded the following additional 

results. In two radio components, A and E, we find good evidence for a fore- 
ground free-free absorber, with column densities of the order of lO^^cm"^. 
The intrinsic (unabsorbed) spectrum of E is extremely steep, with a = — 4.1, 
the best explanation of which is by a non-power-law distribution of electron en- 
ergies. The spatial distribution of peak frequencies found in four components 
in NGC 3079 is suggestive of a foreground absorber with a density gradient. 

We find proper motions of the order of 0.1c among the three strongest 
components. A, B and E. The lack of absolute positions for all but the last 
epoch did not allow us to define a kinematic centre in NGC 3079, but we prefer 
a scenario in which either B is stationary (possibly equvalent to the case in 
which the kinematic centre is not coincident with any of the components) or 
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the components move on helical trajectories. If B is stationary, component E 
is moving towards the AGN, which is unexpected. 

In a statistical analysis of VLBl observations of Seyfert galaxies, we find 
that the spectral and structural properties of NGC 3079 are not uncommon, 
though the spectra are at the positive end of the distribution. 

5. For higher resolution observations of nearby and hence mostly weak 
AGN, we have developed the phase-referencing strategy of fast frequency 
switching, in which the source is observed continuously while rapidly switch- 
ing between the target frequency and a lower reference frequency. Phase solu- 
tions from self-calibration at the low frequency are multiplied by the frequency 
ratio and are used to cahbrate the high-frequency visibilities. We describe the 
strategy and results of an observation of NGC 4261. We found that short- 
term phase fluctuations at 43 GHz and 86 GHz can be calibrated with residual 
rms phase errors of 33° and 70°, respectively using interleaved 15 GHz scans. 
Ionospheric phase fluctuations on timescales of minutes to hours caused resid- 
ual phase drifts of up to a turn of phase and did not allow us to fully calibrate 
the data without self-calibration. We present the first detection of NGC 4261 
at 86 GHz which was not possible with conventional VLBI because it is weak. 
To our knowledge, this is so far the weakest source detected with VLBI at that 
frequency. 
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